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1»  ABSTRACT 


Failure  of  microelectronic  circuits  caused  by  electrical  overstresa  has  been 
investigated  theoretically  and  experimentally .  Computer  calculations  on  heat  flov 
in  silicon  structures  have  pointed  out  the  necessity  of  using  temperature  dependent 
thermal  constats  for  silicon.  The  effect  of  surface  layers  and  pover  distribution 
on  the  peak  temperature  in  silicon  devices  under  transient  conditions  has  been 
investigated. 

Integrated  circuits  representing  several  fabrication  technologies  from 
several  different  manufacturers  were  pulsed  to  failure  over  pulse  lengths  from 
100  nsec  to  10  msec.  The  threshold  pover  per  unit  area  to  produce  permanent 
damage  was  found  to  vary  by  a  factor  of  20  among  the  parts  tested.  Peak  tempera¬ 
tures  to  initiate  current  constriction  vere  estimated  to  be  1000-1200°C. 
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ABSTRACT 


Failure  of  microelectronic  circuit*  caused  by  electrical  over¬ 
stress  has  been  investigated  theoretically  and  experimentally .  Computer 
calculations  on  heat  flow  in  silicon  struc  ores  have  pointed  out  the 
necessity  of  using  temperature  dependent  thermal  constants  for  silicon - 
The  effect  of  surface  layers  and  power  distribution  on  the  peak  temperature 
in  silicon  devices  under  transient  conditions  has  been  investigated. 

Integrated  circuits  representing  several  fabrication  technologies 
from  several  different  manufacturers  were  pulsed  to  failure  over  pulse 
lengths  from  100  nsec  to  10  msec.  The  threshold  power  per  unit  area  to 
produce  permanent  damage  was  found  to  vary  by  a  factor  of  20  among  the 
parts  tested.  Peak  temperatures  to  initiate  current  constriction  were 
estimated  to  be  1000-1200°C. 


EVALUATION 


The  objective  of  this  effort  vas  to  provide  a  physical  basis  for 
assessing  the  reliability  of  microelectronic  devices  when  subjected 
to  electrical  overstress.  The  contractor  performed  a  combined 
analytical  and  experimental  program  to  define  the  significant  parameters 
associated  witn  overstress  failure  of  microelectronic  structures.  While 
considerable  overstress  testing  has  been  done  on  microelectronic 
structures,  this  effort  adds  substantial  understanding  to  the  detailed 
mechanisms  of  failure.  This  was  done  by  the  use  of  well  controlled 
experimental  procedures  and  the  thorough  analysis  of  failed  devices 
and  associated  electrical  waveforms.  Dr.  H.  Domingos  was  able  to 
relate  his  results  on  emitter  base  flashover  shorts  in  integrated 
crcuits  to  recent  work  on  diode  second  breakdown  in  thin  silicon 
filu-  on  sapphire.  This  provides  a  good  physical  understanding 
of  the  surface  flashover  type  of  breakdown  in  microelectronic 
structures.  In  addition,  he  demonstrated  that  differences  in  ability 
to  withstand  overstress  in  the  circuits  tested  were  related  to  the 
distance  from  the  P-15  /junction  to  the  contact  on  the  high  resistivity 
side  of  the  junction  &nd  the  existence  of  multiple  conduction  paths 
in  the  microelectronic  device. 

The  information  obtained  from  this  effort  will  be  used  in  the 
interpretation  of  failures  due  to  overstress  in  microelectronic 
systems.  Additional  in-house  effort  in  the  electrical  overstress 
area  is  being  considered  using  test  structures  consistent  with  LSI 
processing  and  LSI  circuits. 

JOHN  F.  CARROLL 
Reliability  Physics  Section 
Reliability  Branch 
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It  is  a  veil  known  fact  that  voltage,  current,  and  power  ratings  of 
semiconductor  devices  cast  be  scrupulously  observed  in  any  application.  Un¬ 
like  vacuus  tubes,  where  operation  above  rated  values  merely  reduces  the 
operating  lifetime,  overstress  of  semiconductor  devices  often  leads  to  catas¬ 
trophic  failure.  Because  electronic  circuits  must  function  in  environments 
where  voltage  and  current  surges  due  to  power  supply  transients,  line  dis¬ 
turbances,  electromagnetic  pulses,  etc.,  may  be  present,  a  clear  understand¬ 
ing  of  failure  due  to  such  conditions  is  needed.  Apart  from  the  uneasiness 
which  lack  of  understanding  generates,  it  can  also  result  in  systems  which 
are  either  unreliable  or  needlessly  overdesigned.  Better  understanding 
leads  to  better  component  design  as  well  as  more  efficient  engineering  and 
more  reliable  equipment. 

The  type  of  failure  investigated  in  this  study  is  degradation  due 
to  extreme  temperature  rises  resulting  from  excessive  power  dissipation  as 
distinct  from  degradation  due  to  such  mechanisms  as  aluminum  migration, 
purple  plague,  chemical  and  metallurgical  reactions,  etc.  Indeed,  in  any 
device,  semiconductor  or  otherwise,  excessive  power  dissipation  is  of  conse¬ 
quence  only  because  it  raises  tne  temperature:  as  long  as  the  device  can  be 
kept  cool  the  power  dissipation  is  a  secondary  consideration.  In  semicon¬ 
ductor  devices  the  temperature  rise  eventually  leads  to  melting  of  the  sili¬ 
con  or  other  materials  comprising  the  circuit. 

The  earliest  efforts  to  systematically  study  device  failure  were 
undertaken  in  the  late  fifties.  The  recognition  of  the  phenomenon  known  as 
second  breakdown  caused  a  flurry  of  activity  in  the  mid-sixties.  This  has 
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not  yet  abated  as  the  development  of  high  pulse  power  microwave  diodes  has 
entered  the  scene. 

Work  done  at  other  laboratories  with  silicon-on-sapphire  diodes 
while  this  program  was  in  progress  has  led  to  a  much  clearer  understanding  of 
second  breakdown  in  semiconductor  devices,  particularly  single  junctions. 

The  results  of  the  present  study  tend  to  agree  with  these  findings  and  thus 
it  can  be  said  that  a  good  picture  of  second  breakdown  in  diodes  is  now 
available . 

This  report  is  the  final  report  on  a  program  of  theoretical  and 
experimental  investigations  on  breakdown  in  m  c  ^electronic  circuits .  First 
a  study  of  heat  flow  in  silicon  is  presented  If  the  spatial  distribution  of 
the  dissipated  power  is  known,  a  knowledge  of  the  thermal  properties  of  the 
materials  enables  the  temperature  distribution  to  be  calculated.  Unfor¬ 
tunately  the  power  dissipation  depends  on  the  characteristics  of  the  source 
as  well  as  the  electrical  characteristics  of  the  circuit,  which  in  turn 
depends  on  the  temperature.  This  leads  to  a  complex  situation.  Finite- 
difference  jolutions  have  been  obtained  for  temperature  distributions  in 
idealized  one-dimensional  situations  taking  into  account  the  variation  of 
thermal  parameters  with  temperature, 

A  section  on  review  of  second  breakdown  theories  is  followed  by  a 
presentation  of  the  experimental  results.  Several  t^pes  of  microelectronic 
circuits  were  pulsed  into  second  breakdown  to  determine  threshold  energies 
for  damage  as  a  function  of  pulse  width  and  to  determine  the  nature  of  the 
damage,  The  results  are  discussed  in  terms  of  current  models  of  second 
bre&Kdown  and  the  heat  flow  calculations. 


2 


-  'j  '  -  .-  T  ,  -  •  *  '  •  ^ 


SECTION  II 

HEAT  TRANSFER  IN  INTEGRATED  CIRCUITS 


A:  General  Heat  Flow 

When  power  is  applied  to  a  silicon  monolithic  integrated  circuit, 
most  of  the  energy  is  dissipated  in  the  circuit  elements  within  the  chip. 

The  heat  can  be  carried  away  by  conduction,  convection,  and  radiation  to 
other  parts  of  the  chip  and  ultimately  to  the  package  and  to  the  ambient. 
Nearly  all  of  the  heat  flow  is  through  the  silicon  to  the  header  and  hence 
to  the  rest  of  the  package.  Heat  can  also  flow  along  metallization  runs  to 
the  leads  and  then  to  the  package ;  through  the  oxide  and  the  ambient  gas  in 
the  package;  or  be  radiated  to  cooler  parts  of  the  structure.  Heat  transfer 
by  all  three  processes,  radiation,  convection,  and  conduction  was  investiga¬ 
ted;  only  transfer  by  conduction  was  found  to  be  significant  in  typical  micro- 
circuits  . 

Convective  heat  transfer  is  sometimes  interpreted  to  also  include 
conduction  by  the  convecting  medium;  however,  it  differs  from  pure  conduction 
in  that  the  transmission  of  energy  is  augmented  by  movement  of  the  fluid 
itself.  Free  convection  comes  about  as  a  result  of  the  change  in  buoyancy 
of  the  fluid  as  a  result  of  uneven  heating.  If  the  fluid  is  confined  to  a 
very  small  space  the  buoyant  forces  are  not  sufficient  to  overcome  viscosity, 
and  fluid  motion,  hence  convection,  does  not  occur.  This  is  the  case  for 
air  inside  an  integrated  circuit  package. 

Radiation  is  electromagnetic  in  nature  and  radiative  exchanges  are 
governed  by  the  laws  of  electromagnetic  theory  or  optics.  Suppose  we  con¬ 
sider  radiation  from  the  hot  surface  of  a  chip  to  the  gold-plated  lid  of  the 
package.  Treating  the  two  surfaces  as  specularly  reflecting  plates  at  1000^ 
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and  290°K,  respectively,  a  reasonable  value  for  the  heat  transfer  coefficient 
is  0.2  watts /era2.1  For  a  chip  100  mils  on  a  side  radiation  would  account  for 
12.9  mw.  This  is  a  liberal  estimate  since  not  all  of  the  chip  surface  would 
be  at  1000°K,  and  the  thermal  oxide  or  deposited  glass  is  a  good  insulator. 
Even  so,  this  power  level  is  about  three  orders  of  magnitude  less  than  power 
levels  of  concern  in  pulse  testing. 

Radiation  from  the  surface  down  through  the  silicon  to  the  header 
is  of  even  less  importance.  In  the  first  place  silicon  has  a  fundamental 
absorption  edge  at  1:13  ym.  Radiation  at  shorter  wavelengths  would  be  sharp¬ 
ly  attenuated  by  the  excitation  of  electrons  to  the  conduction  band.  The 
dominant  wavelength  of  a  black  body  radiator  at  1000°K  is  2.9  ym.  Radiation 
at  this  wavelength  will  be  attenuated  by  joule  heating  of  the  lossy  silicon. 
For  1  ohm~cm  material  the  penetration  depth  is  only  5  ym. 


B.  Heat  Flow  by  Conduction 

The  differential  equation  which  describes  the  conduction  of  heat  in 
a  solid  is  basically  a  continuity  equation  for  the  conservation  of  heat  energy 
and  can  be  written  in  its  most  general  form  as 


pc 


9T 

3t 


■  * +  ^ + 1  '<I>2+f)2+(f>2) 


‘  9y 


(1) 


where  p  is  the  density  in  grams/cm3,  c  is  the  specific  heat  in  Joules/ 
gram-°C,  T  is  the  temperature  in  °C,  t  is  in  seconds,  q  is  the  power  density 
in  watts/cm3,  k  is  the  thermal  conductivity  in  watts/cm-°C,  and  x,  y,  and  z 
are  space  coordinates  in  cm.  q  can  be  a  function  of  time  and  space  and  k 
can  be  a  function  of  temperature. 

As  written,  Equation  (l)  is  non-linear  and  can  be  solved  only  in 
certain  special  cases.  If  k  is  not  a  function  of  temperature,  the  equation 


3 


becomes  a  linear  partial  differential  equation  for  which  many  solutions  are 
tabulated.2 

Solutions  to  the  heat  flow  equation  must  satisfy  the  initial  con¬ 
dition  that  at  t  =  0 


T  =  f(x,y,z) 


and  one  or  more  boundary  conditions 

i)  temperature  at  one  or  more  surfaces  is  constrained 


Tsurf  =  S(x,y,z,t) 


ii)  a  surface  is  insulated  so  that  no  heat  flows  across  it 

3T 


3n 


=  0 


iii)  miscellaneous  boundary  conditions 

9T 


k  ~  +  H(T-T  )  -  F  =  0  . 
3n  o 


The  first  term  is  heat  flow  across  the  surface  by 
conduction,  the  second  term  is  the  "radiation" 
boundary  condition  and  represents  a  heat  loss  pro¬ 
portional  to  the  difference  between  the  surface  at 


temperature  T  and  another  surface  at  Tq ,  and  the 


last  term  is  the  heat  flux  into  the  surface  from 
an  external  source. 

A  number  of  simple  solutions  to  heat  flc>r  problems  will  now  be 
applied  directly  to  special  situations  oceuring  in  microcircuits.  The  follow¬ 
ing  constants  for  silicon  will  be  used  in  the  calculations  tnat  follow: 

p  ~  2.328  gm/cm3 


c  =  0.7  J/gm-°C 

k  ~  1.1»5  W/cm-°C 

D  =  0.9  cm2 /sec 
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1.  Adiabatic  heat  problem 

Consider  a  small  volume  initially  at  a  uniform  temperature  Tq  to 
which  a  uniform  power  density  Po  watts/cm3  is  supplied  throughout,  with  insu¬ 
lated  boundaries .  The  differential  equation  becomes 

P 


with  the  solution 


8T  _  _o 
3t  pc 


P  t 

T  -  T  =  —  =  0.6ll»  P  t  °C 
o  pc  o 


(2) 


(3) 


where  t  is  in  seconds  and  Pq  is  in  watts/cm3.  This  solution  is  applicable  to 
thermal,  problems  during  a  time  interval  short  enough  that  flow  by  diffusion 
is  negligible,  as,  for  instance,  estimating  the  temperature  rise  in  a  junc¬ 
tion  during  the  first  few  nanoseconds  after  power  is  applied.  The  equation 
gives  an  upper  limit  to  trie  temperature  rise  and  a  reasonable  accurate 
answer  as  long  as  vUtT  is  less  than  the  dimensions  of  the  region. 

2.  Constant  power  into  the  surface  of  a.  semi-infinite  medium 

This  is  a  one-dimensional  problem  where,  the  silicon  extends  from 
x  =  0  to  infinity,  with  a  power  density  Fo  in  watts/cm2  at  the  surface. 


The  differential  equation  is 


3T  _  k  c»2T  _  32T 


3t  pc  3X2 


Dx" 


(U) 


The  solution  is 


6 


^  ?  ■ 


■.\ri  ;!V  '* 


if  ''p'-tr 


T-T 

o 


,  /Dt^  -x2/UDt 
1 l-~)  e 


-  |  erfc 


x 

2(Dt)?5 


} 


(5) 


This  solution  is  applicable  to  the  case  of  uniform  power  dissipation  in  a 
comparatively  thin  layer  at  the  surface  of  a  chip.  The  solution  is  valid  as 
long  as  lateral  heat  flow  and  diffusion  to  the  header  can  be  neglected. 
Equation  (5)  is  plotted  in  normalized  form  in  Figure  1  and  again  in  Figure  2 
as  a  function  of  time  for  x  =  0.  For  a  total  power  of  25  watts  into  an  area 
10  mils  x  10  mils,  F^  =  Uo  x  102  watts/cm2  and  the  maximum  abscissa  in 
Figure  1  is  a  temperature  rise  of  l480°C. 

3.  Constant  power  into  the  surface  of  a  semiconductor  of  finite  thickness, 
with  an  ideal  heat  sink  ’  “ 


This  case  applies  when  the  finite  thickness  of  the  die  must  be  con¬ 
sidered.  It  is  assumed  that  edge  effects  are  still  negligible,  and  that  the 
header  is  able  to  maintain  the  x  =  L  plane  at  constant  temperature  T  .  The 
solution  is 


T-T  = 
o 


F0(l 


-x) 


8F  L 
o 


_  ^  (-x)n 

kTT2  n=0  (2n+i)2 


D(2n+l)~n2t 


L2  .  (2n+l)ir(L-x)  ,  r\ 

Sin  "T2  ( o ) 

k'CIj 


This  is  plotted  in  Figure  3.  Figure  t  shows  the 
function  of  time  at  x  =  0.  For  25  watts  into  an 


temperature  rise  as  a 
area  10  mils  x  10  mils 


=  500  /isec 


Temperature  rise  at  the  surface  of  a  serai-infinite 
solid  due  to  a  power  density  F  watts/cra2  into  the 
surface. 


cm  /  wcitt ) 


x/L  I 

Fig.  3.  Temperature  distribution  due  to  a  constant  power  del 
the  surface  of  a  slab  terminated  at  x  =  L  with  a  pe*' 


t/L2  (sec/cm2) 


Fig.  Temperature  rise  at  the  surface  of  a  slab  terminated 
at  x  =  L  with  a  perfect  heat  sink..  The  asymptotic 
value  is  0.69. 
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and  vitb  L  =  10  ails,  multiply  the  ordinate  by  1000-  t/Ld  =  0.1  then  corre¬ 
sponds  to  62.5  usee.  The  steady  state  temperature  rise  in  this  case  is 
690°C,  which  can  be  verified  by  calculating  the  thermal  resistance.  For  tee 
general  case,  the  presence  of  the  heat  sink  is  important  only  for  t/L2>0.2, 
.and  steady-state  is  reached  for  t/L2>10. 

4.  Semi-infinite  region,  with  power  density  P  vatts/cs-  into  a  finite 
thickness  at  the  surf  ac e  ° 


The  solution  is 
P  *t 

T-T  =  -2-  {i_2  l2erfc  -  2  i2erfc  ~=z  3fcr  0<x<l  (7) 

°  oc  2/Dt  2»/Dt 

2P  t 

=  — —  {  x2erfc  —  pL:  -  i2erfc  x~  ^L:  )  for  x>L 
pc  2m  2/Dt 

where  i2erfc  u  =  h  {(l+2y2)erfc  y  — p  ye  M  } 

/it 

This  solution  is  useful  when  edge  effects  and  finite  length  effects  can  he 
neglected,  but  when  the  heat-dissipating  source  at  the  surface  must  be  con¬ 
sidered  in  more  detail.  The  thickness  of  epitax5,al  layers  and  depths  of 
diffusion  for  example  are  often  important  parameters.  This  model  is  also  the 
basis  of  numerical  calculations  presented  later  in  this  section. 

Equation  (7)  is  plotted  in  Figures  5  and  6.  For  comparison,  we  can 
again  take  25  watts  into  an  area  10  mils  x  in  mils  with  L  =  10  ym.  In  this 
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case  the  ordinate  in  the  two  figures  xiiO  gives  the  temperature  rise,  and  the 
parameter  in  Figure  5  is  read  directly  in  ysec. 

These  figures  agree  with  the  more  general  results  of  Mosekilde, 
et.al.,3vho  have  solved  the  same  problem  but  for  a  semiconductor  of  finite 
thickness  with  an  ideal  heat  sink. 

5 .  Power  into  a  rectangular  area  on  the  surface  of  a  semi-infinite  solid. 


This  represents  the  case  of  power  dissipation  by  an  active  element 
in  a  thin  layer  at  the  surface  of  a  chip.  The  assumptions  again  are  that  the 
edges  of  the  die  are  far  removed  from  the  active  element  and  that  the  die- 
header  interface  is  far  enough  away  during  times  of  interest. 

The  solution  tak«s  into  account  the  flow  of  heat  laterally  as  well 
as  vertically  down  into  the  silicon.  The  mathematics  becomes  complicated, 
but  a  simple  expression  can  be  found  for  the  maximum  temperature  rise,  T^, 
which  occurs  at  the  surface  at  the  center  of  the  active  element.  This  ex¬ 
pression  is 

T  -T  =  ~r~-  {a  sinh-1  ~  +  l  sinh"1  (8) 

m  o  irkaS/  a  & 

For  the  case  a  =  l 

■1  -T  =  3.87  x  10’  (c,) 

m  o  a(ym) 

Steady  state  temperature  distributions  in  a  finite  chip  with  the 
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heat  source  at  the  center  of  the  top  surface  and  an  ideal  heat  sink  at  the 
hottem  surface  have  been  calculated  by  Lindsted  and  Murty , 4  Hor.salixed 
curves  have  been  plotted  for  the  maximum  temperature  rise  for  s  range  of 
geometries . 

6.  Power  into  a  circular  area  on  the  surface  of  a  semi-infinite  solid. 


This  is  similar  to  case  5;  however,  for  this  geometry  it  is  rela¬ 
tively  simple  to  determine  T(z,t)  along  the  axis  of  symmetry  r  =  0.  This  not 
only  gives  the  maximum  temperature  rise,  but  allows  one  to  examine  the 
temperature  distribution  as  a  function  of  time  and  make  a  judgement  as  tc 
when  die  thickness  and  edge  effects  are  no  longer  negligible. 

The  solution  is 


T-T 

o 


2/Dt  P 
irkr  2 


o 


(ierfc  -~ 

2m 


(z2+ro2)h 

ierfc  - - — — } 

2/Dt 


(10) 


1  -  v2 

where  ierfc  p  =  --  e 

YT, 


p  erfo  p  , 


This  equation  has  several  assymptotic  approximations.  When  z  =  0,  T  =  T  and 


T  _T  *  «MJo  for  ^  <  0.U 


m  o 


(11) 


‘n  t-his  equation  and  those  which  follow  t  is  in  seconds,  P  is  in  watts,  and 
spatial  dimensions  are  in  centimeters.  'This  is  the  same  as  the  solution  for 
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case  2  if  x  approaches  zero'  in  Satiation  (5)-  Ibis  means  'that  jEe-a  spices 
chip  -where  D  =  0.9  cz?/sec  the  lateral  heefc  flow  awsgr  from  the  ed^e  cf  the 
circular  region  has  no  effect  m2  if  t<Q.2r  2. 

-  -  ----  j»  -©  ~  - 

In  the  steady  state 

2-1  =  0.22  ,  (12) 
so  r 

o 

valid  for  »^t/r  >5«0. 

o 

For  large  .z  and  t  - 

T-T  %  0.11  —  for  >50.  -  {23 ) 

o  z  r 

© 

ess.  3/r  >  5.0 

°  -j 

SosstiOu  {10}  is  clotted  in  -» 7  gwg  for 

z  =  0  in  Figure  8.  If  p  —  I  watt  and  r  -  1  nil,  the  ordinate  In  Figure  7 
■  o 

multiplied  by  tOO  gives  the  temperature  rise  in  °C.  The  curve  with  parameter 
1  is  the  temperature  distrihution  at  t  =  6.25  vsec  while  the  curve  labelled 
50  corresponds  to  l6  nsec,  at  which  tine  the  distribution  is  essentially 

steady-state.  In  the  steady-state  T  is  83°C  above  ambient. 

.  - 
Steady-state  solutions  for  a  finite  chip  of  cylindrical  geometry 

are  given  by  Kennedy.5  Complete  solutions  for  a  chip  of  very  large  area  but 

finite  thickness  mounted  on  a  perfect  heat  sink  are  given  in  reference  6.  If 

rQ>3I»,  where  L  is  the  chip  thickness,  the  solution  is  approximately  that 

given  by  Equation  (6),  while  if  L'3rQ  the  solution  is  approximately  Eg.  (IC). 

The  solutions  presented  above  are  useful  for  estimating  maximum 

temperatures  and  size  effects  when  the  temperature  rise  is  not  excessive. 

Unfortunately,  trie  thermal  properties  of  silicon,  especially  the  conductivity, 

are  functions  of  temp erasure .  This  raises  serious  questions  about  the 


IT 


(IJDM/UJO  —  Qtf)  d/0j(0l-l) 


Fig.  ?.  Temperature  distribution  in  a  semi-inf inlte  body  due  to  a  eonstant 
power  applied  at  the  surface  over  a  circular  area  of  radios  r0.  T 
curve  is  the  distribution  along  the  perpendieular  axis  of  the 
circular  region. 


(HOM/ujo-Oo)  d/0i(°l-U,l) 


19 


Fig.  8.  Temperature  rice  in  a  aemi-infi-nite  eolid  duo  to  power  applied  at  the  surface  over  a 
circular  area  of  rodiuo  r0.  !l*ho  temperature  rise  lo  the  maximum  value,  at  the  center 
of  the  circular  area.  The  daiihed  line  at  the  right  '.e  the  aeymptotie  value. 


validity  of  using  solutions  to  ths  linear  heat  equation.  In  the  next  sub- 

v 

section  the  variation  of  thermal  conductivity  vith  temperature  will  be  dis¬ 
cussed,  and  some  cf  the  difficulties  associated  vith  the  use  of  an  average 
value  will  be  pointed  out. 

C.  Thermal  Conductivity  of  Silicon 

The  physics  cf  the  thermal  conductivity  of  solid3  is  well  under 
stood7*3  and  good  qualitative  and  even  quantitative  agreement  between  theory 
and  measurements  is  obtain*  for  many  materials.  Any  lack  of  quantitative 
agreement  is  due  to  mathematical  complexities  and  not  to  any  lack  of  under¬ 
standing  of  the  principles. 

Thermal  conduction  in  solids  occurs  by  means  of  phonons,  electrons, 
holes,  electron-hole  pairs,  and  photons.  Lattice  or  phonon  conduction  is  the 
only  important  mechanism  in  silicon  at  temperatures  less  than  room  tempera¬ 
ture.  The  others  become  increasingly  important  at  higher  temperatures . 

Finite  phonon  thermal  conductivity  is  attained  through  the  action 
of  scattering  mechanisms  which  keep  the  phonon  distribution  in  thermal  equi¬ 
librium  with  the  lattice  temperature.  These  mechanisms  are:  diffuse  scat¬ 
tering  from  crystal  boundaries,  phonon-phonon  scattering,  isotope  scattering, 
electron  scattering,  and  scattering  from  crystal  defects  and  imperfections. 
Scattering  establishes  a  mean  free  path  for  phonon  propagation  which  can  be 
related  to  r.he  thermal  conductivity. 

F">r  an  ideal  Maxwell-Bolt zmann  gas  the  thermal  conductivity  is 

given  by 

k  =  |  cv*  (13) 

where  c  is  the  specific  heat  at  constant  volume,  v  is  the  average  velocity, 
and  l  is  the  mean  free  path.  How  do  the  scattering  mechanisms  affect  l  and 
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and  hence  k?  At  very  low  temperatures  the  number  of  phr.aons  available  to 
scatter  other  phonons  is  small.  The  electron  (hole)  population  is  small 
since  there  is  little  thermal  generation  and  even  impurity  carriers  may  be 
frozen  out.  Isotope  effects  are  of  little  importance  because  the  phonons  are 
of  long  wavelength.  This  leaves  only  boundaries  and  defects  as  effective 
scattering  agents.  In  this  case  l  is  fixed,  and,  since  e  «  T3  at  low  temp¬ 
eratures,,  k  «*  T  and  k  approaches  zero  as  T  approaches  zero. 

As  the  temperature  is  increased,  isotopes  can  be  assumed  to  have 

increasing  importance.  If  we  restrict  ourselves  to  pure  material  for  the 

moment,  electron  effects  can  still  be  neglected.  Scattering  by  other  phonons 

0 /PT 

depends  on  the  number  of  other  phonons  available,  which  depends  on  e  , 
where  9  is  the  Debye  temperature.  At  very  high  temperatures  the  phonon 
population  varies  as  T.  c  is  nearly  constant  at  high  temperatures,  hence  k 
varies  as  l/T. 

These  interactions  will  cause  the  phonon  portion  of  the  thermal 
conductivity  to  vary  as  shown  below. 


BOUNDARY 


ISOTOPE  I  PHONON 
t 


»6/2T 


(PHONON) 


Above  room  temperature  other  carriers  besides  phonons  contribute 


to  heat  transport.  Photon  conduction  is  observed  in  pure,  large  bandgap 
semiconductors  at  high  temperatures,  for  example  GaAs  above  B00°K.  It 
#  Natural  silicon  consists  of  about  92%  atomic  weight  28,  9%  29 ,  and  3 t  ^0. 
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actually  should  be  treated  as  conduction  by  radiation,  and  previously  this 
was  shown  to  be  small  in  silicon.  Glassbrenner  and  Slack9  were  unable  to 
detect  a  photon  contribution  up  to  1580°K,  within  100  degrees  of  the  melting 
point. 

Heat  conduction  by  electrons  (holes)  occurs  by  the  same  process  as 
heat  conduction  by  electrons  in  metals.  Contrary  to  metals,  however,  where 
nearly  all  of  the  heat  is  carried  by  the  numerous  free  electrons,  only  about 
6%  of  the  thermal  conductivity  of  silicon  at  the  melting  point  is  due  to 
free  carriers. 

Electron-hole  pair  conduction  is  significant  at  high  temperatures 
because  of  thermal  generation  of  large  densities  of  electron-hole  pairs. 

These  carriers  diffuse  (ambipolar  diffusion)  from  regions  of  high  temperature 
and  concentration  to  regions  of  lower  concentration.  When  recombination 
takes  place  energy  roughly  equal  to  the  band  gap  energy  is  released.  Glass¬ 
brenner  and  Slack  have  found  ambipolar  diffusion  to  comprise  32%  of  the 
thermal  conductivity  at  the  melting  point. 

The  thermal  conductivity  of  pure  silicon  is  plotted  in  Figure  9- 
The  values  are  taken  from  Glassbrenner  and  Slack.  The  variation  is  as  ex¬ 
pected,  and  indeed  Glassbrenner  and  Slack  claim  good  agreement  with  theoreti¬ 
cal  calculations.  The  T3  behavior  at  low  temperatures  is  suppressed  by  the 
scale  on  the  abscissa,  chosen  to  better  display  the  conductiviy  at  higher 
temperatures . 

The  effects  of  doping  with  boron  or  phosphorus  have  been  measured 
by  Slack.10  He  reports  that  the  room  temperature  conductivity  is  reduced 
from  1.56  to  1.2  watts/cm-°K  for  phosphorus  doping  of  2  x  1019/cm3.  Below 
room  temperature  the  relative  change  is  even  greater.  At  the  heaviest 
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Fig.  9.  Thermal  conductivity  of  pure  silicon  as  a  function  of  absolute  temperature 


doping, 'boron  at  3  ^  1020,  k  at  room  temperature  is  reduced  to  about  0.^7 
watts/cm-°K.  He  also  reports  on  earlier  work  by  other  workers  of  boron  dop¬ 
ing  at  k  x  1016/cm3.  In  this  case,  k  was  reduced  only  at  temperatures  less 
than  80°K. 

The  extreme  variation  of  k  in  Figure  9  causes  justified  reluctance 
in  applying  solutions  to  the  linear  heat  flew  equation  even  if  an  average 
value  is  assumed.  It  has  been  common  practice  to  use  about  0.9  watts/cm-°K. 
It  is  worthwhile  to  examine  carefully  a  typical  case  to  see  the  difficulties 
which  can  be  encountered. 

In  one  such  case  an  estimate  of  the  maximum  temperature  rise  under 
certain  conditions  was  needed.  The  power  input  was  2.8  watts  into  a  thin 
layer  at  the  surface  with  area  20  pm  x  kO  pm.  The  model  in  example  5 
applies.  Using  Equation  (8)  the  first  estimate  of  the  steady-state  tempera¬ 
ture  rise  is  31+3°C.  If  we  take  Tq  =  300°K  and  use  an  average  value  of  k, 
that  is 

k  =  3g(k(300)+k(T  ))  (111) 

m 

trial  and  error  yields  Tm  =  ?00°C.  If  instead  we  use  k  at  the  average 
temperature , 

k  =  k[>g(300+T  )]  (15) 

m 

then  T  =  1050°C.  A  graphical  solution  is  shown  in  Figure  10.  Equation  (lH) 
is  plotted  as  the  upper  dashed  line.  It  is  the  average  of  k(300)  =  1.56  and 
k(Tffi)  which  is  simply  Glassbrenner  and  Slack's  data  of  Figure  9  replotted. 

The  vertical  lines  with  arrowheads  indicate  ihe  two  curves  which  are  aver¬ 
aged.  The  intersection  of  the  dotted  line  with  the  solid  curve  k  =  498/T^, 
which  is  a  plot  of  Equation  (8),  gives  the  solution  =  1*97°C. 


2h 


(watt  /  °K-  cm) 


k~  l;(300) 


k-  1/2  (k(300)+  ktT^)) 


498/T. 


,,pk(l/2{300  +  Tm)) 


k  = 


500  700  900  1100  1300 

Tm(«K) 

Fig.  10. 

A  graphical  solution  for  the  temperature  rise  in 
grated  circuit  with  two  different  approximations 
thermal  conductivity.  See  text  for  explanation. 

an  inte 
for  the 

Equation  (15)  states  that  k  should  he  found  at  a  temperature  which  is  the 
average  of  SOO°K  and  5?^-.  Shis  average  is  tee  lower  -dotted  i-Ee*  in  Figure  10. 

The  intersection  with  k  -=  k(T  )  gives  this  solution,  as  ICorC. 

8&  — 

The  disparity  between  the  two  solutions  is.  a  serious  impediatert  to 
the  use  of  solutions  to  the 'linearised  equation  in  general  cases.  Itsnerical 
solutions  mist  he  resorted  to  even  though  these  ere  asre  difficult,  subject 
to  inaccuracies,  and  not  amenable  to  the  usual  corseli ration.  The  next 
section  is  devoted  to  a  discussion  of  numerical  methods  am*  the  presentation 
of  some  results  of  a  computer  program  for  solving  heat  flow  problems . 


D.  Numerical  Solutions  to  Heat  Flow  Prcblecs 

-A  simple  approach  to  tha  solution  of  the  heat  equation  is  to  qnss- 
tize  the  space  and  time  variables  end  write  a  heat  balance  equation  fer  each 
small  volume.  Knowing  the  net  he  it  flow  into  an  incremental  volume  enables 
one  to  calculate  the  temperature  rise  durir?  the  next-  incremental  time  in¬ 
terval.  If  the  space  and  time  inerx-ments  are  small  enough  an  accurate  solu¬ 
tion  to  the  differential  equation  can  be  obtained.  The  thermal  properties 
of  the  material  can  be  adjusted  as  the  temperature  changes,  and  if  the 
temperature  does  not  change  greatly  during  one  time  step  good  approximations 
to  the  physical  problem  can  be  obtained. 

It  is  sometimes  helpful  to  visualize  heat  flov  by  its  analogy  with 
an  electrical  circuit.  The  analagous  quantities  are: 

thermal  resistance  -  electrical  resistance 

thermal  heat  capacity  -  electrical  capacitance 

heat  flow  -  current  flow 

temperature  -  voltage 

For  simplicity  e  one-dimensional  example  will  be  illustrated,  with 
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the  space  variable  x  divided  into  equal  regions.  We  will  speak  of  a  node  as 
the  point  at  the  center  of  each  region.  The  node  temperature  is  measured  in 
degrees  Kelvin  or  Celsius.  The  temperature  difference  between  adjacent  nodes 
will  result  in  a  flow  of  heat  energy,  determined  by  the  thermal  resistance 
between  the  nodes.  The  mass  and  specific  heat  of  each  region  determine  how 
much  the  temperature  will  change  as  a  result  of  net  heat  flow  during  the  next 
time  interval. 

The  electrical  analog  to  the  thermal  model  is  shown  in  Figure  11. 
The  node  temperatures  are  analagous  to  the  node  voltages  in  the  electrical 
circuit.  The  circuit  capacitors  are  connected  between  the  nodes  and  ground. 
The  reason  for  this  is  that  ground  potential  represents  the  ambient  tempera¬ 
ture,  In  the  same  way  that  we  are  really  interested  only  in  the  temperature 
rise  above  ambient,  or  the  temperature  difference  between  nodes,  so  only  the 
node  voltages  with  respect  to  some  arbitrary  reference  or  the  difference  in 
potential  between  nodes  is  of  importance.  Note  that  the  capacitors  are  not 
connected  between  nodes.  Power  dissipation  in  a  region  caused  by  i2R  heat¬ 
ing,  heatingof  a  depletion  layer,  etc.,  is  represented  by  a  current  source 
directed  into  the  node. 

Writing  an  expression  for  the  heat  flow  into  the  nth  node  in  the 
thermal  model  we  obtain 


AX 


U(I -T  ) 

- EL-20- +  lAAx 


where  q  is  the  heat  density  in  watts/cm3  being  supplied  to  the  nth  node. 
Divide  this  equation  by  the  specific  beat  mass  of  the  nth  region  and  in  a 
time  At  the  temperature  rise  will  be 
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Fig.  11.  Analogy  between  the  thermal  model  (top)  and 
the  electrical  model  (bottom). 


(16) 


AT 


y  A 

~  (T  . -2T  +T  )+qAAx 
Ax  n-1  n  n+1 


Since  the  initial  conditions  for  the  problem  are  known.  Equation 
(16)  can  be  immediately  solved  for  each  node.  This  gives  the  new  temperature 
distribution  at  At.  This  technique  is  known  as  the  "explicit"  formulation  of 
the  finite  difference  equation  for  heat  flow,  or  the  Bender-Schmidt  method. 
Equation  (16)  is  derived  in  a  more  formal  way  in  texts  on  finite  difference 
equations . 

For  problems  where  the  heat  source  is  at  one  end  of  the  structure, 
cay  at  node  1,  convergence  is  faster  if  the  temperature  rise  of  the  nth  node 
is  calculated  in  terms  of  the  new  temperature  of  node  n-1 .  As  an  example  oi 
this  approach  consider  a  diode  with  a  cross  section  area  of  10  cm  dissi¬ 
pating  50  watts  in  a  2  urn  layer  at  the  surface.  The  problem  was  solved  in 
one  dimension  only,  with  the  variation  in  thermal  conductivity  approximated 
by  a  simple  analytic  expression.  The  temperature  distribution  as  a  function 
of  time  is  shown  in  Figure  12.  For  comparison  the  dotted  curve  shows  the 
temperature  distribution  with  k  constant  and  equal  to  the  room  temperature 
value.  The  difference  is  appreciable  and  once  again  emphasizes  the  need  for 
realistic  solutions  when  large  temperature  rises  are  encountered. 

This  method  has  been  used  to  investigate  the  temperature  rise  in 
aluminum  metallization  runs  on  integrated  circuits  with  good  results.11*12 
The  technique  is  simple,  yet  very  flexible.  Not  only  can  the  conductivity 
and  specific  heat  be  functions  of  temperature,  but  the  power  sources  can  be 
functions  of  temperature  and  time.  There  is  considerable  freedom  in  choos¬ 
ing  node  sizes  and  different  materials  can  be  accommodated.  The  algorithm  is 
fast.  It  does  suffer  from  one  drawback;  however,  the  solution  is  unstable 


Fig.  12.  Temperature  distribution  in  a  silicon  diode,  solved  by  an 
explicit,  finite  difference  method.  The  dotted  line  is 
the  solution  at  1  nsec,  assuming  a  constant  value  for  the 
thermal  conductivity. 


if  the  time  increment  is  too  large.  Once  the  node  sizes  have  been  selected 


At  must  satisfy  the  condition 

«  Hes. 

2k 


(17) 


If  thermal  oxide  or  metallization  runs  are  to  he  included  in  the  problem  Ax 
is  of  the  order  of  tenths  of  microns.  Then  At  just  be  smaller  than  tenths  of 
nanoseconds.  If  the  final  solution  is  needed  for  times  of  the  order  of 
milliseconds  not  only  are  enormous  computing  times  required,  but  the  solution 
becomes  suspect  due  to  numerical  inaccuracies. 

Thus  in  certain  situations  the  stability  condition  proves  to  be 
such  a  limitation  that  other,  less  intuitive,  but  more  powerful  methods  are 
desirable.  Before  discussing  one  such  method  we  will  digress  for  a  moment  to 
discuss  general  techniques  for  numerical  solutions  to  the  heat  equation. 

The  finite  difference  method  is  the  oldest  and  most  widely  used 
approach  to  the  solution  of  heat  conduction  problems.  Monte  Carlo  methods 
are  occassionally  used,  and  the  finite  element  method  is  a  powerful  way  of 
handling  a  certain  class  of  problems.  Neverthelss,  finite  differences  have 
many  desirable  features  and  are  usually  used.  Many  excellent  texts  are 
available  on  the  subject. 13* lu 

The  heat  diffusion  equation  is  one  type  of  parabolic  differential 
equation.  The  finite  difference  approach  aims  to  approximate  the  first 
order  derivative  and  the  second  order  derivative  by  means  of  Taylor  series 
expansions  and  obtain  algebraic  equations.  Rewriting  the  heat  equation  in 
its  linearized  one -dimensional  form 


3t  Pc 


(18) 


T(x,t)  can  be  expanded  in  the  forward  direction 
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(19) 


T(x+Ax,t)  =  T(x,t)  +  Ax  ||  ♦  0  + 


or  in  tne  backyard  direction 


T(x-Ax,t)  -  T(x,t)  -  Ax||+  .... 


(20) 


If  we  neglect  higher  order  terms.  Equations  (19)  and  (20)  can  be  solved  to 
yield  the  central  difference  approximation  to  the  second  derivative 


—  •  7-~r  [T(x+Ax,t)  -  2T(x,t)  +  T(x-Ax,t)]  . 
3x^  (Ax)^ 

The  forward  difference  approximation  of  the  time  derivative  is 

||=  ~  [T(x,t+At)  -  T(x,t)]  . 


(21) 


(22) 


Equations  (21)  and  (22)  lead  directly  to  the  explicit  heat  flow  equation, 
namely  Equation  (l6). 

The  backward  difference  approximation  of  the  time  derivative  is 
f£  =  "  T(x,t-At)}  •  (23) 


This  leads  to  a  somewhat  different  finite  difference  formulation  of  Equation 
(18),  namely 

Tn(t+4t)-Tn(t>  *  'u  (Vl(tt4t>-2Tn(t+4t)tTn-l(tt4t,lt1Ato)-(2l,> 

The  difference  is  that  the  rise  in  temperature  at  a  node  is  written  in  terms 
of  the  new  node  temperatures,  as  yet  unknown.  We  are  lead  to  write  a  group 
of  n  simultaneous  equations.  This  is  called  the  implicit  method  or  Liebman's 
method.  There  seems  to  be  no  particular  advantage  to  this  procedure  until 
we  realize  that  there  is  no  stability  problem.  The  method  is  inherently 
stable  for  any  time  step  size.  There  are  still  questions  of  accuracy  to  be 
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investigated,  but  not  of  stability.  Using  this  approach  one  can  avoid 
stability  problems  and  choose  a  value  of  At  based  on  other  considerations 
such  as  accuracy. 

A  number  of  variations  of  the  previous  development  exist  for  the 
purpose  of  speeding  convergence  or  increasing  accuracy.  These  typically  con¬ 
sist  of  a  weighted  average  of  forvard  and  backward  differences.  See  refer¬ 
ence  111  for  further  details. 

Convergence  and  accuracy  are  problems  of  which  little  of  a  quanti¬ 
tative  nature  can  be  said.  In  general,  if  At  and  Ax  are  made  sufficiently 
small  the  numerical  solution  will  approach  the  true  solution.  Trial  and 
error  must  often  be  resorted  to.  Having  solved  a  problem  the  solution  can 
be  compared  to  that  of  the  same  problem  by  halving  At  and  Ax.  If  the  change 
in  the  solution  is  small  one  is  usually  Justified  in  assuming  that  the  error 
is  small. 

The  main  drawback  of  the  implicit  method  is  the  need  to  invert  the 
matrix  -  repeatedly  if  variable  coefficients  are  to  be  allowed.  Fortunately 
the  matrix  is  sparse.  Gauss-Siedel  or  some  similar  iterative  method  of 
solution  is  usually  recommended  to  solve  the  simultaneous  equations.  Unfor¬ 
tunately,  round-off  error  can  sometimes  be  a  problem.  Since  the  matrix  has 
many  null  terms  it  is  worthwhile  to  consider  other,  more  direct  methods  than 
iterative  schemes. 

Richtmeyer lu  describes  a  very  fast  way  of  solving  the  set  of 
Equations  ('A).  Each  equation  has  at  the  most  three  unknown  temperatures, 
and  can  be  written  in  the  form 


,  T,  +  A  'I,  +  =  C. 

Jj-l  J-l  J .1  2  J+l  2 
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If  we  write 


;  *  i  T  +  T 
2  2  2*1  *  Fi 


Then  by  substitution  into  (25) 


.  «  k22*l 


T  -  cAhVi 
4  am+am-iej-i 


(26) 


(27a) 
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When  J 


lt  Equation  (25)  has  only  2  unknown  temperatures  and  and 


!L.  tii  T 

*n  2 


so  that 
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and 
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(28) 


How  all  values  of  E^  and  F^  can  be  calculated  from  (27). 

If  we  assume  that  the  nth  node  remains  at  T  «  T  ,  then  the  last 

n  o 

equation  of  (25)  can  be  solved  for  T  , . 


n-l 


E  ,T 
n-l  o 


+  F 


n-l 


How  Equation  (26)  will  give  all  the  unknown  node  temperatures . 

This  algorithm  was  found  to  be  much  more  accurate  than  the  Gauss- 
3eidel  method  for  coefficients  encountered  in  a  realistic  problem.  In 
addition  it  vas  faster  by  a  factor  of  twenty  times  in  certain  problems. 

A  computer  program  was  written  to  solve  one-dimensional  problems 
in  heat  conduction.  A  program  listing  appears  in  Appendix  I.  The  program 
is  called  THAP  -  for  Transient  Heat  Analysis  Program.  It  makes  use  of  the 
implicit  method  and  Kichtmeyer's  algorithm.  As  presently  written  the 
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program  will  solve  problems  with  rp  to  20  nodes  of  arbitrary  sizes  consisting 
of  silicon,  silicon  dioxide,  and  aluminum  in  any  order.  The  initial  tempera¬ 
ture  must  be  uniform,  but  can  have  any  value.  The  thermal,  conductivity  and 
specific  heat  of  all  three  materials  are  functions  of  temperature. 

Examples  of  problems  solved  with  THAP  are  shown  in  Figures  13  and 
lb.  Figure  13  shows  * be  temperature  rise  in  a  silicon  chip  when  power  is 
supplied  uniformly  throughout  a  surface  layer  b  microns  thick.  The  power 
density,  1.5  x  10s  watts/cm2-ym,  corresponds  to  50  watts  over  an  area  about 
90  microns  by  90  microns.  The  peak  temperature  increases  indefinitely  in 
this  example,  in  fact  reaching  the  melting  point  of  silicon  in  about  5  usee. 
In  1  ysec  the  temperature  at  x  ■  10  ym  is  about  25#  of  T  .  This  should  not 
be  taken  to  mean,  however,  that  some  other  boundary  condition  at  x  *  10  ym, 
such  as  a  plated  heat  sink,  would  have  a  pronounced  effect  on  T  .  It  would 
have  some  effect,  of  course,  but  quantitative  answers  can  only  be  obtained 
by  running  the  problem  again. 

One  other  precaution  should  be  pointed  out  when  using  the  results 
of  computer  calculations.  In  Figure  13  the  axes  are  not  normalized.  Since 
the  problem  is  non-linear  the  usual  normalized  parameters  cannot  be  used,  and 
caution  must  be  exercised  when  attempting  to  extrapolate  the  curves  to  fit 
new  situations.  If,  for  example,  the  temperature  rise  is  b0O°C  above  ambient 
in  time  tQ,  the  r*se  is  not  doubled  if  the  input  power  is  doubled.  In  fact 
the  temperature  rise  will  be  more  than  doubled  in  a  typical  case. 

The  same  problem  has  been  rerun  in  Figure  lb,  except  that  the  sur¬ 
face  has  been  covered  by  0.5  microns  of  silicon  dioxide  with  1  micron  of 
aluminum  on  top  of  that.  It  is  evident  that  the  added  thermal  capacitance 
of  the  surface  layers  ha3  lowered  T  appreciably,  except  for  times  less  than 
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Fig.  13.  Temperature  distribution  in  silicon  due  to  power  applied  in  a 
4  ym  layer  at  the  surface.  Solution  is  by  the  implicit  finite 
difference  method  using  Riehtayer’s  algorithm  and  temperature 
dependent  thermal  properties. 
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Fig.  l^t.  Temperature  distribution  under  the  same  conditions  as  Fig,  13* 
but  with  a  0.5  um  layer  of  silicon  dioxide  and  a  1  ym  layer 
of  aluminum  at  the  surface. 


0.5  ysec.  In  the  passive  silicon,  x>5.5  ym,  the  temperature  is  nearly 
unaffected  by  the  added  layers,  at  least  for  the  times  shown. 

A  third  example,  not  plotted  here,  was  exactly  the  same  as  the 
problem  in  Figure  14,  except  that  the  power  was  dissipated  in  a  1  micron 
layer  between  3  and  4  microns  from  the  surface,  simulating  power  dissipation 
in  a  depletion  layer  beneath  the  surface.  The  temperature  distribution  was 
somewhat  different,  of  course,  but  T  was  almost  Identical.  In  other  words, 
fine  detail  in  the  power  distribution  was  not  important.  A  more  quantitative 
estimate  of  the  importance  of  detail  in  the  model  is  obtained  by  calculating 
the  diffusion  length  *  (Dt)  ,  where  D  is  the  diffusion  constant  and  t  is 
the  time.  For  silicon  D  is  about  1  cm2/sec.  If  t  ■  1  ysec.  *>  10  ym.  For 
a  time  period  of  1  usee  heat  spreads  through  distances  of  about  10  urn.  This 
was  observed  in  the  case  of  Figure  13.  It  explains  why  the  exact  distrlbutlcn 
of  power  within  a  layer  less  than  10  urn  thick  is  not  important  for  t  of  the 
order  of  1  usee. 

Further  details  of  the  effect  of  single  surface  layers  on  peak 
temperature  are  given  in  Figure  15.  Here  again  we  have  a  one-dimensional 
problem,  with  a  power  density  of  1.5  x  10s  watts/cm2-ym  in  a  4  ym  layer.  The 
peak  temperature  is  plotted  as  a  function  of  time.  The  dotted  lines  corre¬ 
spond  to  different  thicknesses  of  SiOg  on  the  surface.  The  line  labelled 
5  urn  8i0g  also  holds  for  any  thickness  greater  than  5  ym  in  the  time  scale 
given.  This  implies  that  glassivation  layers  5  ym  or  thicker  insulate  the 
silicon  from  any  possible  heat-sinking  effects  of  metal  deposited  on  top  of 
glass  (for  pulses  less  than  about  10  ysec). 

The  solid  lines,  terminated  at  660°C  shov  the  effects  of  aluminum 
deposited  directly  on  the  silicon.  The  decrease  in  is  quite  pronounced. 
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Fig.  15.  Effect  of  surface  layers  on  peak  temperature  rise  in  silicon. 
Power  is  applied  in  a  h  urn  thickness,  and  the  layers  are 
either  silicon  dioxide  or  aluminum. 
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For  tines  less  than  6  Msec  the  50  Mm  layer  is  the  same  as  an  infinite  layer 
and  gives  a  rough  idea  of  the  effect  of  a  plated  heat  sink  or  a  Schottky 
barrier  on  T  .  Since  copper  and  gold  have  greater  heat  conductivity  than 

Xu 

aluminum  T  would  be  even  lower  than  given  in  the  figure, 
m 

More  detailed  plots  near  the  origin  show  that  the  SiOg  layer  has  a 

negligible  effect -on  T  when  t  is  less  than  100  nsec,  and  the  same  is  true 

m 

for  aluminum  when  t  is  less  than  20  nsec . 

When  both  aluminum  and  SiOg  are  present,  the  results  are  more  com¬ 
plicated.  Figure  l6  shows  the  peak  temperature  after  1.5  Msec  as  a  function 
of  the  oxide  thickness  with  the  aluminum  thickness  as  a  parameter.  As  ex¬ 
pected,  the  effect  of  the  aluminum  becomes  less  and  less  as  the  oxide  thick¬ 
ness  is  increased.  A  1  urn  thickness  of  Si0_  is  nearly  the  same  as  an  in- 

d 

finite  layer,  as  can  be  verified  in  Figure  15.  Of  special  interest  is  the 

fact  that  a  1  pm  layer  of  aluminum  lowers  Tm  to  about  500°C  irrespective  of 

the  oxide  thickness.  This  rather  surprising  result  means  that  a  typical 

metallization  run  has  about  the  same  heat  sinking  capability  whether  it  is 

over  emitter  oxide  or  over  very  thick  oxide!  Again  it  should  be  pointed  out 

that  this  observation  is  only  valid  for  one-dimensional  heat  flow  and  the 

other  particular  circumstances  of  Figure  16. 

In  Figures  17-19  the  power  density  to  reach  a  certain  T  as  a 

m 

function  of  time  is  plotted.  The  curves  are  probably  more  useful  in  in¬ 
verted  form,  that  is,  given  a  power  density,  how  long  does  it  take  to  reach 
a  certain  temperature.  The  three  figures  are  for  power  inputs  in  a  2,  ?,  and 
h  um  layer,  respectively.  This  is  the  range  of  base  diffusions  in  most  of 
the  circuits  tested  experimentally.  A  layer  of  silicon  dioxide  0.!*  um  thick 
Covers  the  surface,  this  being  the  average  between  base  and  emitter  oxide 
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Fig.  16.  Peak  temperature  rise  in  silicon  after  1.5  ysec  with 
variable  thicknesses  of  silicon  dioxide  and  aluminum. 
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Fig.  19.  Seme  as  Figvre  17*  with  the  power  applied  to  a  U  ym  layer. 
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thicknesses  of  typical  microcircuits. 

Comparison  of  the  three  figures  is  obscured  by  the  fact  that  the 
save  pover  density  is  a  different  total  pover  input  in  the  three  cases.  The 
thicker  the  "active”  leyer  the  greater  the  total  pover  input  and  the  greater 
the  temperature  rise.  For  times  less  than  50  nsec  the  heating  is  nearly 
adiabatic  and  T  does  not  depend  on  layer  thickness.  The  three  figures  coin- 
cide  in  this  range  and  the  pover  density  vs  time  curve  varies  as  t”1  as 
ohovn  in  Figure  19.  For  longer  times  the  curves  vould  be  expected  to  vary  as 
t”*5  vhen  heat  transfer  to  the  nearby  silicon  becomes  appreciable.  In  Fig- 
ure  17  the  t  asoymptote  has  been  reached,  vhile  for  thicker  layers  the 
quantity  of  heat  in  the  active  layer  if  still  large  compared  vith  the  heat 
energy  conducted  avay  and  the  situation  is  still  adiabatic  -  like  in  the  time 
scale  of  the  figures. 

These  last  three  curves  are  useful  for  estimating  the  temperature 
vhich  a  Junction  reaches  vhen  it  goes  into  second  breakdown.  Unfortunately, 
the  one-dimensional  solution  is  a  limitation.  In  an  integrated  circuit,  the 
lateral  heat  flov  into  surrounding  silicon  is  important  vhen  t  is  greater 
than  about  100  nsec. 

To  summarize  the  presentation  in  this  section,  we  can  say  that 
conduction  is  the  only  important  mechanism  of  heat  flow  in  integrated  cir¬ 
cuits.  Although  analytic  solutions  to  the  heat  flows  equation  are  available 
for  a  wide  variety  of  relevant  conditions,  they  cannot  be  applied  with  confi¬ 
dence  when  the  temperature  rise  is  large  because  of  the  variation  in  the 
thermal  conductivity  of  silicon.  Accurate  answers  are  obtained  by  numerical 
methods  using  finite  difference  techniques.  Solutions  are  best  expressed  in 
the  form  of  power  versus  time  relationships  which  can  then  be  compared  to 
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experimental  findings. 

In  the  next  section  a  brief  survey  of  second  breakdown  in  semicon¬ 
ductor  devices  vill  be  given.  This  will  serve  as  an  introduction  to  the  ex¬ 
perimental  results  presented  and  discussed  in  later  sections. 
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SECTION  III 


A  SURVEY  OF  SECOND  BREAKDOWN  IN  INTEGRATED  CIRCUITS 

The  effect  in  bipolar  transistors  commonly  referred  to  as  "second 
breakdown"  was  first  reported  in  1958  by  Thornton  and  Simmons.15  They  sug¬ 
gested  that  it  could  be  the  cause  of  certain  failures  occuring  in  electronic 
circuits.  It  appears  as  a  sudden  transition  to  a  low  voltage,  high  current 
mode  of  operation  accompanied  by  current  constriction,  which  may  or  may  not 
cause  permanent  changes  in  transistor  characteristics  or  even  catastrophic 
failure. 

Although  the  term  secord breakdown  was  initially  coined  to  re¬ 
fer  to  bipolar  transistors  operating  in  the  normal,  common-emitter 
configuration  with  the  base  terminal  open,  forward  biased,  or  reverse 
biased,  similar  phenomena  occur  in  diodes,  both  forward  and  reverse  biased, 
MOS  devices,  microwave  diodes,  4-layer  devices,  and,  indeed,  in  all  types  of 
semiconductor  devices.  The  very  large  number  of  papers  published  on  second 
breakdown  and  related  areas,  the  varied  description  of  events  which  accom¬ 
pany  second  breakdown,  and  the  proliferation  of  models  and  explanations  for 
the  phenomenon  attest  to  its  elusiveness  and  complexity.  This  section  will 
describe  briefly  the  historical  development  and  features  of  the  popular 
models,  and  some  of  the  more  recent  experimental  work. 

As  mentioned  above,  the  first  observations  of  second  breakdown 
were  in  1958.  The  previous  year,  Tauc  and  Abraham  reported  the  thermal 
breakdown  behavior  of  silicon  p-n  junctions.16  These  and  later  papers 
offered  a  host  of  explanations.  The  similarity  between  breakdown  in 
Junctions,  transistors,  and  other  devices  led  to  proposals  of  a  common  under¬ 
lying  mechanism.  These  included  the  formation  of  a  central  microscopic  melt, 
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a  critical  temperature,  a  critical  current,  a  lateral  thermal  instability, 
a  pinch  effect,  an  instability  due  to  space  charge  effects,  and  others. 
Special  issues  of  the  IEEE  Transactions  on  Electron  Devices,  August/September 
and  November,  1966,  were  devoted  to  second  breakdown.  In  1967  a  comprehen¬ 
sive  survey  summarized  and  critically  reviewed  the  theories  in  vogue  at  that 
time. 17 

There  are  two  common  methods  of  observing  the  transition  into 
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second  breakdown.  One  is  by  direct  observation  of  the  volt-ampere  charac¬ 
teristic  on  a  curve  tracer.  Figure  20  shows  typical  curves  for  a  transistor 
and  a  diode.  The  transistor  is  operated  in  the  grounded  emitter  connection 
with  the  base  either  forward  (F),  open  (0),  or  reverse  (R)  biased.  Schafft17 
points  out  that  a  criterion  sometimes  used  for  second  breakdown  is  that  the 
voltage  after  breakdown  is  lowest  for  R  and  highest  for  F.  The  diode  in 
part  (b)  of  the  figure  is  assumed  to  be  operating  in  the  reverse  biased 
mode,  although  similar  results  are  obtained  for  a  forward  biased  diode. 

In  the  second  method,  a  voltage  or  current  pulse  is  applied  to  the 
device  under  test.  If  the  amplitude  is  great  enough,  after  a  certain  delay 
time,  which  can  be  as  short  as  nanoseconds,  a  transition  to  a  new  voltage 
and  current  level  will  occur.  Prior  to  the  transition  there  is  often  evi¬ 
dence  of  heating  such  as  a  rise  in  voltage  or  decline  in  current.  After  the 
transition  the  voltage  and  current  levels  are  usually  quite  stable.  Typical 
waveforms  are  depicted  in  Figure  21. 

Following  an  excursion  into  second  breakdown  permanent  changes  in 
the  device  may  be  observed,  depending  on  the  structure,  heat  sinking,  length 

*  During  the  remainder  of  this  report  the  term  second  breakdown  will  be  used 
to  inPo.i  the  sudden  change  to  a  low  voltage,  high  current  mode  of  operation 
of  any  semiconductor  device. 
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of  tine  In  second  breakdown*  and  the  power  dissipation  before  and  after  the 
transition.  A  microscopic  examination  may  reveal  no  observable  damage* 
darkening  of  the  contacts*  a  stress  mark  at  a  junction*  "white  streak" 
damage*  or  aassive  Junction  and  metallization  damage.  Electrical  changes 
may  not  be  detectable,  or  may  include  increased  leakage  current*  soft  junc¬ 
tion  brerJcdown*  changes  in  breakdown  voltage*  changes  in  beta*  a  decrease  in 
delay  time  during  subsequent  pulsing*  or  the  presence  of  a  low  resistance 
short  circuit  across  the  contacts. 

For  devices  where  the  delay  time  depends  on  the  amplitude  of  the 
poise*  it  is  clear  that  a  thermal  mechanism  must  be  Involved.  The  concept  of 
an  intrinsic  temperature  was  proposed  by  Ford18  and  by  Melchior  and  Strutt.19 
In  this  theory,  once  the  junction  temperature  reaches  a  point  where  the 
lightly  doped  side  of  the  junction  becomes  intrinsic*  then  further  increases 
in  temperature  will  cause  a  decrease  in  resistivity.  The  negative  resistance 
results  in  an  instability  which  causes  a  current  constriction  to  form.  The 
low  resistance  of  the  constriction  causes  the  collapse  of  the  Junction  field 
and  the  drop  in  terminal  voltage. 

Another  version  invokes  a  critical  temperature  at  the  point  where 
thermal  generation  of  carriers  is  able  to  satisfy  the  current  requirements 
without  the  need  for  avalanche  multiplication.  When  the  device  becotr  s  hot 
enough  the  current  concentrates  at  the  point  of  highest  temperature  and  the 
extinguishing  of  the  avalanche  produces  the  observed  drop  in  voltage.  If 
this  is  the  case*  the  presence  of  defects  should  have  an  effect  on  the  delay 
time.  Investigations  into  this  point  have  produced  mixed  results  according 
to  Schafft . 1 7  His  conclusion  is  that  second  breakdown  does  not  depend  on 
imperfections  for  its  existence  but  is  enhanced  by  tb'im. 
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If  the  idea  of  a  critical  temper at  lire  is  valid,  changes  in  ambient 
temperature  should  cause  a  difference  in  the  delay  time.  Again  the  results 
are  mixed  as  repo-ted  by  Schafft. 

In  the  microscopic  central  melt  theory20  it  was  postulated  that  a 
small  portion  of  the  junction  melts  during  the  transition  to  second  break¬ 
down.  If  the  melt  or  "mesoplaama"  were  small  enough  it  might  not  be  detected 
during  an  analysis  of  the  device.  A  number  of  other  workers  have  investi¬ 
gated  the  temperature  and  voltage  distribution  in  melted  columns,  the 
stability  of  melted  regions ,  and  the  voltage-current  relationships . 

More  recently  a  number  of  other  important  contributions  to  tha  work 
on  second  breakdown  have  been  made.  Unfortunately,  the  picture  has  not  been 
simplified;  instead,  new  evidence  has  been  presented  for  the  existence  of 
several  distinct  mechanisms. 

In  the  case  of  diodes,  Bowers21  and  Muller  and  Guckel22  have  shown 
that  space  charge  effects  due  to  device  current  can  cause  negative  resistance 
under  certain  conditions.  Apart  from  the  instability  which  accompanies  nega¬ 
tive  resistance,  the  space  charge  will  cause  a  change  in  voltage  across  the 
space  charge  region  which  then  appears  as  a  change  in  terminal  voltage.  With 
high  current  pulses  the  voltage  is  considerably  higher  than  the  junction 
breakdown  voltage  and  it  is  important  to  know  how  much  of  the  excess  appears 
across  the  Junction  and  how  much  is  due  to  parasitic  IR  drops  in  the  bulk 
material,  leads,  etc. 

Hower  and  Reddi23  have  reported  that  second  breakdown  can  be  in- 
duced  in  n  -n-n  diodes  by  a  process  of  avalanche  injection.  The  diodes  were 
sufficiently  small  and  the  pulse  power  and  pulse  time  such  that  the  peak 
temperatures  did  not  exceed  38°C  above  ambient.  Under  these  conditions,  if 
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second  breakdown  is  observed,  it  is  easy  to  argue  that  the  effects  cannot  be 
thermal!  They  observed  a  negative  resistance  region  at  low  currents  which 
was  attributed  to  space  charge  effects.  At  higher  currents  a  discontinuity 
to  a  low  voltage  high  current  state  was  observed.  This  was  correlated  to 
the  onset  of  avalanche  multiplication  at  an  n+-n  junction.  Trariistors 
fabricated  on  the  same  wafer  with  diodes  were  also  found  to  have  second 
breakdown  characteristics  caused  by  avalanche  injection. 

Second  breakdown  triggered  by  a  critical  current  level  was  found 
by  Tasca  in  both  diodes  and  transistor  junctions.24  He  did  not  attempt  to 
fit  hi 8  findings  to  any  particular  mechanism.  On  the  other  hand,  he  also 
reported  that  some  of  his  devices  failed  by  a  thermal  mechanism. 

Chiang  and  Lauritzen25  worked  with  small,  heavily-doped,  guard 
ring  diodes.  They  monitored  the  light  emission  due  to  avalanche  breakdown 
and  incandescence  at  current  constrictions.  Thermal  breakdown  occurred  at 
about  1100°C,  which  was  also  the  temperature  at  which  the  lightly  doped  side 
became  intrinsic.  The  injected  carrier  concentration  at  thermal  breakdown 
was  about  two  orders  of  magnitude  less  than  the  background  doping  concen¬ 
tration. 

Fleming6  and  Wang26  have  argued  the  case  for  the  existence  of  a 
critical  temperature  where  the  leakage  current  has  increased  enough  to 
provide  the  total  device  current.  Avalanche  multiplication  is  no  longer 
required  end  the  avalanche  is  extinguished.  This  point  of  view  is  also 
advanced  by  Sunshine  who  has  performed  a  series  of  detailed  and  comprehen¬ 
sive  experiments  on  breakdown  in  silicon-on-sapphire  diodes.27 

Sunshine's  work  is  an  extension  of  work  by  Dunlin28  who  had  con¬ 
cluded  that  second  breakdown  was  reached  when  the  diode  current  density  was 
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equal  to  the  background  doping  concentration.  Sunshine,  using  SOS  diodes 
supplied  by  Dumin,  devised  careful  experimental  techniques  vhich  enabled  him 
to  ascertain  the  sequence  of  events  which  accompanied  second  breakdown.  In 
the  first  place,  he  reexamined  Dumin 's  work  on  light  emission  during  ava¬ 
lanche  breakdown,  current  constriction  at  the  junction,  and  filament  forma¬ 
tion  in  the  high  resistivity  region.  He  also  assembled  light  transmission 
apparatus  for  determining  the  temperature  of  the  silicon  material  based  on 
the  variation  of  the  optical  absorption  coefficient  with  temperature. 
Finally,  he  used  a  laser  in  a  stroboscopic  technique  which  enabled  him  to 
observe  temperature  variations  in  space  and  time  with  repetitive  pulses. 

The  diodec  were  made  by  growing  a  thin  silicon  film,  either  n  or 
p  type,  about  1  urn  thick  on  r,  sapphire  substrate.  p+  and  n+  regions  were 
diffused  to  form  p  nn  or  p  pn  diodes.  Metel  lands  on  the  p  and  n 
regions  served  as  contacts  for  probes.  The  diodes  were  about  50  ym  wide, 
and  the  high-resistivity  region  was  varied  in  length  from  2  to  300  ym. 
Breakdown  voltage  depended  on  doping  level,  being  typically  18  volts. 

The  dc  experiments  clearly  showed  the  relationship  between  the 
diode  volt-ampere  characteristic  and  the  current  distribution  in  the  diode. 
At  low  current  levels  avalanche  light  emission  showed  nonuniform  breakdown 
across  the  Junction.  The  silicon  had  a  high  dislocation  density,  but  was 
otherwise  free  from  defects.  In  spite  of  the  non-uniform  breakdown  the 
volt-ampere  curves  showed  a  sharp,  clean  breakdown  unlike  that  when  micro- 
plasmas  are  present .  As  the  current  level  incre  ised  the  avalanche  emission 
became  more  uniform  and  more  intense,  shoving  the  levelling  effect  of  the 
negative  temperature  coefficient  of  avalanche  breakdown.  At  still  higher 
currents  the  current  began  to  constrict,  usually  near  the  center  of  the 


diode  where  the  temperature  was  greatest.  At  a  level  corresponding  to  the 
threshold  current  for  the  onset  of  negative  resistance  (second  breakdown)  the 
avalanche  emission  was  extinguished  in  a  portion  of  the  junction.  At  the 
center  of  this  region  a  small  incandescent  spot  appeared.  As  the  current 
was  increased  still  further,  the  spot  grew  in  size  and  intensity  and  ava¬ 
lanche  emission  ceased  everywhere. 

A  constant  current  source  was  also  used  for  pulse  testing.  When 
the  pulse  current  was  less  than  the  threshold  value  the  diode  voltage  would 
show  the  effects  of  diode  heating.  When  the  current  exceeded  the  threshold 
value,  two  characteristic  time  intervals  were  noted:  T^,  the  delay  time,  and 
Tp,  the  filament  formation  time.  When  the  current  pulse  -ms  applied  the 
voltage  would  initially  increase  due  to  heating.  At  time  the  voltage 
reached  a  maximum  value,  then  began  to  decline  rapidly,  After  an  additional 
time  interval  Tp  the  voltage  reached  a  steady,  sustaining  value. 

With  the  stroboscopic  technique  Sunshine  was  able  to  correlate 
these  characteristic  times  with  variations  in  current  and  temperature  in  the 
diode.  At  time  two  things  happened.  First,  avalanche  emission  ceased  in 
a  region  near  the  center  of  the  diode.  Second,  a  hot  spot  formed  at  the 
center  of  the  region,  while  portions  of  the  junction  adjacent  to  the  hot 
spot  began  to  cool.  After  careful  investigation,  Sunshine  concluded  that 
the  extinction  of  the  avalanche  and  the  formation  of  tie  hot  spot  occurred 
when  the  local  temperature  was  high  enough  to  generate  sufficient  thermal 
current  to  supply  the  diode  current  density.  Avalanche  multiplication  being 
no  longer  required,  the  voltage  across  the  Junction  at  that  point  would  drop 
Tne  drop  in  voltage  would  cause  the  nearby  current  to  funnel  into  the  region, 
resulting  in  a  current  constriction. 
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The  constriction  would  elongate  into  a  filament  extending  across 
the  high-resistivity  region  during  the  time  interval  Tj,  if  the  growth  condi¬ 
tions  were  favorable.  This  required  that  the  power  density  he  high  enough 
to  raise  the  temperature  of  the  high-resistivity  material  to  its  intrinsic 
temperature . 

Constriction  occurred  at  Junction  temperatures  of  several  hundred 
degrees  Celcius.  As  the  filament  formed  and  became  incandescent  the  peak 
temperature  increased,  reaching  700-1100°C  before  burnout  occurred. 

In  the  case  of  forward  biased  diodes  the  maximum  temperature  rise 
occurred  near  the  center  of  the  high  resistivity  region  rather  than  at  the 
Junction.  Filament  formation  again  took  place  when  the  temperature  reached 
the  intrinsic  value. 

Sunshine's  work  has  been  extended  and  elaborated  on  by  Budenstein, 
Pontius,  and  Smith.29  Their  conclusions  generally  agree  with  those  of  Sun¬ 
shine.  They  used  SOS  diodes  of  various  geometries  and  resistivities  fabri¬ 
cated  by  the  Autonetics  Division  of  North  American  Rockwell. 

The  diodes  were  pulsed  with  a  constant  current  source.  At  low 
pulse  amplitudes  the  avalanche  emission  was  low  in  intensity  and  beaded  in 
appearance.  There  was  little  indication  from  the  voltage  waveform  or  from 
the  transmitted  light  photograph  of  heating.  At  higher  pulse  amplitudes  the 
avalanche  was  brighter  and  more  uniform  with  evidence  of  heating  at  the 
Junction  and  to  a  lesser  extent  throughout  the  high  resistivity  region.  At 
still  higher  levels,  at  Junction  temperatures  between  600  and  700°C,  current 
constriction  occurred  at  certain  sites  along  the  Junction.  The  constriction 
sites  appeared  in  the  transmission  photographs  as  dark  spots  along  the 
Junction,  with  avalanche  emission  ceasing  in  a  region  1  to  5  mils  on  either 
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side.  The  dark  spots  usually  occurred  at  the  same  places  at  each  pulse.  The 
number  and  spacing  vas  determined  hy  the  pulse  amplitude  and  duration.  The 
short,  high  amplitude  pulses  produced  the  largest  spot  density. 

The  appearance  of  a  current  constriction  vas  accompanied  hy  a 
slight,  often  imperceptible  drop  in  voltage.  As  the  filament  spread  across 
the  high  resistivity  region,  the  voltage  decreased  slovly.  When  the  filament 
finally  bridged  the  entire  region  catastrophic  failure  occurred  as  a  molten 
column  about  1  ym  vide  formed  vithin  the  current  filament.  At  this  point 
the  voltage  dropped  abruptly  to  about  one-half  its  peak  value.  The  damage 
produced  during  this  final  stage  vas  visible  and  also  appeared  as  a  deteri¬ 
oration  of  the  volt-ampere  characteristic.  Budenstein  suggests  that  vhat 
Sunshine  refers  to  as  second  breakdown  is  the  initial  drop  in  voltage  during 
the  nucleation  and  filament  formation  and  that  the  final,  abrupt  change  in 
voltage  accompaniec  by  device  damage  did  not  occur.  Thus  there  is  a  some- 
vhat  different  interpretation  of  second  breakdown  betveen  the  two  groups  of 
researchers . 

Many  other  features  of  Sunshine's  vork  were  corroborated  by  Buden¬ 
stein  et  al.  Nucleation  of  a  filament  was  found  to  occur  when  the  thermally 
generated  carriers  were  able  to  sustain  the  current.  The  filament  spread 
into  the  high-resistivity  region  as  the  temperature  was  raised  to  the  intrin¬ 
sic  temperature. 

In  forward  biased  junctions  or  semiconducting  material  without 
junctions  a  similar  sequence  of  events  took  place.  Heating  occurred  in  the 
high  resistivity  bulk  material  or  around  the  probe  contacts.  A  current 
filament  formed,  and  when  it  finally  bridged  the  contacts  a  molten  filament 
formed  and  voltage  dropped  abruptly.  In  some  reverse  biased  diodes  it 
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was  possible  to  initiate  a  filament  in  the  high  resistivity  region  before 
constriction  at  a  junction  site  took  place. 

The  summary  by  Budenstein  et  al. ,  provides  us  with  a  fairly  clear 
picture  of  second  breakdown  in  junctions.  Two  stages  of  filament  growth  are 
associated  with  second  breakdown.  In  the  first  stage  current  constriction 
occurs  and  a  filament  grows,  accompanied  by  a  relatively  slow,  gradual  drop 
in  voltage.  In  the  final  stage  the  filament  has  bridged  the  contacts  and  a 
melt  column  forms,  accompanied  by  a  precipitous  drop  in  vo^t^ge.  These  two 
stages  cannot  always  be  clearly  distinguished  in  the  voltage  waveform. 

In  the  case  when  a  re*;er»t -biased  junction  is  present  nucleation 
of  a  filament  occurs  at  the  Junction  when  the  thermal  generation  of  carriers 
is  sufficient  to  supply  the  current  requirements.  Otherwise,  a  filament 
forms  when  the  temperature  of  the  bulk  material  reaches  the  peak  of  the 
temperature-resistivity  curve. 

Junction  defects  play  only  a  minor  role  in  the  nucleation  of 
filaments.  The  position  of  junction  dark  spots  depends  mainly  on  the 
temperature  distribution:  When  a  melt  channel  has  been  formed,  a  subsequent 
breakdown  does  not  usually  superimpose  a  second  melt  channel  on  the  first. 

Although  only  thermal  mechanisms  were  considered  it  seems  reasona¬ 
ble  that  other  carrier  injection  mechanisms  would  bring  about  second  break¬ 
down  whenever  the  injected  current  density  was  equal  to  the  local  junction 
current  density. 

This  summarizes  the  present  understanding  of  second  breakdown  in 
semiconductor  junctions.  A  great  deal  of  progress  has  been  made  in  recent 
months.  The  following  section  will  present  the  results  of  pulse  testing  of 
integrated  circuits,  showing  the  relation  with  the  theories  of  this  section. 


SECTION  IV 


EXPERIMENTAL  RESULTS 

In  the  experimental  part,  of  this  program,  several  different  inte¬ 
grated  circuits  vere  stressed  into  second  breakdown  with  single  pulses.  The 
object  was  to  determine  the  relationship  between  threshold  power  and  pulse 
duration.  The  pulse  voltage  and  current  waveforms  and  the  device  damage  will 
be  related  to  the  second  breakdown  models  discussed  in  the  preceding  section. 
Recall  that  second  breakdown  is  defined  here  as  the  transition  to  a  low- 
voltage  high-current  state.  In  the  case  of  reverse  biased  Junctions  this 
transition  almost  always  resulted  in  permanent  electrical  changes  and  visible 
damage.  The  rare  exceptions  will  be  pointed  out.  In  the  case  of  transistors, 
the  results  are  somewhat  more  complex. 

A.  Testing  Procedure 

The  integrated  circuits  were  pulsed  into  second  breakdown  by  the 
application  of  approximately  constant  current  pulses  ranging  from  100  nsec 
to  10  msec  long.  The  device  under  test  was  connected  to  a  Cober  Model  6o4A 
pulse  generator.  Amplitude  of  the  pulser  is  continuously  variable  in  three 
ranges  from  0  to  1500  volts  open  circuit.  Currents  up  to  6A  are  available, 
with  a  peak  power  of  9  KW.  Rise  time  is  adjustable  from  *?0  nsec  to  100  nsec. 
Pulse  width  is  determined  by  an  external  pulse  generator  and  can  be  varied 
from  50  nsec  to  dc.  Output  impedance  for  the  0-150  volt  range,  where  most 
testing  was  done,  is  about  4o  ohms.  The  generator  was  found  to  be  reliable, 
simple,  and  convenient  to  use,  and  was  able  to  supply  sufficient  power  for 
nearly  all  test  conditions.  Voltage  across  the  DUT  was  displayed  on  one 
channel  of  a  Tektronix  555  dual  beam  scope  while  the  other  beam  monitored 
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the  current  by  means  of  a  Tektronix  clamp-on  current  probe.  All  testing  vas 
of  a  single  pulse  nature.  A  Datapulse  101  pulse  generator  operating  in  the 
3ingle  pulse  mode  triggered  the  Cober  pulse  generator  and  fixed  the  pulse 
width  while  simultaneously  triggering  the  oscilloscope  sweep.  A  Polaroid 
camera  recorded  the  current  and  voltage  waveforms  of  each  pulse. 

Figure.  22  shows  the  experimental  arrangement.  Most  tests  were 
performed  at  room  temperature .  In  sane  cases  the  ambient  temperature  vas 
changed  by  immersing  the  device  in  liquid  nitrogen  or  mounting  it  on  a  hot 
stage.  In  some  instances  the  DUT  was  mounted  under  a  microscope  with  a  video 
camera  attached.  The  results  of  single  or  repetitive  pulsing  could  be 
viewed  on  a  monitor  and  simultaneously  recorded  on  an  Ampex  VR78OO  video  tape 
recorder . 

Junctions  were  stressed  in  both  forward  and  reverse  biased  modes. 
Transistors  were  tested  collector-emitter,  collector-base,  and  emitter-base. 
Electrical  characteristics  were  checked  on  a  curve  tracer  before  and  after 
testing.  The  devices  were  delidded  for  visual  inspection  after  testing. 

In  a  few  cases  the  devices  were  decapped  and  photographed  prior  to  testing. 

The  usual  procedure  was  to  select  a  pulse  width  and  gradually  in¬ 
crease  the  pulse  amplitude  ;il  second  breakdown  occurred.  Second  break¬ 
down  was  detected  either  rrom  the  voltage-current  waveform  or  from  changes 
in  the  electrical  characteristics  displayed  on  a  curve  tracer.  To  rule  out 
the  possibility  that  the  initial,  sub-threshold  pulses  could  somehow  alter 
the  device's  susceptibility  to  second  breakdown,  several  Junctions  were 
pulsed  50-100  times  or  were  pulsed  repetitively  at  low  prf  for  several 
minutes.  The  voltage-current  waveforms  were  found  to  be  identical  before 
and  after,  the  electrical  characteristics  s.  ->wed  no  evidence  ofchange,  and 
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Fig.  22.  Experimental  arrangement  for  pulsing  devices. 

The  current  and  voltage  waveforms  were  recorded 
by  a  Polaroid  camera. 


the  second  breakdown  characteristics  were  subsequently  found  to  be  virtually 
identical  with  those  of  other  junctions  on  the  sane  chip. 

After  stress-testing,  most  devices  vere  examined  visually  and 
photographed.  Some  circuits  were  examined  in  greater  detail  by  probing, 
angle-lapping,  cross-sectioning,  and  etching. 

B.  Devices  Tested 

All  devices  vere  in  monolithic  integrated  circuit  form,  either 
commercially  available  or  fabricated  by  standard  techniques'  by  major  semi* 
conductor  manufacturers.  Tests  vere  performed  only  on  parts  where  single 
junctions  or  transistors  could  be  Isolated  from  the  rest  of  the  circuit. 

This  placed  some  restriction  on  the  selection  of  devices;  nevertheless,  a 
reasonable  selection  of  manufacturers  and  fabrication  technologies  was  made. 

A  list  of  devices  is  given  in  Table  1.  With  the  exception  of  the 
core  drivers,  all  circuits  vere  low  power  devices.  Up  to  20  parts  of  each 
kind  vere  tested,  each  part  containing  up  to  6  available  junctions  and 
transistors. 

C.  Results 
1.  Part  A 

This  device  vas  tested  more  exhaustively  than  any  of  the  others. 

It  is  a  dual,  3-input  expander,  mounted  in  a  10  lead  T0-5  can.  The  circuit 
diagram  is  given  in  Figure  23.  The  connected  emitters  and  collectors  pre¬ 
clude  the  possibility  of  testing  individual  transistors  unless  the  metalli¬ 
sation  is  scratched  open.  If  the  collector-base  junction  of  one  transistor 
is  pulsed  it  the  reverse  biased  direction,  instead  of  observing  BV^^q  on 
that  transistor  ve  observe  BV_  ,n  on  one  of  the  adjacent  transistors  plus  the 
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Fig.  23.  Pin  connections  and  schematic  diagram  for  Part  A. 
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emitter-base  bree&dcwn  voltage  of  jbhe  transistor  being  tested.  Consequently 
only  the  emitter-base  Junctions  were  tested. 

Photomicrographs  of  the  circuit  appear  in  Figure  2b.  The  chip 
size  is  3b  x  5b  mils  x  7  mils  thick.  Details  of  the  emitter-base  region  are 
shown  in  part(c)of  the  figure,  where  each  small  division  is  2.5  ym.  Note 
the  misregistration.  The  distance  between  the  emitter  Junction  and  the  base 
contact  was  found  to  'vary  from  5  to  17  pm  in  the  17  parts  examined.  A  corre¬ 
lation  was  found  between  this  spacing  and  the  threshold  power  for  second 
breakdown.  Also  note  that  some  of  the  contact  windows  are  only  partly 
covered  with  aluminum,  the  emitter  contacts  on  one  side  and  the  base  con¬ 
tacts  on  the  other.  This  has  implications  regarding  the  effect  of  contact 
resistance  and  contact  alloying  on  second  breakdown  threshold  energy.  No 
relationship  was  noted  on  these  devices. 

Angle  lapping  and  cross  sectioning  several  samples  yielded  the 
dimensions  shown  in  Figure  25.  There  is  no  epitaxial  layer,  the  collector 
regions  being  formed  from  a  deep  diffusion.  The  oxide  thicknesses  were  of 
the  order  of  0.15  ym  over  the  emitter  regions,  o  ym  over  the  base  dif¬ 
fusions,  and  0.63  ym  over  the  remaining  portions  of  the  die.  Resistors  were 
probed  on  3  chips  to  determine  the  base  resistivity.  It  was  nearly  identical 
on  all  3  chips,  averaging  120  ohms  per  square. 

The  electrical  characteristics  of  all  devices  were  quite  uniform. 
The  breakdown  voltage  of  the  emitter-base  diode  was  6.b  volts  with  no  evi¬ 
dence  of  soft  breakdown,  instability,  etc.  At  77°K  the  breakdown  voltage 
decreased  to  5«b  volts,  while  the  forward  voltage  drop  at  2  raa  increased 

from  0.80  volts  to  l.lb  volts.  BVn„n  was  7  volts,  BV__.  was  2b  volts,  and 

ChO  CBO 

collector-substrate  breakdown  voltage  was  b2  volts.  BV.,  decreased  by 
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Fig.  2h  (cont.)  (c)  View  of  emitter-base  junctions.  Note 
differences  in  Junction-base  contact  distance,  and  metalli¬ 
zation  covering  only  part  of  some  contact  windows .  Bach 
small  division  is  2.5  um. 
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diffused  collector  process.  Dimensions  are  not  to  sc 


5  volts  at  77°K  while  the  collector-substrate  breakdown  decreased  by  6  volts. 

As  explained  earlier,  each  Junction  was  checked  on  a  curve  tracer 
first.  Then  the  device  was  connected  to  the  Cober  pulse  generator,  a  pulse 
width  selected,  and  the  device  pulsed  at  increasing  power  levels  until 
second  breakdown  occurred.  Typical  waveforms  are  shown  in  Figure  26  for 
reverse  biased  diodes  at  room  temperature.  Waveforms  at  -196°C  and  200°C 
were  similar. 

The  upper  waveform  in  each  case  is  the  voltage  across  the  Junction 
and  the  lower  waveform  is  the  current.  The  second  breakdown  transition  is 
readily  apparent  in  each  waveform.  In  parts  (e)  and  (f)  of  the  figure,  the 
current  is  off  scale  after  second  breakdown.  The  pulse  risetime  is  about 
20  nsec.  Current  and  voltage  are  roughly  constant  both  before  and  after 
second  breakdown.  Heating  effects  are  noticeable  during  the  first  part  of 
the  pulse  as  the  voltage  rises  and  the  current  falls.  Eventually  the  slope 
becomes  zero  and  then  changes  sign.  This  is  an  indication  that  second  break¬ 
down  is  imminent.  In  part  (d)  of  the  figure,  the  discontinuity  in  the  slope  is 
particularly  evident  in  the  voltage  waveform  about  1.5  psee  before  second 
breakdown.  This  "glitch"  was  noticeable  in  most  cases,  appearing  as  a  gradu¬ 
al  change  in  slope  in  the  shorter  pulses  (a),  (b),  (c),  and  as  a  "hesitancy" 
in  the  longer  pulses  (e)  and  (f).  If  the  pulse  was  terminated  after  the 
glitch  appeared  but  prior  to  second  breakdown,  no  damage  to  the  diode  could 
be  detected.  Subsequent  pulses  at  the  same  power  level  did  not  induce 
second  breakdown.  The  voltage  and  current  waveform  remained  the  same,  al¬ 
though  the  glit-'h  sometimes  shifted  to  a  slightly  earlier  time  in  the  pulse. 

The  interpretation  of  second  breakdown  by  Pudensrtein  et  al . ,  seems 
o  appl>  here.  It  is  reasonable  to  ascribe  the  discontinuity  in  the  voltage 


Fig.  2 6.  Pulse  waveforms  for  Part  A  in  the  reverse 
bias  mode,  (a)  Upper  waveform  10  v/cm,  lower  waveform 
200  ma/cm,  time  axis  100  nsec/cm.  (b)  5  v/cm,  200  ma/cm 
200  nsec/cm.  (c)  10  v/cm,  200  ma/cm,  500  nsec/cm. 
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Fig.  26  (cont.).  (d)  5  v/cm,  200  ma/cm, 
(e)  5  v/cm,  100  ma/cm,  20  usec/cm.  (f) 
200  usec/cm. 


2  usec/cm. 

5  v/cm,  100  ms, /cm. 
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to  a  change  in  the  conduction  process.  But  what  kind  of  a  change?  The  fact 
that  there  is  no  damage  means  that  it  is  reversible.  The  transition  to 
second  breakdovn  always  caused  damage  in  these  diodes,  hence  must  be  associ¬ 
ated  with  the  formation  of  a  microscopic  melt .  What  stage  lies  between 
normal  avalanche  breakdown  and  melt  formation?  The  only  conclusion  is  cur¬ 
rent  constriction  and  filamentation,  exactly  as  in  the  case  of  SOS  diodes. 


TABLE  II 

DIODE  BREAKDOWN  CALCULATIONS  FROM  FIGURE  26(c ) 


Before  Breakdown 

After  Breakdown 

Average 

V 

15 

5 

volts 

Average 

I 

300 

580 

ma 

Average 

P 

U.50 

2.90 

watts 

Time 

t 

1.87 

1.20 

usee 

Energy 

E 

8.ia  x  10"6 

3.1*8  x  10“6 

joules 

Resistance  R 

50 

8.6 

ohms 

Table  II  lists  the  quantities  of  interest  calculated  from  the 
waveforms  of  Figure  26(c).  The  average  power  is  taken  as  the  average  current 
times  the  average  voltage.  The  diode  resistance  is  the  ratio  of  average 
voltage  and  current.  This  is  important  when  related  to  the  generator  output 
impedance  (38  ohms  in  this  voltage  range).  As  the  diode  shifts  into  the 
second  breakdown  mode  its  resistance  changes  to  a  lower  value.  If  this  shift 
causes  the  diode  to  be  more  closely  matched  to  the  generator  impedance  the 
power  supplied  to  the  diode  will  increase.  For  test  device  A  the  diode  was 
nearly  matched  prior  to  second  breakdown.  After  second  breakdown  occurred 
there  was  a  mismatch  and  the  power  decreased.  This  is  an  important  considerar 


t ion  when  evaluating  damage  to  the  device.  ’ 

Of  primary  interest  is  the  relationship  between  the  delay  time 

a 

before  second  breakdown  and  pulse  power.  This  is  plotted  in  Figure  27.  The 
numbers  near  the  data  points  indicate  the  distance  in  microns  between  the 
emitter  junction  and  the  base  contact  window  (or  base  metallization  if  the 
window  is  not  entirely  covered).  The  dashed  lines  connect  points  from 
junctions  on  the  same  side  of  a  single  chip.  It  is  clear  that  junctions 
with  the  greatest  spacing  have  a  higher  threshold  power.  Unfortunately,  the 
distance  parameter  could  not  be  controlled  and  there  is  not  enough  data  for 
a  quantitative  relationship  with  the  other  variables. 

Seven  junctions  were  pulsed  to  failure  in  the  forward  direction. 

For  this  condition  second  breakdown  was  more  difficult  to  detect  from  the 
voltage-current  traces,  particularly  for  short  pulses.  The  junctions  were 
checked  on  a  curve  tracer  after  each  pulse  to  determine  when  permanent 
demage  took  place. 

Waveforms  of  three  Junctions  are  shown  in  Figure  28.  Part  (a) 
shows  two  successive  pulses  of  short  duration  superimposed.  After  the  first 
pulse  there  was  no  observable  damage  so  the  amplitude  was  increased  slightly. 
The  voltage  and  current  of  the  second  pulse  are  slightly  higher  until  nearly 
the  end  of  the  pulse  when  t,b<-  voltage  changes  slope  and  becomes  less  than  the 
voltage  of  the  first  pulse.  Aft'  r  the  second  pulse  both  electrical  and 
visual  damage  were  observed. 

Part  (b)  of  the  figure  shows  three  pulses  applied  to  the  same 

Junction.  Heating  becomes  more  and  more  pronounced  as  the  amplitude  is 

*  When  pulse  rise  time  and  transition  time  to  second  breakdown  are  not  negli¬ 
gible,  for  example  in  the  case  of  very  short  pulses,  the  pulse  width  is  taker 
between  75  per  cent  amplitude  points  on  the  voltage  waveform  prior  to  second 
breakdown . 
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Fig.  26.  Forward  bias  waveforms  for  Part  A. 
2  A/cm,  .100  nsec/cm.  (b)  5  v/cm,  500  ma/cm, 
(c)  5  v/cm,  500  ma/cm,  500  nsec/cm. 


(a)  10  v/cm, 
500  nsec/cm. 
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increased.  During  the  third,  pulse  the  voltage  reaches  a  maximum  value  then 
decreases.  No  damage  could  ‘be  observed.  Finally  the  pulse  in  part  ,(c)  of 
the  figure  vas  applied.  Near  the  end  of  the  pulse  the  voltage  shows  a 
greater  drop  and  the  rise  in  current  is  more  apparent.  This  pulse  caused 
permanent  damage. 

Part  (d)  shows  two  pulses  of  much  longer  d'lrution.  During  the 
first  pulse  the  voltage  rises  monotonically  while  the  current  remains  con¬ 
stant.  During  the  second  pulse  the  voltage  rises  rapidly,  reaches  a  broad 
maximum,  and  then  drops  abruptly  in  two  steps.  Permanent  damage  was  found . 

As  in  the  case  of  reverse  biased  Junctions  the  change  in  slope  of 
the  voltage  is  an  indication  that  second  breakdown  is  imminent.  However, 
damage  can  occur  without  the  precipitous  drop  in  voltage  which  reverse  biased 
Junctions  show.  Budenstein  et  al,,  noted  similar  behavior  in  their  diodes. 
Their  findings  indicate  that  the  current  distribution  and  heating  are  uniform 
up  to  and  including  the  change  in  slope.  Not  until  second  breakdown  is 
imminent  does  a  filament  form.  Since  the  advent  of  permanent  damage  is  not 
well  defined  in  the  short  pulses  of  Figure  28,  one  concludes  that  the  peak 
voltage  point  must  occur  at  a  very  high  temperature;  when  filamentation 
finally  does  take  place  the  formation  of  a  melted  column  does  not  require 
much  additional  temperature  rise  and  the  voltage  does  not  change  much.  For 
long  pulses  the  temperature  and  current  distributions  are  less  uniform  and 
the  transition  to  a  melted  column  would  be  more  noticeable. 

Junctions  were  also  "zapped"  at  liquid  nitrogen  temperature  and  at 
200°C,  about  10  Junctions  at  each  temperature.  Usually  three  J unci; ions  on 
the  same  side  of  a  chip  were  tested  at  the  three  ambient  temperatures  and 
with  the  same  pulse  width  to  keep  rs  many  variables  as  possible  the  same. 


Voltage-current  waveforms  were  much  the  same  at  each  temperature.  Threshold 
levels  were  shifted,  particularly  at  77°K  and  at  longer  pulse  lengths.  At 
pulse  widths  less  than  1  usee  the  change  in  threshold  power  was  small,  es¬ 
pecially  when  the  ambient  temperature  was  changed  to  200°C.  This  is  further 
evidence  that  filamentation  and  second  breakdown  occurs  at  very  high  tempera¬ 
tures. 

The  results  of  varying  the  ambient  temperature  are  plotted  in 
Figure  29.  The  curves  are  drawn  for  an  average  Junction-base  contact  width 
of  10  microns.  Data  points  are  omitted  to  avoid  confusion  except  for  the 
forward  biased  cases  where  all  seven  Junctions  had  nearly  identical  geome¬ 
tries  . 

«k  1 

The  three  reverse  bias  curves  have  not  yet  reached  a  t  asymptote 
at  the  shortest  pulse  width.  This  is  because  the  devices  are  small,  and  in 

i 

the  reverse  bias  mode  the  power  dissipation  is  highly  localized}  the  heat 
conducted  away  is  an  appreciable  fraction  of  the  total  energy  in  the  small 
heat-producing  volume,  so  adiabatic  conditions  are  not  approximated  with  100 
nsec  pulses.  In  the  forward  bias  mode  the  heat  is  produced  throughout  a  much 
larger  volume.  The  amount  of  heat  which  is  conducted  away  in  100  nsec  is  a 

smaller  fraction  of  the  total  heat  dissipated  and  the  situation  is  more 

nearly  adiabatic. 

One  circ  it  was  used  to  investigate  the  effect  of  heat  sinking. 

The  device  me  decapped  and  half  the  Junctions  were  zapped  at  different  pulse 
widths.  Then  a  drop  of  epoxy  was  placed  on  the  chip.  The  difference  in 

threshold  power  for  pulses  up  to  5  usee  long  was  less  than  10  per  cent.  This 

means  that  the  silicon  dioxide  layer  is  a  relatively  good  thermal  barrier  for 
these  pulse  times. 
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Damage  always  occurred  during  second  breakdown  in  reverse  biases 
diodes,  and  in  forward  biased  diodes  whenever  an  easily  discerned  second 
breakdown  transition  occurred.  The  most  sensitive  indicator  of  damage  was  t 
change  in  the  reverse  breakdown  characteristic  of  a  diode.  The  forward 
characteristics  were  unaffected  unless  massive  damage  had  occurred,  in  which 
case  the  diode  appeared  to  be  shunted  by  a  resistor.  The  damage  typically 
appeared  as  a  change  in  the  knee  of  the  curve  at  breakdown,  with  the  rest  of 
the  characteristic  unaffected.  The  characteristics  of  severed  damaged 
junctions  are  shown  in  Figure  30. 

It  should  be  pointed  out  that  a  criterion  of  minimum  detectable 
change  is  being  used  here.  Damage  does  not  necessarily  imply  a  malfunction 
in  circuit  operation.  The  diodes  in  Figure  30  would  perform  satisfactorily 
in  most  circuit  applications.  The  criterion  of  detectable  damage  however  is 
more  precise,  independent  of  application,  and  is  a  limiting  case. 

Visual  damage  correlates  roughly  with  the  time  spent  in  second 
breakdown  and  the  energy  dissipated  during  that  porticn  of  the  pulse. 

Junction  damage  was  visible  in  all  but  two  cases.  In  these  two  cases  the 
time  spent  in  second  breakdown  was  short  -  less  than  20  nsec  -  but  the  stress 
marks  were  revealed  when  the  metallization  and  oxide  were  removed  and  a  dash 
etch  applied  for  a  few  seconds.  It  should  also  be  noted  that  visible  damage 
might  occur  after  the  pulse  has  been  terminated,  for  instance  when  heat  from 
a  high  temperature  region  diffuses  to  a  contact  and  causes  delayed  melting 
or  contact  damage. 

Examples  of  damage  in  the  reverse  bias  mode  are  shown  in  Figure  31. 
One  of  the  junctions  in  part  (a)  shows  aluminum  in  the  melt  track,  with 
damage  to  the  metallization  on  both  emitter  and  base.  Visible  aluminum  or 
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Fig.  30.  Reverse  bias  breakdown  characteristics 
showing 'damage  due  to  second  breakdown. 
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Fig.  31.  Damage  due  to  reverse  bias  pulses.  In  each  case 
damage  extends  from  the  junction  to  the  base  contact. 
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'•white  streak"  damage  was  rarely  seen  on  these  devices.  The  other  zap  marks 
in  part  (a)  are  typical.  Usually  there  is  a  pronounced  mark  at  the  junction 
with  a  tail  extending  toward  the  hase  contact.  If  the  aluminum  does  not 
cover  the  base  window  entirely,  or  if  the  aluminum  is  removed,  the  track  is 
seen  to  terminate  at  the  base  contact.  There  is  visually  no  evidence  of  dam¬ 
age  extending  to  the  emitter  contact. 

In  part  (b)  of  Figure  31  the  junction  with  the  most  damage  was 
accidentally  pulsed  in  the  repetitive  mcde.  In  cases  where  this  occurred  the 
damaged  region  was  often  mushroom  shaped  with  the  cap  extending  laterally 
alon*.-  the  Junction  and  well  into  the  emitter  contact.  The  stem  extended  to 
the  bast  contact  where  metallization  damage  was  apparent.  The  other  stress 
marks  in  part  (b)  are  typical.  Parts  (c)  and  (d)  show  two  typical  zaps  at 
high  magnification.  In  both  cases  the  zaps  terminate  in  the  base  window,  but 
do  not  extend  to  the  emitter  contact. 

In  rare  instances,  damage  was  not  apparent  under  the  microscope. 

In  Figure  32  one  of  these  cases  is  shown.  In  part  (a)  the  metallization  has 
been  removed.  The  device  in  the  lower  left  has  a  very  faint  track  in  the 
base  region.  The  device  in  the  upper  right  has  no  visible  damage.  In 
part  (b)  of  the  figure  the  emitter  and  base  oxide  have  been  removed.  This 
enhances  somewhat  the  track  at  the  lower  left.  In  part  (c)  a  Dash  etch  has 
been  applied  for  a  few  seconds.  In  the  upper  right  of  the  photomicrograph  a 
faint  mark  is  now  seen  at  the  corner  of  the  Junction.  Thus  it  seems  certain 
that  second  breakdown  in  these  devices  is  always  accompanied  by  electrical 
damage  and  detectable  visual  damage. 

In  some  devices  the  zap  marks  extended  well  into  the  emitter  con¬ 
tact  area.  This  was  the  exception  rather  than  the  rule  and  usually  happened 


Fig.  32.  Search  for  junction  damage,  (a)  Metallization 
removal.  No  damage  visible  at  junction  at  upper  right, 
(b)  Oxide  etch.  Still  nc  visible  damage  at  upper  right. 


Fig.  32  (cont.).  (c)  Dash  etch.  A  faint  mark  now  appears 
at  Junction  at  upper  right.  A  faint  track  extending  to  the 
base  contact  could  be  discerned  under  the  microscope.  The 
other  stress  marks  clearly  extend  from  Junction  to  base 
contact . 


when  the  excinreiohinto  second  hreakdown  'had  a  long  duration »  Thus  it 
appears  that  in  reverse  biased’  Junctions  the  damage  first  takes  place  at  the 
Junction,  spreads  quickly  into  the  base  to  the  base  contact,  and  finally 
extends  through  the  emitter  to  the  emitter  contact, 

.  Six  forward  biased  Junctions  are  shown  in  Fig.  33.  In  the  first 
mlcrophotograph ,  part  (a)  of  the  figure,  several  interesting  features  appear. 
The  device  at  the  upper  left  appears  to  be  damaged  only- near  the  base  coctapt. 
The  device  at  the  upper  right  has  only  barely  disc«rnable  damage  The  damage 
to  the  diode  at  the  lower  left  is  normal,  while  the  device  on  the  lower  right 
seems  to  be  damaged  only  in  the  emitter  region. 

In  the  second  microphotograph,  part  (b)  of  the  figure,  the  metalli¬ 
ze  tion  has  been  removed.  This  reveals  some  details  not  apparent  in  part  (a). 
Finally  in  Fart  (c)  the  oxide  has  been  removed.  All  of  the  devices  show  at 
least  faint  tracks  extending  from  the  base  contact  to  the  Junction. 

In  the  device  at  the  lower  right  corner  the  damage  extends  all  the 
way  to  the  emitter  contact.  The  pulse  which  produced  this  damage  is  the  one 
shown  in  Fig.  28(d),  Referring  to  this  previous  figure,  note  the  final  short 
transition  which  occurs  after  second  breakdown  has  taker  place.  In  a  later 
part  of  this  report  it  will  be  shown  that  there  is  a  strong  correlation 
between  the  presence  of  the  small  voltage  drop  in  the  pulse  waveform  and 
damage  extending  from  the  Junction  to  the  emitter  contact. 
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Fig.  33.  Damage  to  forward  biased  junctions,  (a)  Photo¬ 
micrograph  of  the  junctions,  (b)  Metallization  removal. 
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Fig.  33  (coat.),  (c)  Oxide  removal .  All  Junctions  show 
dbma^e  extending  from  the  base  contact  to  the  Junction. 
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2.  Part  B 


This  circuit  vas  designed  for  a  previous  test  program.  A  number  of 
parts  were  stixl  available  and  have  been  utilized  in  the  present  study.  Sev¬ 
eral  components  were  fabricated  on  the  chip.  The  circuit  connections  for  the 
components  tested  in  this  program  are  given  in  Fig.  34.  The  diode  was  pulsed 
in  both  forward  and  reverse  biased  directions.  The  transistors  were  tested 
base-emitter,  collector-base,  and  emitter-base.  Unfortunately ,  it  was  not 
possible  to  remove  the  IK  resistor  across  the  emitter-base  terminals  without 
opening  the  pcakage  and  scratching  open  the  metallisation. 

The  chip  is  mounted  in  a  l4-pin  metal  flat  pack.  Photomicrographs 
of  the  die  are  given  in  Fig.  35*  Part  (a)  is  an  overall  view.  The  chip 
measures  44  x  58  x  8  mils.  Fabrication  is  by  the  planar  epitaxial 
method  with  a  buried  subcollector,  using  standard  10  MBs  DTL  technology. 

The  diode  is  xocated  at  the  edge  of  the  chip,  pinB  0  and  10  the 
eafcho&s  and  anode,  respectively.  There  is  a  metallization  stripe  exactly 
over  the  active  edge  of  the  Junction.  Three  transistors  are  accessible,  with 
emitters  connected  in  common  to  pin  1.  The  transistors  are  of  roughly  the 
same  dimensions  but  with  slight  differences  in  the  geometry.  The  differences 
con  be  seen  more  clearly  in  part  (b)  of  the  figure.  The  transistor  near  the 
middle  of  the  die  has  a  somewhat  greater  Junction  to  base  contact  distance 
(recall  that  this  made  a  significant  difference  in  the  threshold  power  of 
Device  A,  Fig.  27)  The  other  twc  transistors  have  emitter  windows  of 
slightly  different  area  and  location  within  the  emitter  region.  Details  o t 
the  diode  and  one  transistor  are  more  easily  seen  in  the  last  two  photographs 
in  the  figure.  The  small  divisions  are  6.25  urn. 

Other  derails  are  shown  in  Fig.  36.  The  epi  layer  is  10  ym  thick. 


90 


Fig,  31*,  Schem?,tJj  diagram 


Fig,  35.  Part  B,  special  test  circuit,  (a)  Overall 
view  of  chip,  (b)  View  of  transistors. 
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(d) 


Fig.  35  (cont,),  (c)  Single  transistor.  Each  small 
division  is  6. 25  urn.  (d)  Close-up  of  diode. 
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l—a  base  diffusion  is  2.4  van  deep,  with  a  sheet  resistivity  of  1J0  ohms  per 
square,  The  emitter  diffusion  is  1,7  ym  deep.  The  buried  layer  was  not 
clearly  revealed  by  lapping  and  etching  but  its  presence  is  evident  from  the 
photos  in  Fig.  35.  The  oxide  thicknesses  were  determined  to  be 

collector  -  8,500  A 

isolation  regions  -  6,800  t. 

base  -  4,900  $ 

emitter  -  1,500  t 

The  manufacturer's  specifications  for  the  circuit  are 


20-50 

bvcbo 

38-45 

volts 

^CEO 

13-15 

volts 

bvebo 

6. 3-6. 7  volts 

Electrical  characteristics  are  shown  in  Fig.  37 .  Part  (a)  is  the 
'.lode  characteristic.  Part  (b)  is  the  emitter-base  characteristic  of  the 
three  transistors  superimposed.  In  the  forward  direction  two  of  the  trans¬ 
istors  have  identical  characteristics.  The  third  transistor  (pins  1  and  5) 
has  an  increased  voltage  drop  due  to  the  increased  base  width.  In  the  re¬ 
verse  direction  the  shunting  effect  of  the  IK  resistor  is  apparent.  Again 
two  of  the  transistors  have  identical  traces  while  the  wide  base  device  has 
an  additional  voltage  drop. 

Port  (c)  of  the  figure  shows  the  collector-base  characteristics. 
The  switchback  in  the  reverse  bias  direction  is  due  to  the  emitter  being 
connected  to  the  base.  Part  (d)  is  the  collector-emitter  characteristic. 
The  forward  direction  ^  the  figure  is  under  normal  biasing,  that  is,  the 


95 


l~  Fig,  37.  Electrical  characteristics.  Part  B.  (a)  Diode 

\  characteristic.  Forward,  0.5  ma/cm  vertical  and  0.2  v/ cm 

i  horizontal.  Reverse,  0.5  ma/cm,  2  v/cm.  (b)  Emitter-bass 

|  Junction  characteristics.  Forward,  10  ma/cm  vertical, 

f:  0.2  v/cm  horizontal.  Reverse,  10  ma/cm,  2  v/cm.  See  text 

|  for  explanation  of  multiple  traces. 
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Fig.  37  (cont.).  (c)  Collector-base  Junction  characteristics. 
Forward,  2  ma/cm,  0.2  v/cm.  Reverse,  5  ma/cm,  10  v/cm. 

(d)  Collector-emitter  characteristics.  Forward,  5  ma/cm, 

10  v/cm,  collector  positive  with  respect  to  emitter.  Reverse, 
5  ma/cm,  1  v/cm,  inverted  polarity. 
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collector  is  positive  with  respect  to  the  emitter.  The  switchback  is  more 
complex  than  usual*;,,  again  due to*;the -.conxieet'i'on  between  emitter  and  base.  In 
the  reverse  direction  the  collector  is  negative  with  respect  to  the  emitter. 
The  collector-base  function  is  -then  foihratrdfblasedi  and  -the  IK  resistor 
appears  in  series  with  this,  junction..  When, the  voltage  is  large  enough  to 
cause  emitter-base  breakdown  the  current  begins  to  increase*  with  a  switch- 
back  occurring  due:  to  carrier  injection  from  the  collector  junction. 

The  emitter-base  resistor  wae  disconnected  from  one  transistor  by 
scratching  open  a  portion  of  the  metallization.  The  transistor  character¬ 
istics  were  then  found  to  he  normal,  with  parameters  very  close  to  those 
specified  by  the  manufacturer. 

Second  breakdown  data  on  Device  B  will  be  presented  in  the  follow¬ 
ing  order.  First  the  diode  results  will  he  given.  This  will  be  followed  by 
the  transistor  emitter-base  data,  then  the  collector-hase  results.  Finally 
the  collector-emitter  results  will  conclude  this  portion, 
a)  Part  B  -  Diode  Results 

These  diodes  differed  from  those  of  Device  A  mainly  by  having  a 
significantly  larger  area.  The  linear  dimensions  were  very  nearly  2*s  times 
those  of  Device  A  (except  for  diffusion  depths,  which  were  nearly  identical). 
The  threshold  power  was  consequently  much  greater.  The  base  sheet  resisti¬ 
vity  was  170  ohms  per  square,  compared  to  130  ohms  per  square. 

About  30  diodes  were  tested  -  8  forward  biased,  2  reverse  biased  at 
77°K,  and  the  rest  reverse  biased  at  room  temperature.  Representative  wave¬ 
forms  of  the  latter  testa  are  shown  in  Fig.  38,  parts  (a)  and  (b).  The  wave¬ 
forms  are  similar  to  those  of  the  eraitter-ba.se  junctions  of  Device  A.  There 
are  no  '‘glitches'*  in  Fig.  38,  however,  nor  did  they  ever  occur  in  these 
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Fig.  38.  Reverse  bias  waveforms  for  the  diode  of  Part  B. 

(a)  Upper  trace  20  v/cm,  lower  trace  500  ma/cm,  time  scale 
200  nsec/cm.  (b)  10  v/cm,  100  ma/cm,  100  ysec/cm.  (c)  10  v/cm, 
100  ma/cm,  100  usec/cm,  A  68  ohm  resistor  was  connected  in 
series  and  no  damage  resulted. 
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diodes.  In  most  cases  there  was  a  gradual  reversal  of  slope  prior  to  second 
breakdown.  The  transition  to  second  breakdown  was  usually,  but  not  always, 
abrupt.  In  one  case  the  slope  change  and  second  breakdown  drop  in  voltage 
was  so  gradual  that  the  waveform  was  more  like  that  of  a  forward  biased 
junction.  In  another  case  degradation  of  the  reverse  characteristic  was 
noted  even  though  no  second  breakdown  transition  was  observed.  In  four  cases 
second  breakdown  was  not  accompanied  by  observable  damage,  either  electrical 
or  visual.  Figure  38(c)  shows  one  such  case.  Note  tnat  the  drop  in  voltage 
is  significantly  less  than  in  the  other  two  photographs,  about  30/5  of  the 
maximum  voltage  as  opposed  to  60%.  In  all  such  events  the  pulse  width  was 
relatively  long.  A  68  ohm  resistor  was  connected  in  series  with  the  diode 
in  3  of  the  4  cases  in  order  to  obtain  finer  adjustment  of  the  pulse  ampli¬ 
tude. 

In  spite  of  the  diversity  in  behavior  the  threshold  power  points 
of  all  diodes  formed  a  smooth  curve  over  nearly  the  whole  range.  Only  at 
short  pulse  widths  was  there  significant  scatter  of  the  date  points. 

Forward  bias  waveforms  were  similar  to  those  of  Device  A.  In 
general,  there  was  a  gradual  transition  to  seoond  breakdown.  A  drop  in  volt¬ 
age  did  not  necessarily  signal  the  advent  of  permanent  damage.  Some  of  the 
waveforms  are  given  in  Fig,  39.  Note  in  particular  the  waveforms  in  part 
(b).  Only  the  third,  highest  amplitude  pulse  produced  permanent  damage. 

The  waveform  in  (c)  shows  a  small  but  noticeable  drop  after  the  large  change 
in  voltage.  In  a  later  part  this  will  be  shown  to  be  related  to  damage  pro¬ 
gressing  from  the  junction  through  the  material  on  the  low  resistivity  side 
to  the  contact. 

Power  versus  pulse  duration  curves  are  drawn  in  Fig.  40.  The 
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Fig.  39.  Forward  bias  waveforms  for  the  diode  of 
(a)  20  v/cm,  2  A/cm,  200  nsec/cm.  (b)  10  v/cm,  1 
500  nsec/cm.  (c)  5  v/cm,  500  raa/cm,  ?  usec/cm. 


Part  B, 
A/cm, 
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Fig,  Uo.  Threshold  power  vs  time  curves 
for  the  diode  of  Part  B. 


special  features  of  the  curves  ar«  isa/ie  apparent  by  comparing  *;lth  Pig.  27. 
The  most  obvious  difference  is  the  large  increase  in  threshold  power.  This 
is  due  to  the  much  larger  size,  toother  difference  Li  in  the  more  gradual 
transition  from  the  t*1  asymptote  to  the  t°  asymptote.  This  again  is  ^ 
function  of  the  size,  toother  interesting  feature  is  that  “he  threshold 
power  is  proportional  to  device  area  at  short  pulse  widths  for  both  forward 
and  reverse  biased  Junctions,  while  at  long  pulse  times  the  power  is  pro¬ 
portional  to  linear  dimensions. 

This  can  he  explained  cjs  follows.  At  short  pulse  widths  the  diode 
current  is  high.  About  half  the  power  is  dissipated  at  the  Junction  and  the 
rest  in  the  base  region  between  the  active  portion  of  the  emitter  Junction 
and  the  base  contact  (this  region  has  a  resistance  of  about  35  ohms).  At 
large  pulse  widths  on  the  other  hand,  nearly  all  of  the  power  is  dissipated 
in  the  thin  depiction  layer  of  the  Junction.  This  contrasts  with  the  BOS 
diodes  of  Sunshine  where  nearly  all  of  the  excess  voltage  above  the  -breakdown 
voltage  appeared  across  the  space  charge  resistance  of  the  Junction. 

Visual  damage  was  apparent  in  moBt  cases  >  with  the  exceptions 
already  noted.  Damage  usually  occurred  in  the  base  region  but  not  in  the 
emitter  region.  Damage  often  occurred  to  the  base  metallization  and  occa¬ 
sionally  the  aluminum  stripe  over  the  Junction  was  affected.  On  rare 
occasions  a  suort  appeared  between  the  stripe  and  one  of  the  diode  terminals . 

Two  damaged  Junctions  are  Bhevn  in  Fig.  4l.  Part  (a)  shows  typical 
damage  from  a  reverse  bias  test.  Damage  fextends  from  the  base  contact  to 
the  edge  of  the  aluminum  stripe.  When  the  aluminum  is  removed  in  such  cases 
the  track  can  be  seen  tc  extend  to  the  Junction.  In  part  (b)  the  stripe  has 
melted.  The  base  metallization  also  shows  damage.  There  is  no  damage  in  the 
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emitter  region.  The  pulse  waveforms  for  this  particular  diode  are  those 
shown  in  Fig.  38  part  (b).  The  pulse  width  up  to  second  breakdown  is  520 
ysec  and  the  average  pulse  power  is  8.4  watts.  After  second  breakdown  has 
occurred  the  diode  dissipates  about  6.7  watts  for  380  ysec.  In  is  difficult 
to  ascertain  exactly  when  the  aluminum  melted.  Clearly  the  damage  must  have 
occurred  after  filament  formation.  However,  it  cannot  be  ruled  out  that 
melting  may  have  occurred  after  the  pulse  terminated.  Unfortunately,  the 
one-dimensional  heat  flow  calculations  of  Section  IX  provide  little  insight 
for  events  with  such  long  pulses. 

To  summarize  the  results  of  the  diode  data  for  Device  B,  we  can 
state  that  the  differences  in  behavior  from  Device  A  can  be  attributed  to  the 
difference  in  size.  Anomalous  behavior  of  a  few  devices  was  observed, 
b)  Part  B  -  Emitter-Base  Junctions 

The  emitter-base  junctions  were  stressed  only  in  the  reverse 
direction  and  only  at  room  temperature.  Differences  in  Junction  character¬ 
istics,  breakdown  waveforms,  and  damage  can  be  correlated  with  the  differ¬ 
ences  in  geometry  of  the  transistors.  As  would  be  expected,  there  is  a  close 
relationship  between  the  second  breakdown  characteristics  of  the  emitter-base 
Junctions  and  the  diodes  previously  discussed. 

Current- volt age  waveforms  fell  into  distinct  categories,  some 
occurring  more  commonly  than  others.  The  most  usual  •waveform  showed  a 
slight  rise  in  voltage  to  a  broad  maximum,  a  slight  decrease,  and  then  an 
abrupt  drop  to  a  low  sustaining  level.  In  one  variation  the  transition  was 
first  to  an  intermediate  level  for  a  short  period.  T.j.s  happened  only  on 
rare  occasions  and  always  at  short  pulse  widths.  Degradation  of  the  Junction 
characteristic  did  not  always  occur  in  buch  cases.  In  a  second  variation  a 
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small  drop  in  voltage  (about  2  volts)  took  place  after  the  main  transition. 
This  happened  only  with  longer  pulses, 

Waveforms  are  shown  in  Fig.  h2.  Parts  (a)  through  (d)  of  the  fig¬ 
ure  ere  devices  on  the  same  chip.  Part  (a)  is  the  diode.  Part  (b)  is  the 
"normal”  transistor ,  pins  1-13.  It  differs  from  the  diode  waveforms  only  in 
having  a  lower  voltage.  The  differences  at  the  leading  edge  are  due  to  the 
pulse  generator  and  can  be  ignored.  Part  (c)  of  the  figure  is  the  waveform 
of  translator  1-3,  the  "wida-emitter"  transistor.  The  voltage  and  current 
are  slightly  reduced  because  the  pulse  generator  output  was  reduced.  The 
time  delay  before  second  breakdown  is  correspondingly  increased,  but  other¬ 
wise  the  waveforms  are  the  same  as  (b).  Part  (d)  of  the  figure  is  the  "vi de¬ 
base"  transistor,  pins  1-5.  This  is  an  example  of  the  variation  with  the 
transition  to  an  intermediate  level.  The  Junction  was  damaged  in  this  case. 
One  is  tempted  to  postulate  that  the  first  transition  corresponds  to  a  fila¬ 
ment  forming  in  the  base  region,  followed  by  elongation  to  the  Junction.  The 
evidence  is  too  sparse  to  confirm  this  contention,  however.  Furthermore, 
Budenstein  et  al . ,  found  this  sequence  of  events  occurring  only  with  long 
pulseB  in  their  diodes. 

Figure  k2  (e)  and  (f)  show  breakdown  waveforms  with  longer  pulses. 
The  waveform©  in  (e)  are  be^  jen  pins  1  and  3,  while  those  in  (f)  are  between 
pins  .1  and  5.  Note  the  small  drop  in  voltage  after  second  breakdown  lias 
occurred  in  (f).  This  is  an  illustration  of  the  second  -variation  in  waveform. 
Examination  under  a  microscope  revealed  a  stress  mark  extending  all  the  way 
from  the  base  contact  to  the  emitter  contact  of  the  transistor  In  (f),  but 
only  from  the  base  contact  to  the  Junction  in  (e).  The  curve  tracer  showed 
emitter-base  damage  ir  both  transistors  and  collector-eaittei  degradation  of 
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Fig.  k2.  Emitter-base  pulse  waveforms  for  Part  B,  reverse  bias, 
(a)  20  v/cm,  500  ma/cm,  100  nsec/cm.  (b)  10  v/cm,  500  ma/cm, 

100  nsec/cm.  (c)  10  v/cm,  500  ma/cm,  100  nsec/cra. 


107 


Fig.  ^2  (cont.).  (d)  20  v/cm,  500  ma/cm,  100  nsec/cm. 

(e)  10  v/cm,  200  ma/cm,  1  ysec/cm.  (f)  20  v/cm,  200  ma/cm, 
1  ysec/cm. 


the  transistor  in  (f).  Normally,  second  breakdown  involves  damage  to  the 
Junction  and  the  base  region  only.  If  the  duration  in  second  breakdown  is 
long  enough,  damage  extends  all  the  way  through  the  emitter  to  the  emitter 
window,  accompanied  by  a  further  drop  in  voltage.  Nearly  always  when  this 
variation  wa3  observed  the  distinction  in  visible  and  electrical  damage  could 
be  made. 

The  waveforms  show  some  correlation  with  the  spacing  from  junction 
to  emitter  contact  on  each  of  the  transistors.  The  spacing  wr»s  largest  be¬ 
tween  pins  1  and  3,  and  this  transistor  consistently  showed  the  largest  sus¬ 
taining  voltage  after  second  breakdown.  It  was  higher  by  at  least  2  volts. 
This  transistor  was  also  the  least  likely  to  show  a  second  transition  after 
second  breakdown.  The  explanation  is  straightforward.  Once  second  breakdown 
has  occurred  and  a  melt  column  bridges  the  base  region,  the  emitter  material 
serves  as  a  ballast  resistor.  The  greater  the  spacing  between  the  junction 
and  the  emitter  contact,  the  greater  the  resistance  and  the  greater  the  sus¬ 
taining  voltage.  If  the  distance  is  short  and  the  pulse  length  is  long  the 
filament  is  likely  to  bridge  the  emitter  region,  causing  the  three  observed 
effects:  a  drop  in  voltage,  a  stress  mark  tiidging  the  contacts,  and  changes 
in  the  collector-emitter  breakdown  characteristic . 

The  effects  of  .acreased  base  width  are  evident  prior  to  second 
breakdown  as  an  increase  in  voltage  required  to  induce  second  breakdown  with 
e  given  pulse  width.  Figure  U2  again  illustrates  this  when  part  (d)  is  com¬ 
pared  with  either  (c)  or  (b).  The  current  is  nearly  the  same,  but  the  volt¬ 
age  in  (d )  is  about  50#  greater. 

'’’he  diode  results  verify  thsse  observations.  The  diode  has  about 
the  same  amount  of  active  emitter  periphery  but  larger  emitter  and  base 
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widths  from  Junction  to  contact.  We  would  therefore  expect  a  higher  voltage 
before  second  breakdown  and  a  higher  sustaining  voltage  after  second  break¬ 
down.  Both  these  conclusions  are  substantiated  in  Fig.  42(a). 

In  Fig.  43  the  threshold  power  versus  pulse  duration  curves  are 
plotted.  The  upper  curve  is  for  the  wide-base  transistor  and  makes  clear  the 
greater  power  requirements,  especially  at  short  pulse  widths.  Data  points 
for  the  other  two  transistors  are  plotted  separately  but  there  is  no  clear 
effect  due  to  the  differences  in  emitter  width  and  only  one  curve  is  drawn. 
Both  curves  are  substantially  lower  than  those  of  the  diode,  Fig.  40,  where 
the  base  width  was  considerably  greater.  Recall  also  that  the  Junction  to 
base  contact  width  had  a  pronounced  effect  on  threshold  power  for  Fart  A 
(Fig.  27). 

A  photomicrograph  of  a  typical  emitter-base  zep  is  shown  in  Fig. 
44(a).  Junction  damage  is  clearly  seen,  with  the  track  extending  through  the 
base  to  the  base  contact.  This  type  of  damage  was  observed  in  about  75#  of 
the  cases.  The  breakdown  waveform  for  this  Junction  is  that  shown  in  Fig. 
42(e).  Part  (b)  of  Fig.  44  shows  the  damage  produced  by  the  pulse  of  Fig. 
42(f).  Damage  e?-;ten&8  all  the  .ray  from  contact  to  contact. 

It  is  clear  that  the  emitter-base  pulsing  has  yielded  a  richer 
variety  of  results  than  the  diode  tests.  In  the  first  place  the  important 
role  of  the  base  material  in  increasing  the  power  handling  capability  has 
been  corroborated.  The  effect  of  the  emitter  material  has  also  been  noted. 
Apparently  the  only  appreciable  modification  is  on  post-breakdown  behavior, 
where  a  higher  sustaining  voltage  results.  More  importantly,  increased 
emitter  width  has  the  desirable  effect  of  decreasing  the  transistor’s  sus¬ 
ceptibility  to  emitter  to  base  shorts  during  second  breakdown. 


no 


BASE  JUNCTIONS 


Fig.  M.  Emitter-base  Junction  damage,  (a)  Damage  produced 
by  the  waveform  of  Fig.  U2(e).  (b)  Damage  from  the  pulse 

in  Fig.  42(f). 
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Unfortunately,  the  changes  in  transistor  beta  caused  by  emitter- 
base  second  breakdown  were  net  investigated  in  this  study.  One  is  led  to 
speculate  that  the  effects  would  be  minimal  since  the  junction  forward 
characteristics  are  not  affected  and  because  the  damage  occurs  at  the  sur¬ 
face.  Others  have  found  major  changes  in  oeta; 30,31  however,  the  tests  were 
performed  on  discrete  devices .  Power  levels  were  higher  and  damage  tended 
to  be  more  massive. 

There  is  little  change  in  the  collector-emitter  breakdown 
characteristics  of  either  polarity  as  long  as  damage  is  confined  to  the 
emitter  junction  and  the  base  region.  If  the  damage  bridges  the  emitter- 
base  contacts  noticeable  changes  are  introduced  especially  in  the  inverted 
direction  with  the  collector  negatively  biased  with  respect  to  the  base, 
c)  Part  B  -  Collector-Base  Junctions 

Collector-base  junctions  were  tested  on  lU  parts,  10  in  the  reverse 
bias  direction,  2  in  the  forward  direction,  and  2  in  the  reverse  bias 
direction  at  liquid  nitrogen  temperature.  Good,  consistent  data  were  obtained 
as  far  as  pulse  waveforms  and  threshold  power  were  concerned.  Device  damage 
showed  some  variability.  As  in  previous  tests  the  occurrence  of  second 
breakdown  resulted  in  degradation  of  device  characteristics  in  nearly  every 
instance.  Unlike  previous  tests  the  base  contact  to  junction  distance  had 
only  a  barely  diecernable  effect  on  threshold  energies. 

The  three  transistors  on  each  chip  had  virtually  identical  wave¬ 
forms  for  the  same  pulse  conditions.  Representative  waveforms  are  shown  in 
Fig.  U5.  (a),  (b),  and  (c)  are  reverse  bias  pulses.  For  pulse  durations 

greater  x'mn  1  ysec  the  terminal  voltage  remains  at  a  high  level  for  a 
fraction  of  a  nsec,  then  drops  to  a  slightly  lower  level  as  the  device 
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Fig.  ^5.  Waveforms  for  collector-base  Junctions,  Part  B. 

(a)  20  v/cm,  500  ma/cn,  200  r.sec/cm.  (b)  20  v/cm,  200  ma/cm, 
1  usec/cm.  (c)  10  v/cm,  100  ma/cm,  100  usec/cm. 


114 


«f  WfcSW  tW*t*WM* 


«•  v^iNWvSiW^Ms',  l^vw 


St  *.  «>«<•.  »  >•>"«>„  /. 


« yw  *> -  ,-x pi.-*WV';;>  -  wy  -*•■» 


Fig.  I45  (cont.).  (d)  Forward  bias  mode,  5  v/cm,  500  ma/cm, 
200  ysec/cm.  Permanent  damage  occurs  at  the  small  drop 
in  voltage. 
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begins  to  oscillate.  Whether  or  not  oscillation  takes  place,  the  transistor 
is  damaged  only  when  there  is  a  pronounced  drop  in  voltage  to  about  the  50% 
level. 

With  short  forward  biased  pul  .a  the  voltage  waveform  was  much  like 
that  of  earlier  forward  biased  Junctions .  The  voltage  increased  to  a  broad 
maximum  and  then  began  to  decrease.  Transistors  had  to  be  checked  on  a  curve 
tracer  after  each  pulse  to  ascertain  when  damage  occurred.  Figure  45  (d) 
shows  a  forward  bias  pulse  for  a  long  pulse  width.  The  waveform  is  different 
from  previous  cases.  The  broad  maximum  is  still  present,  but  the  second 
breakdown  transition  is  unusually  small.  It  must  be  labelled  as  second 
breakdown,  however,  since  damage  occurred  in  all  transistors  tested  at  long 
pulse  widths  if  and  only  if  the  transition  was  present. 

Visible  damage  fell  into  two  categories.  Illustrations  appear  in 
Fig.  46.  Part  (a)  shows  the  damage  produced  by  the  pulse  in  Fig.  45  (b). 

For  pulse  widths  less  than  10  ysec  the  damage  was  of  this  type  and  looked 
like  the  damage  seen  in  a  forward  or  reverse  biased  emitter-base  Junction. 
Figure  46  (b)  shows  the  other  type  of  damage  commonly  seen  with  long  pulses . 

A  corner  of  the  base  metallization  shows  clear  evidence  of  heating.  The 
corner  of  the  emitter  metallization  also  shows  signs  of  reaching  high  tempera¬ 
tures.  When  this  type  of  damage  occurred  subsequent  pulses  increased  the 
amount  of  damage  but  no  stress  marks  of  the  type  in  (a)  appeared. 

Electrical  darns*.  *  followed  from  second  breakdown  in  all  but  4  cases. 
All  of  these  were  for  long  pulse  widths  and  3  had  a  68  ohm  resistor  connected 
in  series  with  the  device.  Damage  consisted  of  changes  in  the  reverse  bias 
breakdown  characteristic.*  of  the  emitter  Junction,  the  collector  Junction,  or 
from  collector  to  emitte: ,  either  polarity.  There  was  a  correlation  between 
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Fig.  U6.  Damage 
at  short  pulses. 


due  to  collector-base  pulses,  (a)  Damage 
(b)  Contact  damage  at  long  pulse  widths. 
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electrical  and  visual  damage.  When  the  stress  mark  appeared  at  the  surface, 
as  in  Fig.  46  (a),  the  emitter  junction  was  always  degraded  and  sometimes  the 
collector-emitter  breakdown  in  the  inverted  direction  was  changed.  Sur¬ 
prisingly  enough  the  collector-hase  junction  characteristics  were  often  un¬ 
changed.  Whenever  the  damage  occurred  only  to  the  base  metallization  the 
collector  junction  showed  deterioration  as  well,  as  the  collector-emitter 
breakdown  in  the  normal  direction.  Often  visible  damage  of  either  type  was 
accompanied  by  damage  to  all  junctions . 

The  evidence  indicates  the  presence  of  two  coexisting  current 
carrying  paths.  One  is  from  the  collector  contact  down  to  the  buried  layer, 
aloog  the  buried  subcollector,  back  up  through  the  epi  layer,  across  the 
avalanching  collector  junction,  and  to  the  base  contact.  The  other  path  is 
along  the  surface  from  the  collector  contact  to  the  reverse  biased  collector 
Junction,  through  the  base  region,  into  the  emitter,  across  the  reverse 
biased  emitter  Junction,  and  through  the  base  region  to  the  base  contact. 

At  any  given  terminal  voltage,  each  path  will  conduct  different  amounts  of 
current.  The  heat  produced  along  each  path  will  depend  on  the  current  density 
and  the  electric  field  along  each  portion  of  the  path.  However,  the  paths 
are  close  enough  to  each  other  that  the  thermal  responses  overlap.  (The 
diffusion  length  for  heat  flow  for  a  1  ysec  pulse  is  about  10  um,  the  depth 
of  the  buried  layer. )  Apparently  there  is  a  considerable  amount  of  heat  pro¬ 
duced  near  the  base  contact  due  to  current  flow  along  each  path.  Whether  a 
filament  forms  first  at  the  collector  junction  Just  beneath  the  base  contact 
or  at  the  surface  near  the  emitter  junction  depends  on  the  temperatuie  dis¬ 
tribution.  This,  in  turn,  depends  on  the  current  division  between  the  two 
paths  for  a  particular  applied  voltage.  At  high  pulse  voltages  apparently 
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the  conditions  for  filament  formation  are  reached  first  between  the  base 
contact  and  the  emitter.  Wherever  the  breakdown  occurs,  the  two  junctions 

are  close  enough  near  the  base  contact  that  when  one  is  damaged  the  melted 

\ 

region  is,  very  likely  to  damage  the  other  also. 

Brown ,  Holder,  and  Ruwe30  describe  an  anomalous  mode  of  failure  in 
discrete  transistors  in  terms  of  a  surface  avalanche  which  they  refer  to  as 
"Holder  avalanch^/^  They  also  observed  emitter  Junction  damage,  both  elec¬ 
trical  and  visible,  when  pulsing  between  collector  and  base.  Their  explana¬ 
tion  is  along  the  lines  of  that  given  here  even  though  the  geometry  of  a 
discrete  transistor  is  somewhat  different. 

The  threshold  damage  curve  is  plotted  in  Fig.  1*7.  The  curve  is 
drawn  for  the  data  points  corresponding  to  the  '’normal"  and  "wi de-emitter" 
transistor  geometries,  which  appear  to  have  no  distinct  differences.  A 
separate  curve  is  not  drawn  for  the  "wide-base"  geometry  but  the  data  points 
are  seen  to  lie  slightly  above  the  previous  curve  over  the  entire  range. 

This  again  points  to  the  existence  of  a  path  that  carries  a  portion  of  the 
current  parallel  to  the  surface  at  the  base  contact;  otherwise,  the  increased 
spacing  would  have  no  effect. 

The  curve  shows  nearly  a  t”1  asymptote  at  short  pulses.  This  im¬ 
plies  heat  generation  ever  a  fairly  large  volume  -  despite  the  large  drop  at 
the  junction  -  such  that  the  heat  conducted  away  is  only  a  small  portion  of 
the  total  heat  produced. 

The  power  level  at  1  msec  is  the  same  as  the  emitter-base  functions, 

Fig.  1*3.  At  %  usee  the  power  is  greater  by  a  factor  of  about  2.  A  rather 

surprising  result  is  found  when  the  data  of  Fig.  1*7  is  compared  to  that  of 

the  diode,  Fig.  1*0.  The  results  almost  coincide  for  reverse,  forward,  and 
low  temperature  cases. 
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PART  B 

COLLECTOR -BASE  JUNCTIONS 
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d)  Part  B  -  Collector-emitter  Results 


Collector-emitter  pulsing  vould  "be  expected  to  be  more  complex  bs= 
cause  of  the  interacting  Junctions*  and  indeed  this  proved  to  be  the  case. 

At  first  glance  ve  vould  expect  greater  power  handling  capability  because  of 
the  greater  area  between  the  emitter  and  collector  contacts.  Again  this 
proved  to  be  true.  Pulse  waveforms  were  more  complicated  than  in  previous 
test  devices.  There  were  two  distinct  voltage  drops*  both  of  which  could 
produce  contact  damage  but  only  one  of  which  was  associated  \rith  electrical 
damage.  This  will  be  explained  in  greater  detail  later  on*  with  a  breakdown 
model  postulated  to  explain  the  results. 

Another  distinctive  feature  was  the  common  occurrence  of  demage  to 
metallization  runs  far  from  the  junction.  This  is  a  consequence  of  the  high 
currents  encountered  during  breakdown. 

Seven  parts  were  tested.  Some  underwent  the  normal  procedure  while 
others  were  used  in  multiple  or  repetitive  pulse  testing  in  a  closed-circuit 
TV  system.  All  transistors  were  stressed  with  the  collector  positive  with 
respect  to  the  emitter  and  at  room  temperature  only. 

In  the  initial  tests  a  significant  drop  in  voltage  was  noted  after 
a  certain  delay  time  which  depended  on  pulse  amplitude.  This  was  Interpreted 
as  second  breakdown.  Figure  48  shows  several  waveforms.  At  low  pulse  ampli¬ 
tudes  (parts  (a)  and  (b)  of  the  figure)  the  voltage  change  is  significant  and 
well  defined.  At  high  pulse  amplitudes  (parts  (c)  an"  (d)  of  the  figure) 
although  the  voltage  drop  is  the  same  in  amplitude  it  is  a  much  smaller 
percentage  of  the  toted  voltage.  Also  a -preceding,  much  smaller  step  in 
iroltage  can  be  seen. 

This  response  was  consistently  observed  in  every  transistor.  It 
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Fig,  U8.  Collector-emitter  tretkdown  waveforms  for 
transistors  of  device  B.  (a)  10  v/ cm,  200  ma/cm, 

10  ysec/cm,  (b)  10  v/cm,  500  ma/cm,  1  ysec/cm. 
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Fig.  1*8  (cont.).  (c)  10  v/cm,  500  ma/cm, 

500  ns.,  c/cm.  (d)  20  v/cm,  1  A/cm,  200  nsec/cm. 
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never  resulted  in  electrical  damage  -  all  terminal  pairs  showed  the  same 
characteristics  before  and  after  such  pulses.  Because  of  this  the  transis¬ 
tor  could  be  pulsed  repeatedly  into  "second  breakdown"  over  a  wide  range  of 
pulse  conditions  and  still  yield  exactly  the  same  waveforms.  No  differences 
could  be  detected  among  the  three  types  of  transistors  either  in  threshold 
energies  or  in  waveforms. 

The  origin  of  the  voltage  discontinuity  must  be  thermal  in  origin 
because  the  delay  time  depends  on  pulse  amplitude.  In  the  range  from  3  usee 
to  4  msec  the  threshold  pv.wer  versus  delay  time  plot  was  a  smooth  curve 
slightly  higher  than  the  collector-base  threshold  curve.  Figure  47,  at  long 
pulse  widths  but  greater  by  a  factor  of  2-1/2  at  short  pulse  widths.  Below 
1  nsec  the  curve  seemed  to  level  off  in  a  constant  power  dependence.  This 
is  not  certain,  however,  since  the  transition  is  poorly  defined  at  low  pulse 
widths  and  only  a  few  data  points  were  taken.  If  true,  it  would  point  to 
another  mechanism  becoming  dominant  such  as  current  induced  second  breakdown. 

Although  electrical  damage  was  not  observed  with  this  type  of 
"second  breakdown"  contact  damage  was  produced.  The  damage  consisted  of 
darkening  of  the  emitter  metallization  over  a  sizeable  portion  of  the  con¬ 
tact,  in  contrast  to  the  usual  damage  where  the  aluminum  melted  at  a  specific 
site.  To  determine  the  relationship  between  contact  darkening  and  the  change 
in  voltage  a  video  camera  was  mounted  over  a  microscope  and  connected  to  a 
video  tape  recorder  and  TV  monitor.  Pulses  were  applied  singly  or  at  a  low 
prf  (about  10  pps)  for  tens  of  seconds  at  a  time. 

The  conclusions  were  that  the  voltage  drop  was  not  absolutely 
necessary  for  contact  damage  but  that  damage  was  induced  at  a  distinctly 
accelerated  rate  when  the  voltage  change  was  present.  When  the  pulse 


duration  was  lese  than  the  delay  time ,  some  damage  was 'evident  after  a  suf¬ 
ficient  number  of  pulses  had  been  applied.  When  the  pulse  duration  was  long 
enough  to  include  the  transition  a -single  pulse  was  sufficient  -to  produce 
noticeable  damage.  .  - 

Figure  49(a)  shows  emitter  contact  damage  to  a  transistor  after 
two  pulses.  "Second" breakdown"  occurred  only  during  the  second  pulse.  Un¬ 
fortunately,  the  device  was  not  opened  for  examination  until  after  the 
second  pulse  so  the  extent  of  contact  damage  after  the  first  pulse,  if  any, 
is  not  known.  Notice  that  there  is  no  distinct  evidence  of  melting  and  that 
df image  appears  over  a  fairly  large  region.  Figure  49(b)  shows  the  appearance 
of  another  contact  after  5-10  minutes  of  low  prf  pulsing  at  sub-breakdown 
pulse  widths.  Although  the  contact  is  heavily  damaged,  there  is  still  no 
change  in  Junction  breakdown  characteristics  at  this  point.  Degradation  did 
not  occur  until  further  pulsing  eventually  caused  the  darkened  metal  to  spread 
out  and  reach  the  emitter-base  junction. 

Because  of  the  lack  of  junction  degradation  even  after  many  excur¬ 
sions  into  what  appeared  to  be  second  breakdown,  and  because  the  transition 
at  high  pulse  amplitudes  was  not  pronounced,  it  was  decided  to  pulse  the 
transistors  at  increasingly  higher  levels  until  permanent  junction  damage 
did  occur.  As  the  pulse  energy  was  increased,  either  by  increasing  the  pulse 
width  or  the  pulse  amplitude,  a  new  breakdown  phase  came  into  existence, 
heralded  by  a  further  reduction  in  voltage.  A  single  pulse  incorporating 
this  new  feature  was  sufficient  to  produce  dramatic  effects.  These  included 
permanent  device  degradation  as  well  as  visible,  often  massive,  contact 
damage. 

Examples  of  pulse  waveforms  under  these  conditions  are  shown  in 
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Fig.  *49.  Contact  damage  produced  by  pulses  at/or 
just  below  first  transition  levels.  No  electrical 
degradation  resulted,  (a)  After  2  pulses,  (b)  A 
different  transistor  after  5-10  minutes  of  low  prf 
pulsing. 
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Figure  50.  The  first  photo  is  for  very  long  pulse  lengths.  The  initial 
voltage  peak  and  two  transitions  are  completed  in  about  750  usee .  The 
residual  voltage  is  only  about  5  volts  and  has  several  minor  fluctuations. 

The  s^ond  photograph  in  the  figure  is  at  intermediate  pulse  lengths.  Here 
the  details  are  much  clearer.  The  voltage  initially  rises  from  device 
heating  to  a  maximum,  then  declines  by  5  volts  in  the  first  transition. 

The  voltage  continues  to  decline  slowly  at  this  intermediate  stage  until 
final  breakdown  takes  place  at  about  19  usee.  This  final  drop  is  compara¬ 
tively  gradual  and  occurs  in  two  stages.  The  residual  voltage  is  again 
about  5  volts. 

The  last  waveform  in  the  figure  is  representative  of  short  pulses. 
The  voltage  overshoots  at  the  beginning,  then  rises  again  due  to  heating. 

The  first  break  occurs  about  300  nsec  after  the  beginning  of  the  pulse.  This 
is  followed  by  a  gradual  decline  until  final  breakdown  takes  place  about 
800  nsec  from  the  start  of  the  pulse .  Just  before  the  pulse  terminates ,  the 
voltage  rises  sharply  -  this  is  associated  with  melting  of  the  metalliza¬ 
tion  run. 

Several  damaged  transistors  can  be  seen  in  the  photomicrographs  of 
Figure  51.  The  first  photo  shows  a  failure  induced  by  a  short  pulse.  Alumi¬ 
num  at  the  corner  of  the  contact  has  melted  and  the  stress  mark  extends  into 
the  emitter  region.  This  is  the  only  visible  damage.  There  is  no  evidence 
of  damage  to  the  emitter  Junction,  the  collector  Junction,  or  the  collector 
contact.  A  curve  tracer  did  show  junction  damage,  however. 

The  second  photomicrograph  is  the  same  transistor  shown  in  Fig. 
l*9(a).  Several  more  pulses  have  been  applied,  the  last  one  sufficient  to 
induce  the  second  voltage  transition. 
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Fig.  50.  Collector-emitter  breakdown  waveforms  showing 
the  second  transition,  (a)  5  v/cm,  500  ma/'cm,  1  msec/cm. 
(b)  10  v/cm,  500  ma/cm,  5  psec/cm.  (c)  20  v/cm,  1  A/cm, 
200  nsec/om. 
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Fig.  51  (cont.).  Additional  examples 
Note  the  damage  to  metallization  runs 


of  contact  damage, 
in  each  case. 


130 


Figure  51(c)  shows  damage  from  a  15  usee  pulse  terminated  just  after 


the  final  voltage  drop.  A  corner  of  the  metal  contact  has  melted  and  the 

track  ertends  all  the  way  to  the  collector.  Of  special  interest  is  the 

meltec.  metallization  stripe.  Melting  occurs  only  where  the  oxide  is  thickest 

and  probably  was  initiated  at  an  oxide  step.  On  this  particular  device  the 

o 

aluminum  thickness  is  8,000  A. 

The  last  transistor  in  the  figure  has  the  aluminum  completely 
melted  and  opened.  The  damage  was  produced  by  the  pulse  of  Fig.  50(c). 

Damage  can  be  seen  at  the  collector  and  emitter  contacts  and  in  the  emitter 
region.  Again,  there  is  no  sign  of  damage  to  the  junctions  in  the  picture. 

The  threshold  power  curve  for  collector-emitter  tests  is  plotted  in 
Fig.  52.  The  average  power  at  each  point  was  computed  graphically  from  the 
voltage  and  current  waveforms.  There  was  no  dependence  on  transistor  geometjy. 
In  some  events  the  device  failed  too  early  in  the  pulse  to  make  an  accurate 
estimate  of  the  delay  time.  In  such  cases  the  power  was  computed  from  the 
preceding  photograph  and  the  point  entered  as  "not  failed”  in  the  figure. 

The  threshold  power  is  greater  for  the  collector-emitter  mode  than  any  other, 
exceeding  collector-base  pulsing  by  50!?  at  1  msec  and  by  hOOJ?  at  1  nsec. 

In  interpreting  t.he  results  of  collector-emitter  pulsing  the 
pertinent  questions  are:  what  is  the  significance  of  each  of  the  two  volt¬ 
age  drops,  and  how  does  each  relate  to  device  damage?  The  first  transition, 
although  thermal  in  nature,  does  not  result  in  a  microscopic  melt;  if  it  did 
tnere  would  be  permanent  damage.  It  is  always  reproducible.  Contact  damage 
does  result  after  such  a  transition,  so  it  indicates  that  a  localized,  high 
temperature  site  has  formed.  The  mere  fact  that  the  voltage  is  discontinuous 
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Threshold  damage  curve  for  collector- 
emitter  pulsing.  Part  B. 


implies  a  change  in  the  conduction  mechanism 


The  hypothesis  is  that  the  voltage  drop  signals  the  formation  of  a 
current  constriction  at  the  collector  Junction  near  the  corner  of  the  emitter 
contact  and  the  formation  of  a  broad  filament  extending  from  the  constriction 
to  the  collector  contact.  The  local  temperature  rises,  but  not  to  the 
eutectic  point  of  the  silicon-aluminum  interface.  This  configuration  appears 
to  be  stable  for  short  periods  of  time  if  the  total  power  dissipation  is 
less  than  15  watts.  At  levels  less  than  this  the  voltage  remains  constant  or 
rises  almost  imperceptibly  for  periods  ranging  up  to  200  vsec.  When  the  power 
level  is  greater  than  15  watts  the  voltage  continues  to  fall  slowly.  After  a 
delay  period  the  process  reaches  the  second  phase. 

The  final  transition  introduces  permanent  damage.  This  means 
localized  melting  has  occurred.  The  proposition  is  that  the  filament  grows 
from  the  collector  junction  through  the  base  region,  across  the  forward 
biased  emitter  Junction  and  to  the  emitter  contact.  Then  a  melt  forms. 

The  extent  of  the  damage  would  depend  primarily  on  how  long  the 
pulse  lasted  after  breakdown  occurred  and  on  the  power  levels.  In  every  case 
the  damage  at  the  emitter  contact  vas  apparent.  Occasionally,  localized 
damage  appeared  at  both  collector  and  emitter  without  a  stress  mark  in  between 
(as  in  Fig.  51(d)).  This  is  an  indication  that  a  filament  extended  from  con¬ 
tact  to  contact.  If  the  oxide  and  aluminum  had  been  removed  on  these  devices 
it  seems  highly  likely  that  a  continuous  damage  streak  would  be  found. 

This  model  can  only  be  regarded  as  speculative  since  no  etching  or 
cross-sectioning  of  these  devices  was  done  to  determine  the  location  and 
extent  of  the  resolidified  melt.  Nor  is  there  any  way  to  determine  the 
sequence  of  steps  as  in  the  SOS  diode  experiments.  It  does  fit  the  observed 
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facts  for  this  transistor,  however,  and  is  a  reasonable  extension  of  second 
breakdown  in  simple  junctions.  Budenstein,  Pontius,  and  Smith  have  shown 
that  broad  filaments  can  exist  and  remain  stable  without  melting  under  cer¬ 
tain  situations.  In  the  results  of  the  next  device  to  be  presented  there  is 
evidence  that  transitions  similar  no  the  final  voltage  drop  are  associated 
with  damage  in  the  emitter  region.  Thus,  although  the  model  cannot  be  re¬ 
garded  as  conclusive,  there  is  good  evidence  to  support  it. 

Another  theory  that  has  been  proposed  is  that  the  final,  destructive 
phase  is  due  to  the  melting  of  aluminum  at  the  contact  which  then  alloys  into 
the  silicon  to  form  a  melting,  conducting  channel.  This  does  not  appear  to 
be  the  case  here.  There  is  evidence  that  the  voltage  transition  must  be 
essentially  complete  before  permanent  damage  is  encountered.  If  melting 
aluminum  initiated  the  final  phase,  some  form  of  damage  would  appear  even  if 
the  transition  was  not  complete. 

3.  Part  C 

This  circuit  is  a  standard  DTL  NAND  gate.  Photomicrographs  of  the 
circuit  appear  in  Fig.  53.  The  chip  is  32  mils  square  and  is  glassivated. 

The  input  diodes,  in  the  upper  right  and  left  corners  of  the  picture  are  the 
only  components  which  can  be  individually  tested.  Part  (b)  shows  details  of 
one  input  group.  Each  division  is  2.5  nm. 

The  h  diodes  are  formed  by  emitter  diffusions  into  two  separate 
base  regions .  Th  two  p-rogions  are  connected  by  an  overlapping  n+  diffus¬ 
ion,  and  an  aluminum  anode  stripe  makes  contact  to  the  n+  region  and  the 
p-type  regions  on  either  side.  A  diagram  of  this  construction  and  the  pin 
connections  are  given  in  Fig.  5^< 

Three  different  diode  configurations  were  tested,  all  at  room 
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Fig.  53.  Photomicrographs  of  the  circuit  of  Part  C. 
(a)  Each  division  is  12.5  pm.  (b)  Close-up  of  the 
diode  cluster.  Each  division  is  2.5  urn. 
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Fig.  54.  Pin  connections  for  the  input  diodes 
and  details  of  the  construction. 


temperature  and  all  in  either  the  reverse  bias  or  the  back-to-back  mode.  The 
standard  configuration  is  simply  an  ordinary  diode  connection  between  pin  8 
or  9  and  pin  7.  The  second  combination  is  two  diodes  on  opposite  sides  of 
the  anode  stripe  connected  back-to-back  (pins  1  and  3).  For  the  third  com¬ 
bination  two  diodes  in  the  same  p-region  (pins  2  and  3)  were  pulsed  back- 
tc-back. 

The  diode  cha.'icteristics  are  shown  in  fig.  55.  Part  (a)  is  the 
characteristic  between  p:.ns  9  and  7.  In  the  forward  direction  the  diode 
shows  signs  of  reaching  high  level  injection  when  the  current  is  only  15  ma. 

In  the  reverse  direction  a  discontinuity  is  seen  as  a  result  of  the  inter¬ 
action  between  the  emitter  and  collector  junctions  of  the  transistor  formed 
by  the  n+  diffusion,  the  p  diffusion,  and  the  n-type  is'-  ited  region.  In 
other  words  the  reverse  characteristic  shows  the  collector-emitter  breakdown 
of  the  transistor  operating  in  the  inverted  mode. 

Part  (b)  is  a  back-to-back  diode  characteristic.  Part  (c)  shows  a 
portion  of  the  reverse  breakdown  characteristic  on  an  expanded  scale  to  show 
the  difference  in  voltage  drops  across  different  pairs  of  diodes.  The 
characteristic  with  the  smallest  voltage  drop  is  for  two  diodes  in  the  same 
p-region,  while  the  other  trace  belongs  to  two  diodes  on  opposite  sides  of 
the  center  contact.  The  increased  voltage  drop  had  little  efi'ect  on  thresh¬ 
old  power  except  ac  short  pulse  lengths. 

Representative  se(  .id  breakdown  waveforms  are  shown  in  Fig.  56. 

The  pulse  in  the  first  photograph  was  applied  between  pins  1  and  2.  Weveforms 
for  this  connection  were  typically  of  this  form  -  a  relatively  small,  slow 
voltage  drop  often  in  two  step3,  with  a  large  residual  voltage.  Visible 
damage  was  not  often  appa'.  ent  in  this  mode  of  pulsing.  When  damage  could  be 
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Fig.  55.  Diode  characteristics,  (a)  Pins  9  and  7.  Forward; 

0.5  v/crn,  5  ma/cm.  Reverse;  2  v/ era,  5  ma/ern.  (b)  Pins  9  and  8. 
2  v/cra,  5  aia/cm.  (c)  C.5  v/ctn,  5  ma/cm.  Left-hand  characteris¬ 
tic,  pins  2-3.  Right-hand  characteristic,  pins  1-3. 


seen  it  appeared  at  the  contact  and  was  localized.  -The  .glass  layer  and 
aluminum  were  removed  on  one  part.  When  a  Dash  etch  was  applied  the  zap  mark 
could  be  easily  discerned  extending  downward  and  into  the  emitter  material . 

The  exact  extent  of  the  melt  column  could  not  be  determined.  However,  there 
was  never  any  indication  of  damage  extending  laterally  into  the  p-region 
except  on  one  part  which  failed  from  repetitive  pulsing. 

Microsecond-length  pulses  applied  to  any  of  the  three  diode  con¬ 
figurations  resulted  in  breakdown  waveforms  very  much  like  the  one  in  part 
(a)  of  the  figure.  Any  visible  damage  was  always  confined  to  the  n+  emitter 
region  in  such  cases. 

Longer  pulses  applied  to  pins  1  and  2  produced  the  same  type  of 
waveform  and  the  same  type  of  damage.  Longer  pulses  applied  to  the  other  two 
diode  configurations  produced  different  results.  The  waveform  in  part  (b), 
Fig.  56  is  typically  seen  between  pins  8  and  7.  The  first  transition  does 
not  produce  permanent  damage.  The  second  transition  does.  Furthermore, 
whenever  both  transitions  appear  the  zap  mark  extends  past  the  emitter 
junction,  through  the  p-regicn,  and  to  the  aluminum  contact.  In  other  words 
whenever  a  pulse  produces  permanent  damage  a  zap  mark  can  always  be  found  in 
the  emitter  material.  Whenever  damage  is  found  in  the  base  region,  both 
transitions  are  present  in  the  pulse  waveform.  It  appears  then  that  the 
first  transition  in  part  (b)  of  the  figure  is  caused  by  filament  formation 
in  the  p-region.  The  second  voltage  drop  comes  about  when  the  filament 
extends  through  the  emitter  region  to  the  emitter  contact  and  a  melt  forms. 
Thus  the  breakdown  and  damage  process  is  similar  to  that  of  the  transistors 
of  Part  B. 

The  third  photograph  in  Figure  56  is  the  waveform  between  pins  2 


and  3.  The  emitter  regions  are  lined  up  end  to  end  in  this  configuration 
and- the  power  levels  are  lowest.  Again  the  double  transition  is  associated 
with  damage  to  both  the  n+  and  the  p-region  while  a  single  transition  is 
associated  with  damage  in  only  the  emitter  region.  It  appears  possible  that 
the  order  of  events  could  be  reversed  in  part  (c),  but  there  is  no  way  of  , 
determining  this  from  the  data  which  was  taken. 

Figure  57  shows  the  threshold  damage  curves  for  this  circuit .  Due 
to  the  small  size  the  power  levels  are  low,  and  the  t°  asymptote  is  reached 
in  about  100  ysec.  The  lowest  curve  is  for  pulses  applied  between  pins  2  and 
3.  In  the  other  configurations  the  emitter  diffusions  are  broadside  to  each 
other  and  the  power  levels  are  50%  greater.  This  points  out  that  the  power 
required  to  damage  a  circuit  depends  on  the  length  of  the  active  portion  of 
the  emitter  junction  and  not  on  the  area  of  the  emitter  diffusion. 

The  figure  also  shows  that  the  threshold  power  is  the  same  for  con¬ 
nections  1-2  and  8-7  for  t  greater  than  10  ysec.  One  would  expect  at  least 
a  small  difference.  No  matter  whether  the  current  path  between  1  and  2  is 
along  the  surface  or  via  a  buried  subcollector  there  are  extra  parasitic 
losses  as  the  characteristics  of  Fig.  55(c)  show.  The  voltage  change  is  only  • 
about  h%  at  low  current  levels  and  if  it  is  about  the  same  for  low  power 
pulses,  it  would  be  difficult  to  see  in  the  normal  \0~15%  scatter  of  the  data 
points.  At  high  pulse  power  levels  the  increase  in  parasitic  losses  is 
readily  apparent.  The  damage  curve  rapidly  approaches  a  t_1  variation, 
indicating  power  dissipation  throughout  a  broad  region. 

In  most  cases  damage  to  the  devices  -  both  electrical  and  visible  - 
was  slight.  Electrical  damage  often  consisted  only  of  small  changes  in  the 
knee  of  the  breakdown  curve.  Long  pulses  caused  greater  damage  but  not  so 


much  as  to  give  the  appearance  of  a  low  resistance  short  across  the  junction. 
Zap  marks  were  similar  in  appearance  to  those  of  other  base-emitter  junctions 
except  that  damage  always  existed  in  the  emitter  region  and  sometimes  only  in 
the  emitter  region.  There  was  rarely  any  damage  to  the  metallization.  One 
reason  was  the  small  amount  of  heat-sinking  by  the  glass  layer  (about  15  pm 
thick).  A  more  important  reason  was  the  small  size  of  the  devices  and  the 
lower  total  power  levels  involved. 

The  upshot  of  all  this  is  to  point  out  the  differences  in  second 
breakdown  between  simple  p-n  junctions,  integrated  circuit  transistors,  and 
three-region  structures  such  as  those  discussed  here.  In  all  cases,  a  con¬ 
striction  forms  at  the  reverse  biased  junction,  with  the  exception  of 
Part  C  under  high  power  pulses,  constriction  is  followed  by  filament 
formation  in  the  high  resistivity  material.  In  a  simple  p-n  junction  the 
next  stage  is  melt  formation  and  permanent  damage.  In  other  devices  the 
filament  grows  all  the  way  to  the  emitter  contact  before  damage  occurs.  At 
high  power  levels  in  the  Part  C  structures  the  filament  forms  only  in  the 
emitter  region  before  failure  occurs. 

This  model  cannot  be  regarded  as  conclusive,  but  there  is  good 
correlation  among  the  various  devices  tested  in  spite  of  the  differences  in 
size  and  construction. 

it.  Part  D 

Part  D  is  a  dielectrically  isolated  dual  DTL  expandable  NAND  gate. 

U  diodes  with  a  common  anode  terminal  are  fabricated  on  each  side  of  the 
chip.  A  photomicrograph  of  two  of  the  input  diodes  is  shown  in  Fig.  58. 

The  diodes  are  fabricated  in  pairs  in  four  separate  isolated  n~ 
regions.  All  diodes  have  the  same  dimensions  but  each  pair  has  a  different 
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Fig.  58.  Part  D  input  diodes.  The  chip  utilizes  dielectric 
isolation.  Each  division  is  8.25  pm. 
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arrangement  within  the  n-region.  Each  diode  is  formed' from  a  p-type  diffusion 
with  an  n+  diffusion  serving  as  the  cathode.  A  separate  n+  diffusion  over¬ 
laps  part  of  the  p  diffusion  and  the  n-region.  Aluminum  deposited  on  top  of 
this  serves  as  the  anode  contact.  The  base  sheet  resistance  is'  110  ohms  per 
square  and  the  breakdown  voltage  is  7.2  volts.  Once  again  it  is  really  the 
breakdown  of  an  n+-p»n  structure  operated  in  he  inverse  mode.  A  discontinu¬ 
ity  in  the  breakdown  characteristic  occurs  at  10  ma. 

Second  breakdown  waveforms  are  shown  in  Fig.  59.  They  are  dis¬ 
tinctive  because  of  the  relatively  small  voltage  drop  during  second  break¬ 
down.  For  pulse  widths  less  than  1  ysec  the  diodes  were  checked  on  a  curve 
tracer  after  each  pulse  to  determine  whether  or  not  damage  had  occurred. 

For  pulse  widths  greater  than  100  ysec  the  voltage  waveform  usually  showed  a 
second  small  d? ^continuity  such  as  appears  in  part  (c)  of  the  figure. 

In  a  few  instances  the  second  breakdown  transition  was  not  accom¬ 
panied  by  electrical  or  visual  damage.  Whenever  this  occurred  a  succeeding 
pulse  of  higher  amplitude  caused  breakdown  at  a  level  comparable  to  that  of 
an  undamaged  junction. 

Visible  damage  was  either  not  apparent  or  appeared  as  a  small  mark 
in  the  emitter  region.  Strangely  enough  the  mark  sometimes  appeared  at  the 
edge  of  the  diode  rather  than  in  the  region  between  the  contacts.  No  damage 
was  ever  seen  to  the  metallization  during  single  pulse  testing.  Two  diodes 
were  stressed  in  the  repetitive  pulse  mode  while  being  monitored  with  the 
closed  circuit  TV  system.  These  were  the  only  diodes  where  well  defined 
stress  marks  could  be  seen  extending  from  contact  to  contact  and  where 
damage  was  seen  to  the  contacts  themselves, 

The  second  breakdown  power  versus  delay  time  curve  is  plotted  in 
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Figure  60.  Note  the  well  defined  t”1  dependence  -for  t  less  than  1  usee. 

There  is  little  scatter  over  the  whole  range.  In  spite  of  the  fact  that  these 
diodes  sure  the  smallest  tested  thus  far  the  power  levels  are  not  the  lowest. 

5.  Part  E 

This  circuit  is  a  dual  U  input  diode  expander  fabricated  on  a 
dielectrically  isolated  and  glassivated  chip.  The  diodes  are  emitter-base 
diodes,  arranged  in  groups  of  U,  two  each  in  a  separate  p-type  region.  An 
aluminum  metallization  run  makes  contact  to  two  adjacent  p-regions  and  also 
connects  the  n-type  region.  A  layout  of  a  diode  group  is  given  in  Fig.  6l. 

The  p-diffusion  is  about  2.6  pm  deep  and  the  n-type  tub  is  about  13  pm  deep. 
The  diodes  have  a  reverse  breakdown  voltage  of  7.2  volts  and  show  the 
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familiar  discontinuity  due  to  n  -p-n  action. 

About  30  diodes  were  tested,  all  in  the  reverse  direction  and  all 
at  room  temperat  The  only  configuration  which  was  tested  was  between  the 
emitter  diffusions  and  the  common  aluminum  contact. 

Several  typical  breakdown  waveforms  are  given  in  Fig.  62.  These 
are  quite  similar  to  the  waveforms  of  previous  circuits.  The  voltage  drop 
is  comparatively  small  in  every  case.  The  third  photograph  shows  three 
successively  higher  pulses.  The  third  pulse  causes  breakdown  and  shows  the 
small  post-breakdown  drop  in  voltage  commonly  seen  in  previous  devices. 

In  all  but  three  cases  the  second  breakdown  transition  produced 
permanent  changes  in  the  diode  characteristics.  A  systematic  investigation 
to  determine  the  nature  of  the  damage  was  not  undertaken,  but  on  two  pack¬ 
ages  which  were  opened  the  damage  was  much  like  that  of  diodes  tested 
earlier.  Such  damage  as  could  be  seen  was  minor,  consisting  of  only  a  small 
mark  in  the  emitter  region  and  no  contact  damage. 
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Fig.  62.  Breakdown  waveforms  for  the  diodes  of  Part  E. 

(a)  10  '/cm,  500  ma/cm,  100  nsec/cm.  (b)  5  v/cm,  200  ma/om, 
lusec/cm.  (c)  5  v/cm,  100  ma/cm,  'jOysec/cm. 


The  threshold  damage  versus  pulse  width  curve  appears  in  Fig.  63. 
The  power  approaches  a  constant  value  for  t  greater  than  100  Msec  as  expected 
for  a  small  device .  At  the  other  extreme,  for  t  of  the  order  of  100  nsec, 
the  expected  asymptote  is  not  yet  reached.  This  is  probably  due  to  normal 
variations  in  threshold  power  and  is  not  necessarily  significant  in  itself. 

6.  Fart  F 

The  final  circuit  tested  was  a  dielectrically  isolated  quad  core 
driver.  There  were  four  large  interdigitated  transistors  on  each  chip.  All 
terminals  were  individually  connected,  making  it  possible  to  perform  tests  on 
all  three  configurations.  Each  transistor  measured  roughtly  20  mils  by  1*0 
mils.  The  manufacturer's  specifications  are  as  follows,  with  the  measured 
values  in  parentheses : 


bvceo 

1*0  volts  (1*3) 

bvcbo 

80  volts  (120) 

bvebo 

6  volts  (7) 

50  (90) 

~D 

800  mW  at  25°C 

The  data  obtained  during  this  portion  of  the  experimental  phase  is 
somewhat  meager  for  several  reasons.  Only  7  parts  were  available.  The  large 
size  of  the  transistors  indicated  at  the  outset  that  sufficient  pulse  power 
'-as  probably  not  available  at  short  pulse  widths.  This  was  found  to  be  true 
although  some  devices  failed  at  pulse  widths  less  than  1  yoec.  Even  more 
important,  the  threshold  power  was  erratic  from  one  transistor  to  another. 
Fhen  a  failure  point  was  determined  for  one  transistor,  others  pulsed  under 

the  same  conditions  would  fail  prematurely  resulting  in  little  useable  data. 
Finclly,  when  a  failure  was  induced,  it  was  nearly  impossible  to  locate  the 
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site  of  damage  in  the  vast  expanse  of  emitter  and  hase  fingers. 


Contrary  to  earlier  testing  procedure  the  method  used  here  was  to 
select  a  pulse  amplitude  and  then  increase  the  pulse  width  until  damage 
occurred.  Sixteen  transistors  were  pulsed  to  failure  in  this  manner  in  the 
collector-base  mode.  Eight  others  were  pulsed  collector-emitter  and  four 
emitter-base . 

Breakdown  waveforms  appear  in  Fig.  6U .  The  first  picture  shows 
collector-base  breakdown  consistently  observed  for  i ilses  longer  than  10 
ysec.  The  voltage  rises  slowly  and  the  current  drops  correspondingly  with  a 
"time  constant"  of  roughly  1*0  ysec.  Second  breakdown  manifests  itself  in  a 
sharp  drop  to  a  low  voltage.  Damage  always  followed  this  transition. 

Three  transistors  on  one  chip  failed  in  less  than  200  nsec.  All 
three  waveforms  showed  a  voltage  drop  from  130  volts  to  50  volts,  then  a 
further  drop  to  20  volts. 

Collector-emitter  pulsing  yielded  a  variety  of  waveforms.  At 
short  pulse  widths  second  breakdown  was  not  apparent  from  the  waveforms  and 
the  transistors  had  to  be  checked  after  each  pulse.  At  intermediate  pulse 
widths  breakdown  sometimes  occurred  in  tvo  separate,  distinct  steps  of  equal 
magnitude.  This  is  reminescent  of  collector-emitter  breakdown  in  the 
transistors  of  Part  B. 

At  long  pulse  widths  the  waveforms  had  a  different  appearance. 

An  example  appears  in  part  (b)  of  Fig.  64.  Two  different  pulses  are  shown. 
The  voltage  in  each  case  rises  to  a  peak,  decreases  slowly,  then  drops 
abruptly.  More  important,  both  pulses  were  applied  to  the  same  transistor. 
This  means  that  damage  was  not  detected  on  the  curve  tracer  after  the  first 
pulse.  In  the  collector-emitter  connection  it  was  often  possible  to  pulse 
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Fig.  6!*.  Breakdown  waveforms  for  core  driver  transistors. 

(a)  Collector-base  breakdown,  50  v/cm,  200  raa/cm,  10  usee /cm. 

(b)  Collector-emitter.  20  v/cm,  200  ma/cm,  20  usec/cm. 

(c)  Emitter-base.  5  v/cm,  2  A/cm,  1  usec/cm. 
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the  transistor  more  than  once  into  second  breakdown  before  degradation  was 
apparent.  Conceivably,  a  small  damage  site  might  go  undetected  because  of 
the  large  size  and  current  carrying  capability.  This  is  not  likely  though 
because  the  threshold  on  subsequent  pulses  was  in  accord  with  an  undamaged 
Junction. 

Part  (c)  of  Fig.  64  is  the  waveform  of  emitter-base  breakdown.  All 
four  transistors  on  one  die  showed  the  same  type  of  waveform  for  a  total 
pulse  width  in  the  range  from  4  to  80  usee .  Damage  was  observed  in  each  case 
following  the  appearance  of  this  unusual  waveform. 

Threshold  power  points  are  plotted  in  Fig.  65.  The  data  points 
are  so  dispersed  as  to  make  it  impossible  to  draw  meaningful  curves.  Power 
levels  are  lowest  for  the  collector-emitter  points,  contrary  to  the  transis¬ 
tor  results  of  device'  B.  The  only  conclusion  which  can  be  drawn  from  the 
figure  is  that  large  area  components  have  a  wide  variation  in  breakdown 
thresholds . 
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PART 


Threshold  damage  points  for  the  ^ore 
driver  transistors  of  Part  F. 


SECTION  V 


DISCUSSION 

Since  current  constriction,  filamentation,  and  second  'breakdown 
occur  at  temperatures  where  either  the  material  becomes  intrinsic  or  the 
thermally  generated  currents  quench  the  avalanche  breakdown,  the  temperature 
rise  of  the  silicon  is  of  major  importance  in  determining  the  outcome  of  a 
particular  experiment.  If  the  power  distribution  is  known  then  the  tempera¬ 
ture  rise  throughout  the  chip  can  be  calculated  by  techniques  such  as  those 
described  in  Section  II.  Unfortunately,  a  change  in  temperature  causes  a 
redistribution  of  the  dissipated  power  and  an  exact  analysis  becomes  very 
complex. 

In  a  reverse  biased  diode  the  greatest  power  density  occurs  at 
the  depletion  layer.  Some  of  the  pertinent  questions  relating  to  this  por¬ 
tion  of  the  total  power  are:  how  wide  is  the  depletion  layer;  ho,;  is  the 
voltage  related  to  doping  levels  and  junction  curvature;  how  does  the  volt¬ 
age  change  with  temperature ;  and  how  is  the  voltage  affected  by  high  current 
densities? 

For  emitter  Junctions  most  of  these  questions  can  be  readily 
answered.  The  high  doping  levels  result  in  a  low  breakdown  voltage,  typica¬ 
lly  between  6  and  7  volts.  Calculations  not  included  in  this  report  indicate 
that  the  depletion  layer  in  such  cases  is  very  narrow  —  of  the  order  of 
0.1  pm.  This  thickness  is  much  smaller  than  the  radius  of  curvature  of  a 
typical  emitter  diffusion;  hence,  the  breakdown  voltage  is  not  limited  by 
Junction  curvature.  Avalanche  breakdown  usually  occurs  at  the  surface 
because  of  the  heavier  base  doping  at  the  surface.  This  localization  of 
breakdown  is  further  enhanced  by  lateral  IR  drops  in  the  base  region  which 


157 


produce  an  effect  like  current  crowding.  The  upshot  of  all  this  is  that  the 
current  flows  across  the  Junction  at  the  edge  of  the  emitter  closest  to  the 
base  contact. 

Temperature  variations  of  the  breakdown  voltage  should  be  nearly 
zero  in  emitter-base  Junctions,  unlike  high  voltage  Junctions  where  the 
change  in  ionization  coefficients  produces  a  positive  temperature  coefficient. 
Consequently,  temperature  changes  would  be  expected  to  have  little  effect  in 
themselves  on  the  power  dissipated  at  the  Junction.  This  was  not  found  to 
be  entirely  true  as  will  be  explained  below. 

Changes  in  Junction  voltage  and  depletion  layer  width  can  occur  at 
high  current  densities  when  the  carrier  densities  in  the  space  charge  layer 
become  comparable  to  the  densities  of  donors  and  acceptors.  The  space  charge 
resistance  which  arises  under  these  circumstances  can  be  either  positive  or 
negative,  depending  on  the  impurity  distribution.  For  an  abrupt  Junction  the 
space  charge  resistance  is  positive,  but  for  a  p-i-n  diode  it  is  negative. 

The  remainder  of  the  applied  power  is  dissipated  in  the  diode  bulk 
resistance,  contact  resistance,  lead  wire  resistance,  etc.  If  the  geometry 
and  doping  levels  are  known,  the  parasitic  resistance  can  be  readily 
estimated. 

When  the  above  considerations  are  directed  toward  an  analysis  of 
the  diodes  of  Device  A  it  becomes  apparent  that  the  greatest  power  density 
does  indeed  occur  at  the  Junction  but  that  the  largest  portion  of  the  total 
power  i  -  dissipated  in  the  parasitic  base  resistance.  The  series  resistance 
calculated  from  the  sheet  resistivity  is  about  25  ohms.  Voltage  and  current 
were  read  from  the  leading  edge  of  pulse  waveforms  for  voltages  ranging  up 

to  30  volts  and  current  up  to  1  amp.  For  diodes  with  a  10  urn  base  region 
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these  points  fall  on  a  straight  line  with  an  incremental  resistance  between 
25  and  30  ohms.  Furthermore,  the  resistance  depends  on  the  junction  to  base 
contact  distance.  For  a  diode  with  a  base  width  of  l4  um  the  incremental 
resistance  is  between  35  and  1+0  ohms,  decreasing  tc  15-20  ohms  for  a  5  Um 
base  width. 

At  low  current  densities  the  change  in  breakdown  voltage  between 
77°K  and  U73°K  is  2  volts.  This  is  somewhat  larger  than  expected  based  on 
TVagi's  measurements.32  At  a  diode  current  of  500  ma  the  difference  in 
terminal  voltage  increases  to  k  volts.  If  we  attribute  2  volts  to  the 
change  in  breakdown  voltage  the  additional  2  volts  would  have  to  be  assigned 
to  changes  in  mobility  in  the  base  region,  a  value  which  is  not  unrealistic 
at  the  high  base  doping  levels. 

The  space  charge  resistance  appears  to  be  negligible  in  these 
diodes.  First  of  all  the  voltage-current  plot  is  a  straight  line  over  the 
entire  current  range.  Secondly,  the  Junction  is  a  linear  graded  type  with 
J'.\e  electric  field  a  parabolic  function  of  distance.  This  situation  is 
somewhere  between  that  of  an  abrupt  junction  and  a  p-i-n  junction;  hence, 
a  value  close  to  zero  is  not  unrealistic.  Finally,  the  current  density  at 
1  amp  is  an  order  of  magnitude  too  low  to  affect  the  space  charge  at  the 
surface  of  the  base  region.  The  high  doping  level  in  effect  pins  the  edge 
of  the  depletion  layer  and  prevents  it  from  expanding  into  the  base  region. 

In  summary,  it  appears  that  the  diode  resistance  in  the  avalanche 
breakdown  region  is  accounted  for  mostly  by  the  series  base  resistance.  The 
"thermal  resistance"  is  much  smaller  and  is  divided  equally  between  changes 
in  ionization  coefficients  and  changes  in  mobility.  Space  charge  resistance 
can  be  neglected. 

We  are  now  in  a  posicion  to  describe  quantitatively  the  power 
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density  throughout  the  diode.  Consider  a  reverse  biased  diode  at  room 
temperature  carrying  500  ma.  About  12  volts  appears  across  the  base  region, 
6.5  volts  across  the  depletion  layer  roughly  0.1  van  wide,  and  about  0.5  volts 
across  the  emitter  material.  If  the  power  is  dissipated  in  a  2  ym  layer  the 
power  densities  in  the  emitter,  depletion  layer,  and  base  regions  are 
3  x  10“4,  3  x  10“ A,  and  6  x  10“3  vatts/ym3,  respectively. 

The  current,  power,  and  temperature  distribution  arc  illustrated  in 
Figure  66.  As  the  temperature  rises  the  power  distribution  changes  somewhat 
due  to  the  thermal  resistance.  The  changes  also  depend  on  the  pulse  gen¬ 
erator  characteristics .  Eventually  the  temperature  reaches  a  critical  point 
where  a  constriction  occurs  either  at  the  Junction  or  in  the  base  material. 

At  a  still  later  time  second  breakdown  occurs  (at  about  600  nsec  at  these 
power  levels). 

Estimates  of  the  Junction  temperature  when  constriction  occurs  are 
difficult  to  make  for  diodes  in  the  reverse  bias  mode.  Temperatures  based 
on  solutions  to  the  linear  heat  equation  are  subject  to  the  limitations 
pointed  out  in  Section  II.  The  solutions  by  numerical  techniques  in  that 
section  are  inapplicable  because  they  are  based  on  one-dimensional  problems. 
It  is  clear  from  Figure  66  that  two  and  preferab3y  three  dimensional  models 
are  required. 

Junction  temperatures  were  estimated  from  extrapolated  steady  state 
threshold  power  versus  ambient  temperature  curves.  The  results  showed  that 
second  breakdown  does  not  occur  until  a  temperature  of  about  1100°C  is 
reached.  This  is  unexpectedly  high,  but  other  evidence  presented  later  also 
attests  to  such  high  temperatures . 

When  these  same  diodes  are  forward  biased  a  completely  different 


160 


Fig.  66.  Current,  power,  and  temperature  distribution  in  Part  A 
for  the  reverse  bias  case.  The  power  density  distribution  is 
calculated  for  a  current  of  500  roa. 


situation  holds .  In  the  first  place  the  voltage  across  the  function  can  be 
ignored.  Secondly,  conduction  is  not  restricted  solely  to  the  active  edge 
of  the  Junction  but  includes  much  more  Junction  area.  Thirdly,  the  incremen¬ 
tal  diode  resistance  at  currents  greater  than  a  few  hundred  milliamps  is 
only  about  4.5  ohms,  indicating  a  high  degree  of  conductivity  modulation  by 
injected  crrriers.  All  of  this  means  a  relatively  uniform  distribution  of 
power  over  a  fairly  broad  area  of  the  diode,  and  that  the  computer  solutions 
in  Section  II  are  valid  over  short  time  intervals  where  lateral  heat  flow  can 
be  ignored. 

Figure  29  shows  that  second  breakdown  occurs  in  100  nsec  at  a 
pulse  power  of  180  watts  in  the  forward  bias  condition.  If  the  forward  bias 
curve  is  replotted  to  indicate  threshold  power  for  constriction  formation 
(assumed  to  occur  when  the  pulse  voltage  reaches  its  maximum  value)  the  power 
level  is  150  watts  at  100  nsec.  Taking  the  power  dissipating  region  as  ex¬ 
tending  over  most  of  the  emitter  and  base  regions  between  centerlines  of 
the  contacts,  and  assuming  a  2  pm  depth  the  power  density  is  3.8  x  106  watts/ 
cm2-pm.  Figure  17  shows  that  the  temperature  rises  to  about  1300°C  when  the 
constriction  forms. 

The  threshold  curve  for  constriction  formation  should  coincide 
with  one  of  the  curves  in  Figures  17-19  at  short  pulse  lengths.  The  best  fit 
occurred  with  the  1200°C  curve  in  Figure  18  for  t  up  to  2  ysec .  The  confi¬ 
dence  in  snch  curve  fitting  is  limited  by  the  small  amount  of  data  for  the 
forward  biased  direction,  but  the  result  reinforces  the  contention  that  the 
peak  temperatures  are  very  high  at  the  onset  of  current  constriction. 

The  temperature  rise  in  the  diodes  cf  Part  B,  the  special  test 
device,  can  be  estimated  by  the  same  techniques.  In  the  reverse  bias  mode 
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the  temperature  reaches  about  1000°C  before  a  constriction  forms.  In  the 
forward  bias  direction  the  critical  temperature  is  about  1100°C. 

The  power  density  distribution  in  the  transistors  and  other  three 
layer  structures  is  complicated  by  the  interaction  between  junctions,  and  the 
multiple  conduction  paths.  Since  the  distribution  is  non-uniform  the  tempera¬ 
ture  rise  cannot  be  determined  from  the  curves  of  Section  II.  A  tabulation 
of  threshold  power  and  po(>  .  per  unit  area  allows  for  some  interesting  com¬ 
parisons  and  is  presented  in  Table  III. 

The  numbers  in  the  Device  column  give  the  pin  connections  for 
different  configurations  or  different  geometries.  The  areas  are  calculated 
in  most  cases  from  the  distance  between  contact  centerlines  and  a  distance 
slightly  greater  than  the  length  of  the  active  edge  of  the  Junction.  Some 
question  can  be  rightfully  raised  about  the  significance  of  the  area  calcu¬ 
lated  in  this  manner.  The  length  of  the  active  periphery  is  more  important 
in  some  cases.  On  the  other  hand  base  width  has  been  shown  to  be  an  important 
factor  as  has  junction  depth.  When  two  junctions  interact  the  concept  of  an 
active  perimeter  diminishes  in  importance.  For  these  reasons  the  device 
area  was  chosen  as  the  normalizing  factor  even  though  it  can  be  misleading  in 
some  cases.  Threshold  power  was  taken  directly  from  the  experimental  curves, 
extrapolated  to  steady-state  conditions  and  to  0.1  ysec  where  necessary. 

With  due  regard  for  the  vagaries  introduced  by  the  device  area  it 
is  still  possible  to  discern  trends  in  the  figures  of  Table  III,  In  the 
first  place  the  largest  devices  tend  bo  have  the  lowest  power/area  ratio  even 
though  the  total  power  may  be  greater.  This  is  expected  in  a  steady-state 
situation  because  the  thermal  spreading  resistance  is  a  nonlinear  function  of 
area.  However,  even  at  0.1  ysec  the  disparity  has  not  diminished.  It  is 
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difficiit  to  see  why  this  should  depend  on  device  area,  hence  can  only  oe 
attributed  to  differences  in  device  "hardness . " 

T;.  '  other  conclusion  derived  from  Table  III  is  that  glassivated 
end  dielectrically  isolated  devices  have  higher  threshold  power/area. 

Whether  this  -virtue  is  a  generic  property  of  these  two  technologies  cannot  be 
firmly  ostrblished  on  the  basis  of  these  limited  results.  The  uniformity  and 
consistency  of  the  data  in  this  regard  are  too  strong  to  be  ignored,  but  argu¬ 
ments  presented  below  suggest  that  the  cause  is  due  to  construction  features 
and  not  to  advantages  inherent  in  any  one  technology. 

The  question  remains,  why  is  there  a  twenty- fold  difference  in 
power/nrea  among  the  circuits  tested?  Why  is  one  circuit  harder  than 
mother  ?  The  answer  is  contained  in  the  two  general  concepts  of  power  dis- 
ribvtion  and  critical  temperature. 

The  more  uniform  the  power  density,  other  things  being  equal,  the 
greater  the  power  handling  capability.  This  is  the  main  reason  that  a  diode 
requires  more  power  in  the  forward  direction  to  reach  second  breakdown.  In 
the  reverse  direction  the  high  power  density  at  the  junction  produces  a  local 
peak  in  the  temperature  distribution  (Figure  66).  If  the  peak  were  ironed  out 
it  vo’u.--  take  longer  to  reach  the  critical  temperature. 

Thih  is  also  the  reason  that  a  low  voltage  junction  is  harder  than 
a  high  voltage  junction.  Note  the  low  threshold  power  per  unit  area  for 
collector-base  junctions.  Again,  the  power  density  is  much  greater  in  the 
depletion  layer  of  the  high  'oltagc  junction.  For  a  given  device  area  any 
modification  which  smooths  out  the  power  distribution  will  lead  to  greater 
powe>  har.dlint.  capability.  Practices  which  increase  the  parasitic  resis¬ 
tance,  decrease  breakdown  voltages,  increase  active  periphery,  increase 
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Junction  depths,  increase  contact  areas,  etc.,  will  all  lead  to  increased 
threshold  power. 

The  second  important  concept  is  that  of  a  critical  temperature. 

The  device  which  has  the  highest  critical  temperature,  other  things  being 
equal,  will  have  the  highest  threshold  power.  Any  fabrication  technique 
which  increases  the  critical  temperature  will  increase  the  threshold  power. 

In  this  category  should  also  be  included  heat  sinking  improvements  wnich 
restrain  the  rise  in  temperature . 

Within  this  framework  how  can  we  explain  the  relative  hardness  of 
certain  circuits?  All  of  the  circuits  have  roughly  the  same  breakdown  volt¬ 
ages,  sheet  resistivities,  junction  depths,  etc.  The  only  distinguishing 
feature  of  the  hard  circuits  is  the  possibility  of  multiple  current  paths,  in 
particular  conduction  via  collector  material  or  a  buried  layer  as  well  as 
along  the  surface.  The  low  threshold  circuits.  Parts  A  and  B,  have  generally 
only  a  surface  conduction  path.  Device  C,  C-B  junction,  has  multiple  paths. 
The  low  threshold  can  be  explained  however  by  the  high  breakdown  voltage  and 
hence  high  power  density.  The  collector-emitter  connection  can  also  conduct 
along  multiple  paths.  Significantly,  this  connection  can  be  classified  as 
hard  for  short  pulses.  Thus  there  is  a  certain  correlation  between  hardness 
and  current  flow  along  subsurface  paths. 

What  about  the  effects  of  heat  sinking  by  glass  and  oxide  layers? 
The  influence  of  surface  layers  of  thermal  oxide  or  deposited  glass  does  not 
have  a  significant  effect  on  the  peak  temperature.  This  was  demonstrated  in 
Section  II  and  has  been  pointed  out  earlier. 12  The  oxide  isolation  in 
dielectrically  isolated  circuits  has  been  shown  to  have  a  negligible  effect 
on  peak  temperature  for  pulses  up  to  10  used2  The  present  results  indicate 


166 


an  insignificant  effect  for  up  to  10  msec.  The  thermal  resistance  of  the 
oxide  layer  is  relatively  small  because  of  the  large  area  to  thickness  ratio 
furthermore,  its  effect  would  be  to  decrease  rather  than  increase  the  thresh 
old  power. 

Finally,  one  should  consider  the  possibility  that  the  critical 
temperature  of  the  glassivated  and  dielectrically  isolated  circuits  is  sig¬ 
nificantly  higher  for  some  reason  associated  with  the  fabrication  process. 
The  critical  temperature  of  Parts  A  and  B  was  estimated  to  be  in  the  range 
1000-1200°C.  The  maximum  increase  in  critical  temperature  could  only  be 
200-l»00°C.  The  threshold  curves  for  Device  A,  Figure  29,  show  that  this  in¬ 
crease  only  accounts  for  a  threefold  increase  in  power  levels,  even  in  a 
temperature  range  where  the  thermal  conductivity  is  changing  rapidly. 

Thus  it  appears  that  differences  in  circuit  performance  are  due 
primarily  to  differences  in  component  area  and  differences  in  the  uniformity 
of  the  power  distribution. 
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SECTION  VI 


CONCLUSION 

This  research  program  has  resulted  in  a  number  of  interesting 
observations  which  are  summarized  in  this  section.  Some  of  these  are  new, 
others  corroborate  findings  of  other  workers.  Some  are  based  on  firm  infor¬ 
mation,  others  must  be  regarded  as  tentative.  They  are  listed  here  in  sum¬ 
mary  form  without  elaboration. 

A.  Heat  Conduction 

1.  The  thermal  conductivity  of  silicon  varies  so  drastically  with 
temperature  that  the  assumption  of  a  constant  value  leads  to  large  errors  in 
heat  flow  calculations  unless  the  temperature  range  is  small. 

2.  Surface  layers  of  thermal  oxide  or  deposited  glass  have  little  heat 
sinking  capability.  Furthermore  layers  thicker  than  1  ym  prevent  heat  sink¬ 
ing  by  deposited  metal  films  for  pulse  lengths  of  1  usee  or  less. 

3.  Plated  heat  sinks  adjacent  to  heat  producing  junctions  have  little 
effect  for  pulses  less  than  10  nsec. 

B.  Second  Breakdown  in  Diodes 

1.  The  second  breakdown  model  of  Budenstein  et  al.,  applied  to  the  in¬ 
tegrated  circuits  tested  here.  In  the  reverse  direction  a  current  constric¬ 
tion  forms  at  the  junction,  followed  by  growth  of  a  current  filament  into  the 
high  resistivity  material.  When  the  filament  reaches  the  base  contact  second 
breakdown  occurs.  A  melted  column  is  produced  and  permanent  damage  occurs. 

2.  In  the  forward  bias  direction  the  filament  forms  in  the  high  resis¬ 
tivity  region  and  grows  to  entirely  span  the  region.  When  this  is  completed 


second  breakdown  occurs. 


3.  Temperatures  above  1000°C  are  reached  before  constriction  occurs. 

4.  Variations  in  threshold  power  for  a  given  part  are  caused  by  vari¬ 
ations  in  the  series  bulk  resistance. 

5.  The  threshold  level  can  be  changed  by  simple  changes  in  construction. 
If  the  base  and  collector  regions  are  shorted  together  to  form  the  anode 
some  of  the  current  flows  along  a  subsurface  path.  The  broader  distribution 
of  power  produces  a  harder  junction. 

C.  Second  Breakdown  in  Transistors 

1.  Emitter-base  breakdown  is  similar  to  ordinary  p-n  junction  breakdown. 

2.  In  collector-base  breakdown  a  dual  path  for  current  flow  exists.  At 
short  pulse  lengths  failure  occurs  along  the  surface  with  degradation  of  the 
emitter-base  Junction.  At  long  pulse  lengths  damage  occurs  beneath  the  sur¬ 
face  and  the  collector  base  junction  is  degraded. 

3.  The  threshold  power  per  unit  is  smallest  for  collector-base  pulsing 
because  of  the  high  power  density  at  the  junction. 

4.  Collector-emitter  breakdown  occurs  in  two  stages.  The  first  stage 
is  nondestructive  and  takes  place  when  a  current  constriction  elongates  from 
the  collector  junction  to  the  collector  contact.  In  the  second  stage  the 
filament  grows  through  the  base  region  and  to  the  emitter  contact.  Perma¬ 
nent  damage  then  results. 

5.  The  collector-emitter  configuration  is  the  hardest,  both  in  terms  of 
total  power  and  power  per  unit  area. 

D .  Areas  for  Further  Study 

1.  A  three-dimensional  analysis  of  the  heat  flow  problem  would  be  ex¬ 
tremely  valuable  in  analyzing  second  breakdown  phenomenon.  The  lack  of 
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definitive  temperature  distributions  in  reverse  biased  diodes  was  a  handi¬ 
cap  in  this  study. 

2.  Further  testing  and  analysis  of  transistors  and  certain  diode  struc¬ 
tures  is  needed  to  verify  details  of  the  models  proposed  here. 
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APPENDIX 


The  computer  program  described  here  is  written  in  Fortran  for  use 
on  a  Honeywell  635  time-sharing  system.  Most  of  the  calculations  in  Sec-*- 
tion  II  were  performed  with  this  program. 

As  shown  here  the  program  solves  transient  heat  flow  problems  con¬ 
taining  up  to  20  nodes  of  arbitrary  sizes  consisting  of  silicon,  silicon 
dioxide,  or  aluminum  in  any  order.  The  initial  temperature  can  have  any 
value,  but  must  be  uniform.  Power  cr.n  be  supplied  in  any  region  but  cannot 
be  a  function  of  time.  The  last  node  is  fixed  at  the  ambient  temperature. 

To  use  the  program  the  operator  divides  the  device  into  regions 
containing  any  number  of  nodes  but  only  one  kind  of  material.  For  each 
region  the  dimension  of  the  region  in  microns,  the  number  of  nodes,  the 
power  density  in  watts/cm2-ym,  and  the  kind  of  material  must  be  specified 
(1  for  aluminum,  2  for  silicon  dioxide,  3  for  silicon).  Also  the  ambient 
temperature,  time  interval  between  iterations,  final  time  for  the  problem, 
and  the  number  of  iterations  between  printout  must  be  typed  in.  The  computer 
then  lists  the  nodes,  their  location,  the  time  and  the  node  temperatures . 

Statement  numbers  150-470  read  In  sll  the  input  data.  Statements 
480-640  calculate  the  size  and  location  of  each  node.  A  portion  of  the 
printout  is  handled  by  statements  650-720. 

Statements  730-1770  are  a  tabulation  of  the  thermal  conductivity 
and  specific  heat  of  aluminum,  silicon  dioxide,  and  silicon  from  0  to  1700°K 
at  50  degree  intervals,  A  straight  line  interpolation  is  used  at  inter¬ 
mediate  temperatures . 

The  heat  flow  equations  are  set  up  in  statements  1800-2230.  New 
coefficients  are  calculated  at  the  beginning  of  each  iteration.  The 


equations  are  solved  in  statements  221*0-2330. 

An  example  of  a  problem  solution  is  included.  The  device  is 
divided  into  11  regions.  The  first  region  is  an  aluminum  layer  1  pm  thick 
containing  1  node.  The  second  region  is  silicon  dioxide,  0.5  urn  thick,  with 
1  node.  The  third  region  is  silicon,  3  ym  thick,  with  3  nodes.  The  fourth 
region  represents  the  depletion  layer  of  a  junction  and  contains  the  source 
of  heat.  This  region  is  1  ym  thick,  has  1  node,  and  has  a  power  density  of 
6  x  10s  watts/cm2-ym.  This  is  equivalent  to  about  50  watts  in  a  1  ym  layer 
of  a  device  measuring  70  x  100  ym. 

Regions  5  through  10  are  all  silicon  with  gradually  increasing 
node  size  to  provide  a  smooth  transition  to  a  relatively  thick  substrate. 

For  region  11  a  group  of  zeros  is  typed  in  to  terminate  the  input  data. 

The  ambient  temperature  is  27°C.  The  time  interval  is  10  nsec 
with  a  final  time  of  2.1  ysec. 
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10******THAP  -  TRANSIENT  HEAT  ANALYSIS  PROGRAM.  MAXIMUM  NUMBER  OF 
20******NODES  =  SO*. NODE  1  HAS  ZERO  HEAT  TRANSFER  OUT  OF  THE  SURFACE. 
3f)******LAST  NODE  CONSTRAINED  TO  AMBIENT  TEMPERATURE. 

40  DIMENSION  S ILICONK(35  >*  S ILICONC(35 )*  SI02K(35>* 

504  SI02C  (35  >  *  ALK(35 )*  ALC(35)*  X(20)*  DX(20>*  0(20)* 

604  T  (  2'0 )  *  NODES  (20)*  XRC20 )  *  OR<20)*  NR<20)*  MR(20)* 

704  MAT (20 ) *  R(20)*  RHOCDD(20)*  A(20*21  )*  S(20>* 

804  F(20 ) 

90  PRINT  99 

100  99  FORMAT (  ////"  ####♦  THAP  -  TRANSIENT 

1104  HEAT  ANALYSIS  PROGRAM  #####") 

120  PRINT  100 

130  100  FOR  MAT  (  ///••  DIMENSIONS*  NODES*  HEAT  IN 

14 04 PUT*  MATERIAL"/) 

150******TYPE  IN  DIMENSIONS  OF  EACH  REGION  IN  MICRONS*  NUMBER  OF 
160******NODES  IN  THAT  REGION*  POWER  DENSITY  IN  WATTS  /C  M2  /MICRON* 
170**%***AND  TYPE  OF  MATERIALS  ALUMINUMS*  SILICON  DI0XIDS=2* 
180******AND  S I LIC0N=3 • 

190  DO  102  NI  *  1*20 

200  PRINT  103*  NI 

210  103  FORMAT (  "  REGION  "12) 

220  READS  XR(NI )*NR(NI )*QR(NI )*MR(NI ) 

230  Nil  =  NI  -1 

240  N  =  N  +.NRCNI) 

250  IF(XR(NI))  102*104*102 

260  102  CONTINUE 

270  104  CONTINUE 

280  M  =  N  +  1 

290  L  *  N  -  1 

300  PRINT  110 

310  110  FORMAT(  //"  AMBIENT  TEMPERATURE") 

320******TY?S  IN  AMBIENT  TEMPERATURE  IN  DEGREES  CENTIGRADE 

330  R? ADs  TA 

340  DO  1 1 1  IN  =  1*20 

350  111  T  ( I N )  «  TA 

360  PRINT  150 

370  150  FORMAT (  //"  TIMS  INTERVAL" ) 

380******TYPS  IN  TIME  INTERVAL! ALSO  INITIAL  TIME)  IN  SECONDS* 

390******FINAL  TIME  IN  SECONDS*  AND  NUMBER  OF  TIME  STEPS  BEFORE 

400******SACH  PRINTOUT  OCCURS 

410  READ:  DT 

420  PRINT  151 

430  151  FOR  MAT (  "  FINAL  TIME") 

440  READS  FT 

450  PRINT  152 

460  152  FORMAT!  "  STEPS  BETWEEN  PRINTOUT") 

470  READS  NS 

480******CALCULATE  GEOMETRIC  FACTORS 
490  IB  =  : 

500  IC  =  NR ( 1 > 
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500 

IC  =  NR( 1 ) 

510 

DO  125  IA  r  1 *NI ! 

520 

XX  =  XR? IA)/NKC IA) 

530 

QQ  =  QR(IA) 

540 

MM  =  MR?IA) 

550 

DO  126  I  =  IB* IC 

550 

MAT<I>  =  MM 

570 

DX< I >  a  XX 

580 

QC I >  =  QQ 

590 

126 

NODES? I>  =  I 

600 

IB  =  IC+1 

610 

125 

IC  =  IC+NR< I A+l  ) 

620 

XC1)  =  0.5  *DX? 1 > 

630 

DO  127  I  =  2*N 

640 

127 

XCI>  =  X?I-1 )+0.5*(DXCI-l)+DX(I)) 

650 

PRINT  170 

660 

170 

FORMAT?  ///" 

SOLUTION  ") 

670 

PRINT  171*  ? NODES? I )*  I  =  1*N> 

680 

171 

FORMAT?  //"  NODES  "  817/1 IX  8I7/12X 

417) 

690 

PRINT  172*  ?X?I)*  I  =  l *  N) 

700 

172 

FORMAT?  "  POSITION  "8F7.2/13X  8F7.2/14X  4F7.8) 

710 

PRINT  173 

780 

173 

FORMAT?  /"  TIME 

NODE  TEMPERATURES" > 

730******TA8LE  CF  THERMAL  PROPERTIES*  K  IS  THERMAL  CONDUCTIVITY 
740  ******  IN  WATTS /C  M/DE  GRSE  AND  C  IS  SPECIFIC  HEAT  IN  CALORIES 
750******/GRAM/DSGREE 
760  ALK?  1  )  =50* 

770  ALK? 2)  =  8* 

780  ALK?3)  =2*9 

790  ALK? 4)  =2*5 

800  ALK<5  >  =  2*4 

810  ALK<6>  =  2*39 

820  ALK?  7)  =  2.38 

830  ALK?8>  =  2.38 

849  ALK<9>  =  2.38 

850  ALK<  10)  =  2.38 

860  ALK< 1 1 )  =2.38 

870  DO  10  IND  =  12*35 

880  10  ALKCIND)  =  ALKUND-l)  -  0.03 

890  ALCCi)  =  0. 

900  ALC?  2 )  =  .04 

910  ALC  <  3  >  =  .12 

920  ALC  ?  4 )  =  .165 

930  ALC  (5)  =  .193 

943  ALC  <  6>  =  .207 

950  ALC<7>  =  .217 

960  ALC  ?  8  )  =  .225 

970  ALC  ?  9  )  *  .233 

980  DO  12  IND  =  10*35 

990  1  2  ALC?IND>  =  ALC(IND-l)  +  .005 

1006  SI02  K<  1  )  =  0. 

1010  DO  14  IND  =  2*6 
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1080  14  SI02K<  iND)  =  S I02KC  IND-1  > 

1030  DO  16  IND  *  7*35 

1040  16  SI02KCIND)  =  SI.02KC  IND-1  > 

1050  SI02CC1  )  «  0. 

1050  SI02C  (2  )  =  »034 

1070  SI02CC3)  *  *07 

1080  SI02CM)  =  *097 

1090  SI02f,  <5)  *  .125 

1100  SI02CC6)  =  .15 

1110  SI02C<7)  =  .175 

1120  S I02C  <8 >  »  .19 

1130  SI02CC9)  =  *205 

1140  S 1 02  C  (  1 0 )  «  .82 


ii5o; 

SI02CC11)  = 

.231 

1163 

SI02C  <  1 2 )  a 

.242 

1170 

SI02CC13)  = 

.25 

1180 

SI02CC14)  » 

.256 

1190 

SI02CU5)  « 

.262 

1200 

SI02CC16)  « 

.267 

1210 

SX02CC17)  a 

.271 

1220 

SI02CC18)  a 

.275 

1230 

SI02CC19)  a 

.279 

1240 

SIO2CC20)  » 

.282 

1250 

SI02CC21  )  » 

.285 

1260 

S I02C  <22  )  = 

.288 

1270 

SI02CC23  >  a 

.891 

1280 

DO  18  IND  « 

24.35 

1890 

18  SI02CCIND)  = 

SI02C 

1300 

SILICONKC  1) 

a  0. 

1310 

SILIC0NK<2) 

=  24. 

1320 

S ILICONKC3  > 

a  8.8 

1330 

SILIC0NKC4) 

a  3.9 

1340 

S ILIC0NKC5  ) 

a  2  «6 

1350 

SILIC0NKC6) 

=  1.9 

1360 

SILIC0NKC7) 

=  1  *5 

1370 

SI  LICQNKC  8  > 

a  1.2 

1580 

SILIC0NKC9) 

=  1  .0 

1390 

SILICONKC10) 

=  .85 

1400 

SILICONKCl 1  > 

=  .75 

1410 

SILIC0NKC12) 

=  .67 

1420 

SILICONKC 13  > 

=  .61 

1430 

SILICONKC  14) 

=  .55 

1440 

SILICONKC  15  > 

=  .50 

1450 

SILIC0NKC16) 

=  .45 

1460 

SILICONKC17) 

=  .42 

1470 

SILIC0NKC18) 

a  .38 

1480 

SILICONKC  19) 

=  .35 

’490 

SILICONKC20) 

=  .33 

1500 

S  ILIC0NK(2 1  > 

=  .31 

1510 

SILICONKC  22) 

=  .30 

1520 

SILIC0NKC23) 

=  .29 

+  .00272 
+  -.0005  82 


+  .002 
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./....HWHiMt-WW*  *.«*v*m*  •***•* *■ cjy. 


15  30 

S  I L I C  ON  K  ( 24  ) 

=  .88 

1540 

S ILICONKI25 ) 

=•  .28 

15  50 

SILIC0NK<26> 

=  .27 

156  0 

SILIC0NK(87> 

=  .27 

15  70 

S ILIC0NKI28  > 

=  .27 

15  80 

SILIC0MKC29) 

=  .26 

15  90 

S I  LICONK(  30 ) 

=  .26 

160.0 

SILIC0NKC31 > 

=  .26 

1610 

S I LICONKC 38  ) 

=  .26 

1520 

SILIC0NKC33 ) 

=  .26 

1630 

S  I-L1C0MK<34  > 

=  .25 

1640 

SILIC0NK(35) 

=  .25 

165  0 

S ILICONC  Cl)  = 

0. 

166  0 

S ILICONC <2  )  = 

.03  5 

16  70 

S  ILICONC  (3  = 

.070 

1680 

S I L I C  UNC  <  4  )  = 

.110 

15  90 

3 ILIC0NC(5  )  = 

.130 

1700 

S ILICONC (6  )  = 

.155 

1710 

S ILICONC (7)  - 

.175 

1720 

S ILICONC  <  8 )  = 

.  1  85 

1730 

S ILICONC (9)  = 

.195 

1740 

S ILICONC < 10) 

=  .200 

175  0 

SILICONCI 1 1 > 

=  .205 

1760 

DO  21  KK  =  12 

.35 

1770 

21  S I LICONC ( KK) 

=  S ILICONC (KK~1 >  +  .002125 

1780 

TIME  =  0. 

1790 

IS  =  FT/DT/NS 

1800******solvs  heat  flow  squat  ions 


1810  DO  71  IT  =  l#  IS 

1820  DO  70  ITT  =  1 > NS 

1830******CALCULATE  THEHMAL  PROPERTIES  AT  THE  NODE  TEMPERATURES 
1840  DO  80  I  =  1«N 

1850  ENTER  =  (TCI)  +  2730/50.  +  1. 

1860  II  =  ENTER 

1870  AA  »  II 

1880  BB  =  ENTER  -  AA 

1890  K  «  MAT  (  I  ) 

1900  GO  TOC62>61*60)>K 

1910  60  AK  =  SILICON Kill) +BB*<SILIC ON K(II+1>-SILIC ON K(II>) 

1920  RC  =  SILIC0NCUI>+33*(SILIC0NCCII+l  > -S ILICONC ( 1 1 > > 

1930 **♦*♦>•<0 ONUS R T  SPECIFIC  HEAT  FROM  CALGRISS/GRAM/DEGRSS  TO 


1940****** JOULES /C M3 /DEGREE 
1950  RHOC  =  4.1 84*2. 3S8*R0 

1960  GO  TO  63 

1970  61  AK  =  SI02K(II  )+B9*CSI02K(1 1+1  >-SI02K(ID) 

1980  RC  =  SI02C(II)+i33*(SI02C(II+l>-SI02C(II>) 

1990  RHOC  =  4*184*2  »20*RC 

8000  GO  TO  63 

2010  62  AK  =  ALK< 1 1 >+BB*CALK< 1 1+1 >-ALKC II ) > 

2020  RC  =  ALC< II )+3B*<ALC( 1 1+1 > -ALC (I  I  )  > 

2030  RIIOC  =  4  •  1  84*2  »  707*RC 


179 


2040 

GO  TO  63 

2050 

63 

■RC  I  )  =  DXC  I  >/AK*l  .E-4 

206  0 

20 

RHOCDDCI)  =  1  .5  -4*RH0C*DXC  I  > /DT 

2070******CALCULATE  COEFFICIENTS  IN  EQUATIONS 

2030 

DO  30  I  -  1 #N 

2090 

IM  =  1-1 

2100 

IP  =  1+1 

2110 

IF  C IM)3 1 #3 1 #32 

2120 

32 

AC  I » IM)  =  2./CRC IM)+RC I  )) 

2130 

IFC IP-N)36#36#31 

2140 

36 

AC  I  #  I  )  =  -2./CRCIM)+RCI )>-2./CRCI)+RCIP>>- 

21504 

RHOC DDC I ) 

2160 

31 

AC  I #  M)  =  -DXC I  ) *Q  C  I  ) -RHOCDDC I  )  *T  C  I  ) 

2170 

IFC IP-N>34#34»35 

2180 

34 

AC  I  #  IP)  =  2./CRCD+RCIP)) 

2190 

35 

CONTINUE 

2200 

30 

CONTINUE 

2210 

AC  1  #  1  )  =  -2./CRC1 )+R<2) )-RH0CDDCl ) 

2220 

L  ■  N-l 

2230 

ACN#N)  =  -2»/CRCL)+RCN) ) -RHOCDDC N) 

2240******SOLVS  EQUATIONS  USING  RICHfMYER'S  ALGORITHM 
2250  FC1  )  =  AC  1  #M)/AC  1/1  ) 

2260  EC1)  =  -AC1#2)/AC1#1 ) 

2270  DO  50  J  n  2#N 

2280  FCJ)  =  <A< J#M)-AC J* J-1)*FCJ-1 >>/ 

22904  CAC  J#  J)+AC  J#  J-l  >*EC  J-l  )) 

2300  50  ECJ)  «  -AC  J#  J+l  ) /C  AC  J#  J) +AC  J#  J- 1  )  *E  C  J-l  >  ) 

2310  TCN-1)  =  TA+ECN-1)  +  FCN-l) 

2320  DO  5  1  J  =  2#L 

2330  51  TCN-J)  *  E C N~ J) *T C N- J+l )  +  FCN-J) 

23 A 0  TIME  *  TIME  +DT 
2350  70  CONTINUE 

2360  PRINT  9#  TIME#  CTCJJ)#  JJ  *  1#N> 

2370  9  FOR  MAT  C  S10.3#2X  8F7.1/13X  8F7.1/14X  4F7.1) 

2380  71  CONTINUE 

2390  STOP 
2400  END 


0**0*  THAP  -  TRANSIENT  HEAT  ANALYSIS  PROGRAM  00000 


DIMENSIONS#  NODES#  HEAT  INPUT#  MATERIAL 

REGION  1 
°  1»#1»0*#1 
REGION  2 
=  *5  »  1  # 0  •  #2 
REGION  3 
n  3 • #3  #0  •  #3 
REGION  4 
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»  1 »>1>6 «E5 >3 
REGION  5 
3  i •> 1 >0«>3 
R2GION  6 
B  4.>2>0.> 3 
REGION  7 
a  8»>2.»0».»3 
REGION  8 
a  !6»>2»0»#3 
H2GION  9 
3  30 ♦ #  3  » 0 •  *  3 
REGION  10 
a  60  *  *4  >0  *  *3 
REGION  11 
=  0>0*0*0 


AMBIENT  TEMPERATURE 
=  27.  1'0 


TIME  INTERVAL  . 
a  1  0  .E  -9 
FINAL  TIME 
a  2« IE  -6 

STEPS  BETWEEN  PRINTOUT 
a  50 

SOLUTION 


NODES 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

13 

14 

15 

1.6 

17 

18 

19 

20 

POSITION 

0  »5  0 

1  *25 

2 .00 

3.00 

4  *00 

5*00 

6  *00 

7.50 

9.50 

12.50 

16.50 

22.50 

30.50 

39*50 

49.50 

5  9.5  0 

72  .00 

87.00 

102  .00  1 1  7.00 

TIME 

NODE 

TEMPERATURES 

0*5  00E  ”06 

56  »6 

105.4 

183.7 

195.2 

213.6 

240.4 

196.4 

147.9 

106  .2 

68.1 

45  .0 

31  .0 

27.6 

27.1 

27.0 

27.0 

27.0 

27.0 

27. 

0  27.0 

0. 1  002  -05 

.128.8 

196.3 

296*7 

312.9 

336.9 

3  70.2 

310.0 

241.8 

179.2 

119.5 

77.3 

45  .2 

31.9 

27.9 

27.1 

27.0 

27.0 

27.0 

27. 

0  27.0 

0.150S-05 

214.7 

286.8 

393.6 

413-9 

4*3. 3 

483.9 

408*1 

323.8 

244.8 

167.3 

109.4 

63.2 

39.7 

30.3 

27-7 

27.1 

27.0 

27.0 

27. 

0  27.0 

0.200E-05 

304.0 

3  75  *5 

482  .5 

507.2 

542-1 

5  90.2 

498.9 

398.6 

305*0 

212.5 

140.7 

82.1 

49.3 

34.3 

29.0 

27.5 

27.1 

• 

t- 

co 

27. 

0  27.0 
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